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FOREWORD 


This  report  details  the  work  performed  under  the  direction  of 
Dr.  S.  V.  Shelton  at  the  School  of  Mechanical  Engineering,  Georgia 
Institute  of  Technology,  Atlanta,  Georgia,  under  Contract  No. 
F08635-70-C-0129  with  the  Air  Force  Armament  Laboratory,  Eglin  air 
Force  Base,  Florida  32542.  Lieutenant  William  S.  Bulpitt  (DLDL) 
was  program  manager  for  the  Armament  Laboratory.  This  effort  was 
conducted  during  the  period  from  15  July  1970  to  15  July  1972.  The 
first  12  months  were  devoted  to  modeling  the  heat  transfer  mechanism 
in  gun  barrels;  during  the  second  12  months  this  model  and  its 
results  were  used  to  aid  in  understanding  erosion  in  gun  barrels 
and  its  reduction  by  particle  seeding  (TiC^)  of  the  propellants. 

This  report  is  divided  into  two  parts.  Part  I  -  Analytical 
Model  of  Heat  Transfer  in  Gun  Barrels  consists  of  Section  II  through 
VI,  and  Part  II  -  Particle  Seeding  Effects  on  Gun  Barrel  Heat  Transfer 
and  Erosion  consists  of  Sections  VII  through  XI. 

This  technical  report  has  been  reviewed  and  .;s  approved. 

DALE  M.  DAVIS 

Director,  Guns  and  Rockets  Division 
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ABSTRACT 


This  study  models  tne  combustion,  solid  propellant  move¬ 
ment,  gas  dynamics,  projectile  dynamics,  transient  boundary  layer,  heat 
transfer  and  barrel  metal  temperatures  in  a  gun  barrel.  This  is  accom¬ 
plished  by  rigorous  development  of  the  coupled  partial  differential  equa¬ 
tions  and  carrying  out  a  detailed  numerical  solution  to  these  equations. 
Comparison  of  barrel  temperature  solutions  to  experimental  data  is  shown. 
Erosion  mechanisms  are  discussed  in  the  light  of  these  solutions.  This 
model  is  then  used  to  study  hypotheses  concerning  the  heat  transfer,  tem¬ 
perature,  and  erosion  effects  of  submicron  size  solid  particle  additives 
(TiO  )  to  8un  propellants.  A  mechanism  not  previously  studied  offers 
excellent  theoretical  results  in  explaining  the  reduced  erosion. 
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SECTION  I 


SUMMARY 

The  objective  of  this  study  is  to  (1)  provide  a  mathematical  model  which 
can  be  used  to  predict  the  performance  and  heat  transfer  characteristics  of 
guns,  or  devices  which  produce  high  pressure  and  temperatures  in  an  enclosed 
but  expanding  volume  by  burning  solid  propellant  and  (2)  utilize  this  model 
to  study  the  heat  transfer  and  erosion  effects  of  adding  solid  micron  size 
particles,  such  as  titanium  dioxide,  to  the  propellant. 

The  internal  ballistic  model  is  developed  by  deriving  the  time-dependent 
one-dimensional  conservation  equations,  including  the  effect  of  skin  friction 
and  wall  heat  transfer.  A  boundary  layer  analysis  is  carried  out  by  deriving 
the  boundary  layer  momentum  equation  for  a  non-steady,  developing  compressible 
flow  in  a  tube.  As  a  first  approximation,  the  profile  shape  fraction  (ratio 
of  the  displacement  thickness  and  the  momentum  thickness)  is  assumed  to  be 
a  constant.  The  conduction  heat  transfer  equations  for  the  tube  itself  are 
written  and  coupled  with  this  boundary  layer  heat  transfer  analysis.  Two 
limiting  cases  of  solids  velocity,  namely  (1)  the  same  velocity  as  the  com¬ 
bustion  gases  and  (2)  zero  velocity,  are  considered. 

Numerical  techniques  are  used  to  solve  the  above  equations  together 
with  the  equation  of  state  of  the  combustion  gas  and  to  determine  all  the 
ballistic  properties,  namely  pressure,  velocity,  gas  density,  gas  temperature, 
volume  fraction  of  solids,  and  the  boundary  layer  thickness  at  each  point 
along  the  length  of  the  tube  at  every  time  step.  The  heat  transfer  co¬ 
efficient  at  the  inside  surface  of  the  tube  is  obtained  from  Colburn's 
analogy  and  the  tube  wall  temperature  is  determined  simultaneously,  and 
coupled  with,  the  interior  ballistic  solution. 


1 


’1 


Results  are  obtained  for  a  typical  set  of  input  data  which  show  that 

the  final  projectile  velocity  and  the  time  of  travel  for  the  two  extreme 

cases  of  solids  velocity  are  quite  close.  The  total  heat  loss  to  the  tube 

wall  is  found  to  be  five  to  six  percent  of  the  input  energy.  Very  high 

2 

values  of  the  heat  transfer  coefficient  (50  kcal/m  -sec-°K  on  the  average) 
are  found,  which  produce  peak  tube  wall  surface  temperatures  of  800  to 
1000°K  near  the  projectile  standing  position. 

Excellent  agreement  between  predicted  and  experimentally  measured 
single  shot  tube  wall  temperatures  are  found.  This  agreement  is  produced 
without  the  use  of  experimentally  measured  interior  ballistics  data  such 
as  chamber  pressure  versus  time. 

Repetitive  firing  calculations  were  made  for  ten  to  fifty  repeated 
firings  with  the  tube  heat  transfer  analysis  yielding  tube  wall  temperatures. 
For  fifty  firings  over  a  three  second  period  (0.060  second/firing)  the  inside 
tube  wall  peak  temperature  reaches  1500°K.  Repeated  firing  bursts  would  pro¬ 
duce  temperatures  high  enough  to  cause  melting  of  a  very  thin  layer  of  metal. 
This  verifies  che  importance  of  heat  transfer  in  the  gun  tube  erosion  process. 

The  possible  effects  of  solid  micron  size  particles  on  the  heat  transfer 
from  the  combustion  gases  to  the  tube  wall  are  investigated.  A  review  of 
the  existing  literature  concerning  experimental  measurements  on  dusty 
gases  reveal  that  a  reduction  in  wall  friction  and  convective  heat  transfer 
is  possible  in  many  cases.  In  fact,  reductions  up  to  30  percent  have  been 
measured,  but  no  experimental  heat  transfer  or  friction  data  at  the  Reynolds 
number,  particle  loading,  and  particle  sizes  applicable  to  the  gun  tube 
problem  are  available.  If  it  is  assumed  that  the  particles  cause  a  reduc¬ 
tion  of  the  convective  heat  transfer  coefficient  of  20  percent,  the  inside 
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tube  wall  temperature  is  reduced  by  about  100°K,  which  appears  as  insufficient 
to  cause  a  significant  reduction  in  erosion. 

All  data  on  wear-reducing  gun  propellant  additives  reveal  that  a  coating 
remains  on  the  tube  wall  after  numerous  firings.  A  study  of  the  insulation 
effect  of  micron  size  particles  forming  a  packed  bed  in  the  surface  rough¬ 
ness  crevices  of  the  tube  wall  showed  that  the  tube  metal  temperatures  would 
be  reduced  by  up  to  300°K.  This  occurs  for  a  surface  roughness  of  ten 
microns.  It  is  found  that  sufficient  particle  loadings  are  contained  in 
propellants  used  in  actual  practice  to  supply  this  coating  material. 

This  insulation  mechanism  hypothesis  is  upheld  by  all  calculations  and 


is  concluded  to  be  the  most  probable  predominate  wear-reducing  mechanism 
of  wax  with  Titanium  dioxide  particles.  This  deposition  of  particles  should 
therefore  be  optimized. 
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SECTION  II 


INTRODUCTION 

Definition  of  the  Problem 

Devices  which  produce  high  pressure  in  an  enclosed  but  expanding 
volume  by  burning  combustible  mixture  of  gases  or  solid  propellant  with 
the  objective  of  performing  work  are  common  in  practice.  Internal 
ballistics  of  these  devices,  for  example  the  problem  of  the  gun,  have  been 
solved  experimentally  since  the  fourteenth  century  when  gunpowder  first 
came  into  use  [1]*.  But  surprisingly  enough,  an  analytical  solution  which 
may  be  used  to  accurately  predict  the  performance  of  such  devices  has 
not  been  found.  This  lack  of  mathematical  model  compels  a  designer  to 
choose  the  comparatively  expensive  path  of  experimentation,  although 
only  limited  information  can  be  obtained  from  these  experiments.  More¬ 
over,  a  large  number  of  experiments  have  to  be  performed  before  a  set 
of  optimum  design  parameters  can  be  determined  for  a  particular  purpose, 
and  still  the  final  result  remains  in  question  as  to  whether  a  tiuly 
optimum  condition  has  been  achieved. 

The  problem  of  internal  ballistics  requires  a  modeling  of  the 
fluid  flow  phenomena  and  heat  transfer  to  the  wall  inside  the  expanding 
volume.  For  simplicity,  throughout  this  work  we.  shall  be  restricted  to 
the  special  geometry  of  a  closed  cylindrical  tube  with  a  sliding  piston 


^Number  in  [  ]  refers  to  the  references  in  Bibliography. 
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at  one  end,as  shown  in  Figure  1.  The  combustible  mixture  is  burnt 
inside  the  enclosed  volume  whereby  the  pressure  is  increased  and  the 
piston  is  set  into  motion.  The  products  of  combustion  which  flow  down 
the  cylinder  behind  the  piston  impart  a  considerable  amount  of  its 
energy  to  the  piston  and  a  fraction  is  lost  to  the  tube  wail.  This 
cools  the  combustion  gases  and  modifies  the  pressure  and  flow  conditions. 

while  heat  transfer  has  some  effect  on  the  ballistic  properties, 
this  is  probably  more  important  with  respect  to  the  material  properties 
of  the  tube.  Since  the  combustion  gases  are  usually  at  a  temperature 
of  2000-3000°K,  after  repeated  use  of  the  device  at  high  frequency  the 
wall  temperature  of  the  tube  may  reach  a  value  high  enough  to  cause 
appreciable  wear  as  the  piston  slides  down  the  tube.  A  model  of  heat 
transfer,  which  can  be  used  to  predict  the  wall  temperature,  will  help 
a  designer  to  choose  the  optimum  design  parameters  which  will  minimize 
the  erosion  rate. 

The  purpose  of  the  present  research  is,  therefore,  to  provide  a 
working  analytical  model  which  shall  be  able  to  predict  all  the  ballistic 
properties,  namely  velocity,  pressure,  temperature  and  density  of  the 
combustion  gas  mixture  as  a  function  of  space  and  time.  The  heat  loss 
to  the  tube  wall  shall  be  considered  and  the  temperature  distribution 
at  the  wall  shall  be  determined.  This  model  will  then  allow  study  and 
optimization  of  various  parameters  without  expensive  trial  and  error 
experimentation. 


7 


Related  tfork 


Theoretical  solutions  to  the  problem  of  interior  ballistics  have 
been  attempted  since  the  days  of  Lagrange  who  in  1'93  first  tried  to 
determine  the  spatial  distribution  of  pressure,  density  and  gas  velocity 
in  the  tube  at  air  times  after  the  combustion.  The  work  available  until 
now  can  be  divided  into  two  broad  categories: 

(1)  Semierapirical  solutions  which  may  have  practical  utility  in 
the  study  of  familiar  devices. 

(2)  Exact  theories  which  attempt  to  include  the  predominate 
phenomena  up  to  a  certain  order  of  magnitude  by  formulating  a  simple 
mathematical  model  of  the  flow. 

Semiempirical  Solutions 

The  major  works  in  this  area  with  special  application  to  the 
guns  using  solid  propellant  are  described  in  references  1  and  2 
The  main  purpose  of  these  works  is  to  obtain  a  solution  which  matches 
with  the  experimental  values  of  peak  chamber  pressure  and  muzzle  velocity 
of  the  projectile.  Only  a  few  of  the  number  of  solutions  shall  be 
discussed  here. 

Isothermal  Solution.  The  solution  as  described  by  Corner  [2]  is 
based  on  the  following  assumptions: 

(1)  The  propellant  stays  in  the  chamber  burning  under  the  tube 
head  end  (breech)  pressure  and  the  rate  of  burning  is  proportional  to 
that  pressure. 

(2)  During  the  period  of  burning  of  the  propellant,  the  progres¬ 
sive  cooling  of  the  combustion  gases  due  to  the  work  done  on  the  projec¬ 
tile  can  be  approximated  by  taking  a  mean  gas  temperature  over  this  time 
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interval,  corresponding  to  an  effective  mean  force  constant  X. 

(3)  Uniform  gas  density  and  linear  velocity  distribution  in  the 
space  between  the  tube  head  end  and  piston  base. 

(4)  Resistance  to  motion  of  the  projectile  can  be  taken  into 
account  by  introducing  an  increased  effective  projectile  mass  instead 
of  actual  mass. 

(5)  The  covolume  q  (volume  correcting  term  in  the  equation  of 
state  of  the  combustion  gas)  is  equal  to  the  specific  volume  of  the 
propellant  material. 

The  expressions  for  breech  pressure  P,  projectile  velocit>  V, 
and  projectile  distance  from  breech  face  x,  are  given  as  a  function  of 
"convenient  variable"  f,  the  form  factor  0,  the  force  constant  X, 
burning  rate  6,  and  central  ballistic  parameter  M.  The  central  ballistic 
parameter  M  itself  is  a  function  of  X,  3,  initial  mass  and  web  size  of 
propellant,  effective  projectile  mass,  and  tube  diameter.  The  form 
factor  9  depends  on  the  geometrical  shape  of  the  propellant  and  the 
variable  "f"  goes  from  one  to  zero  as  the  propellant  is  burnt.  Other 
parameters,  namely  X,  M  and  3  are  chosen  following  a  trial  and  error 
procedure  until  good  agreement  is  obtained  with  the  experimental  values 
of  peak  pressure  and  muzzle  velocity.  The  solution,  however,  does  not 
take  into  account  the  heat  loss  to  the  tube  wall. 

Coppock’s  Solution  [2] .  This  is  an  extension  to  the  isothermal 
solution  described  above  with  the  following  modifications: 

(1)  Instead  of  taking  a  mean  gas  temperature  during  burning,  the 
analysis  takes  into  account  the  kinetic  energy  of  the  projectile  and  that 
of  the  gases,  assuming  that  the  combustion  gases  are  uniform  in  density 
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between  the  breech  and  the  projectile  and  that  their  velocity  at  any 
point  is  proportional  to  the  distance  from  the  breech  face.  The  total 
heat  loss  to  the  tube  wall  up  to  a  particular  instance  of  time  is 
assumed  to  be  a  certain  fraction  of  the  total  kinetic  energy  of  the 
projectile  and  the  gases  at  that  instant,  In  practice,  the  effect  of 
heat  loss  is  incorporated  in  the  energy  equation  by  a  proper  choice  of 
Y  (ratio  of  specific  heats  at  constant  pressure  and  constant  volume). 

(2)  The  gases  have  a  constant  covolume  n,  not  necessarily  equal 
to  the  specific  volume  of  the  propellant  material. 

From  the  observed  peak  pressure  it  is  possible  to  back-calculate 
the  central  ballistic  parameter  M,  and  thence  the  burning  rate  B. 

The  solution  is  superior  to  the  isothermal  solution  because  there 
is  only  one  arbitrary  parameter,  namely  the  burning  rate  B,  whose  value 
is  selected  so  that  the  peak  pressure  matches  the  experimental  data. 
Moreover,  the  model  takes  into  account  the  heat  ]oss  to  the  tube  wall, 
though  in  a  crude  fashion. 

Goldie^  Solution  [2].  The  solution  follows  Coppock's  solution 
described  above  with  the  only  modification  that  the  projectile  is  assumed 
to  be  motionless  until  a  "shot-start  pressure"  is  produced  inside  the 
chamber.  If  there  is  any  resistance  to  motion  at  later  times,  the  effect 
is  simulated  by  a  change  in  effective  shot  weight. 

Apart  from  these  solutions,  there  are  solutions  which  attempt  to 
use  a  better  relationship  between  the  burning  rate  and  the  corresponding 
pressure.  But  the  solutions  '•till  need  trial  and  error  of  one  or  more 
variables  to  match  experimental  data.  Besides,  there  is  no  guarantee 
as  to  how  good  the  solutions  will  be  when  prediction  of  performance  of 


11 


a  new  device  is  desired.  Also  no  information  regarding  the  ballistic 
properties  in  between  the  breech  face  and  the  projectile  is  available 
from  any  of  these  models.  Even  a  recent  publication  [3]  fails  to  pro¬ 
vide  such  informations. 

Exact  Theories 

As  mentioned  earlier,  Lagrange  took  the  initiative  to  solve  the 
one-dimensional  problem  of  interior  ballistics  in  1793.  He  introduced 
the  "Lagrange  approximation"  which  assumes  that  the  gas  velocity  at  any 
instant  increases  linearly  with  distance  along  the  tube,  from  zero  at 
the  tube  head  end  to  the  full  projectile  velocity  at  the  back  of  the 
piston.  It  is  further  assumed  that  all  the  propellant  charge  is  in 
gaseous  form  from  the  start  and  at  any  time  the  gas  density  is  the  same 
at  all  points.  It  can  be  shown  from  the  equation  of  continuity  that  if 
gas  density  is  independent  of  position,  the  velocity  distribution  is 
linear;  but  the  converse  is  not  necessarily  true. 

In  other  work,  Hugoniot  in  1889  used  the  theory  of  waves  of 
finite  amplitude  developed  by  Riemann  in  1858,  with  the  assumption  that 
all  the  propellant  was  completely  burnt  when  the  piston  began  to  move. 
He  followed  the  resulting  wave  of  rarefaction  on  its  journey  to  the 
tube  head  end.  The  method  was  extended  by  Gossot  and  Liouville  to 
follow  the  wave  as  it  travels  back  to  the  piston  after  being  reflected 
from  the  tube  head  end.  Finally,  Love  [4]  carried  the  analysis  as  far 
as  the  third  wave  traveling  toward  the  breech  and  Pidduck  [4]  applied 
Love's  solution  in  the  special  case  of  internal  ballistics.  But  all 
these  solutions,  though  completely  analytical,  hold  good  under  two 
important  assumptions: 
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(a)  Instantaneous  combustion. 

(b)  Adiabatic  expansion  of  each  element  of  gas. 

The  assumptions  may  be  applicable  for  the  devices  which  use  gaseous  fuel 
as  propellant,  say  automobile  engines,  but  for  the  devices  using  solid 
propellant  the  assumptions  are  far  from  the  real  situation.  In  this 
case,  gradual  burning  of  the  propellant  must  be  considered. 

Analytical  work  based  on  most  realistic  assumptions  has  been  done 
by  Carriere  [5].  For  simplicity  he  assumed  the  propellant  to  be  station¬ 
ary  in  the  combustion  chamber  at  the  time  of  burning  which  is  a  good 
assumption  for  cast  propellant  in  a  rocket-motor.  From  the  basic  con¬ 
cept  of  conservation  of  mass,  momentum  and  energy,  he  derived  three 
partial  differential  equations  expressing  gas  density,  gas  velocity  and 
entropy  as  a  function  of  time  and  distance.  He  transformed  those 
equations  into  three  ordinary  differential  equations  along  three  char¬ 
acteristic  directions  in  the  time-space  co-ordinate.  Then  with  proper 
choice  of  the  equation  of  state  for  the  combustion  gas,  he  followed  what 
is  commonly  known  as  the  "method  of  characteristics"  to  determine  the 
gas  properties  at  any  time  and  position.  The  effect  of  frictional  losses 
and  heat  loss  to  the  tube  wall  were  disregarded  in  the  analysis. 

The  problem  of  heat  loss  to  the  tube  wall  has  been  studied  by 
Hicks  and  Thornhill  in  England.  A  fairly  elaborate  description  of  their 
method  has  been  given  in  both  references  [l]  and  [2].  This  work  is  also 
based  on  the  Lagrange  approximation  of  linear  velocity  distribution  and 
uniform  gas  density  in  between  the  breech  face  and  the  piston. 

It  can  be  shown  that  at  high  velocity,  heat  is  mainly  transferred 
to  the  tube  wall  by  convection.  It  is  also  evident  that  a  boundary  layer 
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is  formed  at  the  inner  surface  of  the  tube.  The  heat  transfer  rate  per 

unit  area  through  the  boundary  layer  ;an  be  given  as  h(T^-Ts),  where  h, 

T  and  T  are  the  film  heat  transfer  coefficient,  temperature  of  the  gas, 
g  s 

and  temperature  of  the  inner  surface  of  the  tube  respectively.  All 

three  quantities  depend  on  time  as  well  as  position  along  the  tube. 

Hicks  and  Thornhill  considered  the  flow  in  the  boundary  layer 

to  be  the  same  as  the  flow  over  a  flat  plate.  In  internal  ballistic 

applications  the  flow  is  in  the  turbulent  region  most  of  the  time. 

Therefore,  they  used  the  analogy  solution,  as  extended  by  Von  Karman  to 

cover  Prandtl  number  other  than  unity,  to  obtain  a  relation  between  the 

heat  transfer  coefficient  h  and  wall  shear  stress  t  .  To  get  the  wall 

w 

shear  stress  they  first  found  a  "best"  power  law  for  the  velocity  profile 

(non-dimensionalized  with  respect  to  the  shear  velocity  /t  /p)  inside 

the  boundary  layer  which  was  capable  of  giving  the  local  wall  shear  stress 

r  within  three  per  cent  of  the  value  that  could  be  obtained  by  using 
w 

more  rigorous  logarithmic  form  of  the  velocity  p*  file  when  applied  to 
steady  and  uniform  flow  situations.  Then  they  used  the  boundary  layer 
momentum  integral,  including  the  terms  due  to  non-steady  and  non-uni  form 
nature  of  the  flow,  and  used  the  "best"  power  law  found  earlier  to  obtain 
the  local  wall  shear  stress  at  all  points.  The  heat  transfer  coefficient 
h  is  then  easily  calculated  from  the  analogy  solution.  They,  however, 
omitted  one  boundary  condition  that  the  boundary  layer  thickness  at  the 
base  of  the  piston  be  zero  at  all  times. 

The  heat  transfer  in  the  tube  ’/all  has  been  calculated  by  using 


the  differential  equation  for  unsteady  heat  conduction  with  proper 
boundary  conditions.  For  the  case  studied  by  Hicks  and  Thornhill,  i.e. 
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the  first  round  of  firing  from  a  cold  gun,  the  curvature  effect  of  the 
wall  was  neglected  as  the  temperature  rise  was  onfined  within  one 
millimeter  of  the  inside  surface.  Consequently,  there  was  no  heat  loss 
from  the  outer  surface  of  the  tube  which  remained  at  ambient  temperature. 
The  heat  conduction  along  the  length  of  the  barrel  was  also  neglected. 
Knowing  the  tube  material  properties,  namely  thermal  conductivity  and 
diffusivity,  it  was  possible  to  obtain  the  temperature  distribution  at 
the  inner  surface  of  the  tube  along  the  length  at  all  times.  The  free 
stream  values  of  the  gas  velocity,  density  and  temperature  were  taken 
from  the  onc-dimensional  ballistic,  solution. 

It  has  been  indicated  in  reference  [2]  that  frictional  pressure 
drop  is  small  compared  to  the  inertia  pressure  drop  needed  to  accelerate 
the  gas.  But  no  analysis  until  now  indicates  quantitatively  the  effect 
of  skin  friction  on  the  ballistic  properties.  Even  the  heat  transfer 
solution  has  not  been  fed  back  to  study  its  effect  on  the  one-dimensional 
solution. 


Present  Investigation 

In  the  light  of  available  theories,  it  is  clear  that  a  good  one¬ 
dimensional  solution  is  first  required  to  replace  the  Lagrange  approxi¬ 
mation,  or  at  least  check  its  validity  for  the  particular  problem.  The 
first  and  most  formidable  difficulty  in  writing  down  the  one-dimensional 
continuity,  momentum  and  energy  equations  during  the  burning  of  the  solid 
propellant  is  due  to  the  uncertainty  of  the  relative  velocity  between 
the  gas  pha^e  and  the  solid  phase.  It  extremely  difficult  to  estimate 
the  drag  exerted  on  the  burning  solid  particles  by  the  accelerating 
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combustion  gases.  Therefore,  two  limiting  cases  of  the  solids  velocity 
have  been  considered  in  the  present  work: 

Case  I.  The  solid  particles  move  at  the  same  velocity  as  the 
gas  phase. 

Case  II.  The  solid  particles  remain  at  their  initial  positions 
throughout  the  period  of  burning. 

For  both  cases  the  conservation  of  mass,  momentum  and  energy 

results  in  four  coupled  partial  differential  equations  expressing  volume 

fraction  of  solid  v  ,  gas  density  p  ,  gas  velocity  U,  and  pressure  P  as 
s  g 

a  function  of  axial  distance  x  and  time  t.  The  heat  release  due  to 
gradual  burning  of  the  propellant  is  ta.^en  into  account.  Atypical 
propellant  geometry,  namely  a  hollow  cylinder,  is  considered  whereby 
the  total  burning  surface  remains  constant,  although  this  assumption  is 
not  essential. 

The  ballistic  properties  at  the  internal  points  ate  calculated 
from  these  equations  after  writing  the  same  in  finite  difference  form. 

But  to  obtain  the  properties  at  the  two  ends,  namely  the  tube  head  end 
and  the  piston  base,  the  equations  are  transformed  into  ordinary  differ¬ 
ential  equations  along  the  characteristic  directions.  The  covolume  of 
the  gas  is  assumed  to  be  constant,  and  experimental  data  for  burning 
rate  is  used.  As  one  of  the  initial  conditions,  it  is  assumed  that  the 
piston  does  not  start  until  a  certain  specified  pressure  is  reached 
inside  the  chamber  and  thereafter  the  piston  does  not  experience  any 
resistance  to  motion. 

The  boundary  layer  momentum  integral  for  a  non-steady,  non-uniform, 
developing  flow  inside  a  tube  is  derived.  The  profile  shape  factor  H 
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(ratio  between  the  displacement  thickness  6  and  momentum  thickness  9) 
is  introduced  and  the  Ludwieg-Tillmann  [6]  friction  factor  is  used.  As 
a  first  approximation,  the  shape  factor  is  assumed  to  be  constant  in 
the  present  work.  The  flow  is  in  the  high  Reynolds  number  region  for 
which  the  boundary  layer  thickness  is  small  compared  to  the  tube  radius. 
It  is  therefore  legitimate  to  replace  the  free  stream  values  of  gas 
density  and  velocity  by  the  values  obtained  from  the  one-dimensional 
solution  neglecting  the  boundary  layer  thickness. 

The  local  heat  transfer  coefficient  h  is  calculated  by  using 
Colburn’s  analogy  [7]  between  heat  and  momentum  transfer.  It  covers 
Prandtl  numbers  other  than  unity  and  is  simple  to  use.  The  values  of 
viscosity  and  gas  density  at  the  film  temperature  are  used.  The  heat 
transfer  in  the  tube  wall  is  computed  from  the  unsteady  one-dimensional 
(radial)  heat  conduction  equation  with  appropriate  boundary  conditions. 
The  wall  temperature  is  also  found  as  a  function  of  axial  distance  and 
time. 

The  heat  loss  term  is  entered  into  the  one-dimensional  energy 
equation  and  a  comparison  of  ballistic  properties  is  made  with  the  solu¬ 
tion  without  heat  loss.  Effect  of  wall  shear  stress  is  also  included. 
The  ballistic  efficiency  of  the  piston-cylinder  arrangement  is  compared 
by  varying  different  design  parameters. 
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SECTION  III 

MATHEMATICAL  ANALYSIS 

The  mathematical  analysis  consists  of  two  major  parts: 

0-)  One-dimensional  analysis  with  gradual  burning  of  the  solid 
propellant,  including  the  effect  of  heat  transfer  and  skin 
friction. 

(2)  Formulation  of  the  boundary  layer  problem  and  determination 
of  heat  transfer  to  the  tube  wall. 

As  outlined  in  the  previous  chapter,  the  present  analysis  is 
carried  out  for  two  extreme  cases  of  solid  velocity.  In  the  first  case, 
it  is  assumed  that  a  burning  solid  particle  moves  with  the  same  velocity 
as  the  combustion  gases.  In  the  second  case,  however,  the  solid  particles 
are  assumed  to  be  stationary  at  their  initial  positions  throughout  the 
period  of  burning.  Henceforth  these  two  cases  are  referred  as  Case  I 
and  Case  II,  respectively. 

One-Dimensional  Analysis  Including 
Heat  Transfer  and  Skin  Friction 

Case  I 

The  assumptions,  other  than  that  regarding  the  solids  velocity, 
which  are  made  to  simplify  the  model  are  as  follows: 

(1)  At  any  instance  of  time,  the  linear  speed  of  burning  r  £s 

b 

same  for  all  the  solid  particles  and  it  is  a  function  of  the 
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average  pressure  in  the  chamber  (space  in  between  the  tube 
head  end  and  the  piston  base). 

(2)  The  solid  propellants  are  single  perforated  circular  cylin¬ 
ders  in  shape  whereby  the  total  burning  surface  remains 
constant  during  the  whole  period  of  burning. 

(3)  The  burning  rate  is  fast  enough  to  consider  that  the  temper¬ 
ature  of  the  remaining  solids  at  any  instance  of  time  remains 
constant  at  the  initial  temperature. 

(4)  The  propellant  material  is  incompressible  and  its  coefficient 
of  thermal  expansion  is  negligible. 

(5)  The  piston  starts  to  move  only  when  the  chamber  pressure 

reaches  a  certain  value  P  ,  and  thereafter  the  resistance  to 

o 

its  motion  is  negligible  compared  to  the  pressure  force 
exerted  on  it  by  the  combustion  gas~s  in  the  chamber. 

The  conservation  equations  are  as  follows  (for  derivation  see 
Appendix  I ) : 

Solid  continuity: 


0 


(1) 


Gas  continuity: 


3t 


A 

+  U3^  + 


(l-,8) 


3U  =  (ps~pg) 

9x  (1— v  ) 

s 


(2) 


Momentum: 


(3) 


l£  +  u~-_  —  w 

3t  3x  p  3x  p  R 
m  m 


Energy: 


Dh  Dh 

s  .  / i  s  g  DP 

v  o  -  4*  (1— v  )  d  — “•  —  — 

s  s  Dt  v  V  pg  Dt  Dt 


p  2h  2t  U 

=  p  (W  +  —  -  h  )v  , - ~(T-T  . .  +  -~ 

3  PD  g  d  R  w,i)  R 
s  s 


where  is  the  volume  rate  of  decrease  of  solids  per  unit  cylinder 
s 

volume  and  is  given  by; 


. bt,  ',s(x't)rb 

vd  <*,t>  ■  (— yjrr - 

s  P  /  p 

!  vg(x,t)dx 
Jo 


The  equation  of  state  of  the  gas  is; 


P(v  -  n)  -  R  T 
g  g 


P(—  -  n)  -  R  T 
P  g 

g  8 


where  the  gas  constant  is  obtained  from  the  ratio  of  the  universal 
gas  constant  R^  and  the  molecular  weight  of  the  gas  M. 

It  has  been  shown  in  Appendix  I  that  under  assumptions  three  and 


four  as  stated  earlier,  the  differential  of  enthalpy  of  solids  per  un;Lt 


20 


mass  hf  and  the  differential  of  enthalpy  of  gases  per  unit  mass  h  can 
J  8 

be  given  by: 


dh  -  —  dP 
3  ps 


(y-np  ) 

dh  =  - - fr  dP  -  - rr - 7  dp 

g  p  (y-1)  (y-1)p  /  g 

S  8 


Substituting  equations  (7)  and  (8)  into  the  energy  equation  (A); 


(l-vs)(l-nPg)  op  _  <1-^s)yP  ^ 

(y-1)  Dt  (y-1)p  Dt 


-  p  (W  + - h  )v  , 

s'  p  g  d 
s  6  s 


2h.  2t  U 

(T-T  , )  +  -~— 
R  w,i/  R 


Using  gas  continuity,  i.e.  equation  (2)  to  replace  in  equation 
(9)  the  final  form  cf  the  energy  equation  becomes: 


q-v^a-nPg)  DP  yF  an  vffe.-o,)  , 

(Y-1)  Dt  (v-1)  3x  (Y-i)p  Vd 

g  s 


p  (W  +  —  -  h  )v.  -  — ”(T-T  ,) 

s  p  g  d  R  w,i 
s  s 
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or, 


3P 

3t 


+ 


"r  + 


3JJ 
I  3x 


2h. 

C  v  -  E  — -  (T-T  .) 
Id  I  R  w,i 

s 


+  E 


I 


2t  U 
w 

'  R 


(10) 


where 


BI  (1-v  ) (1-np  ) 
s  g 

YP(p  -p  )  +  (Y-I)p  P  (W  +  —  -  h  ) 

s  g  S  g  Ps  g 

CX  =  Pg{(l-vs)(l-nPg)} 


(11) 


(12) 


=  (Y-l) 

I  (1-v  ) (1-np  ) 
s  g 


(13) 


The  initial  conditions  of  the  conservation  equations  are: 
Position  of  the  piston,  1.^(0)  “  Lq 


U (x,0)  =  0  ;  P(x,0)  =  P  ;  T(x,0)  =  T  ;  p  (x,0)  =  P 

w  O  5- 


(14) 


and  v  (x,0)  **  v 
s  ’  s 


at  0  <  x  <  L 

—  —  o 


where  P  is  the  pressure  at  which  the  piston  starts  to  move,  and  T  is  the 
o  r  o 

adiabatic,  stagnation  flame  temperature  of  the  propellant.  By  knowing  PQ 
and  T0  it  is  possible  to  determine  p^from  the  equation  of  state  (2.6): 
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1 


(15) 
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(19) 


The  unknowns,  and  the  corresponding  equations  from  which  they 
can  be  calculated  are  listed  below: 

Unknown  Equation 

Volume  fraction  of  solids,  Solid  continuity,  (1) 

Gas  density,  p  Gas  continuity,  (2) 

S 

Velocity,  U  Momentum  equation,  (3) 

Pressure,  P  Energy  equation,  (10) 

Gas  temperature,  T  Equation  of  state,  (6) 

The  conservation  equations,  i.e.  (1),  (2),  (3)  and  (10)  are 

written  in  finite  difference  form,  and  a  numerical  scheme  which  takes 

into  account  both  forward  and  backward  space  derivatives  are  used  to 

calculate  the  corresponding  unknowns,  i.e.  v  ,  p  ,  U  and  P,  at  all  the 

s  g 

interior  points  at  an  advanced  time  by  knowing  the  present  values  at  and 
around  those  points.  The  gas  temperature,  T,  is  then  calculated  from 
the  equation  of  state  (6),  The  details  of  the  solution  technique 
shall  be  discussed  in  Chapter  IV. 

The  above  solution  technique,  however,  is  not  applicable  to  the 
boundary  points,  i.e.  the  piston  base  end  and  the  tube  head  end,  as  space 
derivatives  on  both  sides  of  these  two  points  are  not  available.  This 
necessitates  the  transformation  of  the  conservation  equations  tr,  ordinary 
differential  equations  along  characteristic  directions,  i.e.  to  follow 
the  "method  of  characteristics"  [8]. 
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P-U  Charac  .ristic.  The  energy  equation  (10)  is: 


ii  +  ul-  +  B  •  =  C  v  -  ET  (T-T  .) 
3 1  ax  I  3x  Id  I  R  w,i 

s 


Multiplying  the  raon'entum  equation  (3)  by  an  arbitrary  constant  X: 


,  3P  ,  ,  3U  .  ,  lt3U  t  vj 

X— —  +  Xp  - —  +  Xp  lb? —  =  —  X  — — 
3x  m  3t  m  3x  R 


Adding  equation  (20)  to  equation  (10): 


[l  +  <U+X>  f]  +  [Xp.  17  +  <BI  +  XpmU>77_ 


CI  -  EI  T  ^w.i) 
s 


r  [V-X] 


To  obtain  the  characteristic  directions,  the  value  of  X  shall  be  such 
that: 


dt  ,  _1_ 

dx  "  U+X  “  3  +Xp  U 
I  m 


Ai,2  =  1 
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Dividing  equation  (2.21)  by 


and  using 


1  3_  (J+A)  _  d 

4+(u+a)z  3t  /i+(u+a) 2  3x  dn^ 


where  n  corresponds  to  \\  ,  i.e.  positive  sign  of  x>  and  £  corresponds 
to  \2  ,  i.e.  negative  sign  of  X,  the  equation  (21)  becomes. 


dP  dU  _  1  r  . 

dn,f.  l,2Pm  dn,5  '  /JVYu+j^  L  Ivds 


Zll.  ZT 

-  h  -r<T-I„(1> +  ir(Eiu-x)J  <23> 


An,S  “  /(Ax)z  +  (At)^  3  At  /ET(U+X)‘ 


Therefore,  along  n-characteristic,  i.e.  3  - 

dX  b  +  /B  /p 
I  m 


_ _ _  r  2h 

AP  +  p  /OfT  AU  =  C  v  -  Et  — -  (T-T  ,) 
mlra  Lid  I  R  w,i 


+  T  (Ei"  '  'TO]  at 


and  along  ^-characteristic,  i.e.  —  -  - — 


U  -  /B  /p 
1  m 
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r  2h , 

-  p  v^'r/p  AU  =  C  v  -  Et  — -  (T-T  .) 
m  I  m  Lid  IR  w,i 


+  T  <ex"  +  "W] 


P-p  Characteristic.  By  rearranging  equation  (9) f 


3p  3p 

- S.  4.  u — & 

at  3x 


r»tI)»i  +  (T-D»sp.(hg  -  ps  - u) 


(l-vg)yP 


(y-1)po  f2h  2t  U 

q. - . — £L  — i(T_T  ) - — 

(1-vs)yP  L  r  w,i'  r 


Dividing  this  equation  by  /l+Uz  and  using, 


1  3_  U  3_  =  d_ 

/I«F  dt  /IW  9x  =  d; 


As  ■*  / (Ax) z+(At)^  =  At  /I+U7 


,  ,  ,  ,  .  dt  1 

along  a  particle  path,  i.e.  ^  ^ 


„  n  ,  n  (h - W) 

P,.(l-hP  )  (y-1)p  p  g  p 

AP  +  - 5 -  v  At 

YP  (l-vs)yP  dg 


(y-l)p  r 2h  2x  D 

+  TPofF  L-r  "-Li’-f”  At 

s 
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vs  pg  Characteristic.  From  equation  (2) 


3U 

3x 


(1-vJ  f3p  3p  1 

-  £  +  Ur— “ 
3t  3x 


(P  ~pJ 

+  - 


P  Vd 
g  s 


(29) 


9  U 

Substituting  this  expression  for  —  in  solid  continuity  (1)  : 


3v 

3v 

v  (1-v  ) 

'£p  3p  i 

r  v  p  +(i-v  )p  -i 

s  . 

_ s 

s  s 

— &  +  u~~& 

L  at  3x 

i 

S  S  S  B 

at 

3x 

p 

T 

p 

g 

J 

g  J 

•=  0 


(30) 


Proceeding  in  the  same  fashion  as  for  the  P-p 
obtains: 

Along  a  particle  path,  i.e.  —■  »  -jj  : 


g 


characteristic, 


one 


v  (1-v  )  p 

Av  »  ~  Ap  -  —  v.  At  (31) 

s  P„  g  p  d 

g  g  s 

The  procedure  of  solving  the  above  characteristic  equations  are  discussed 
in  Chapter  IV. 

Case  II 

Ir.  this  case  the  solid  propellant  particles  are  assumed  to  be 
stationary  at  their  initial  positions  throughout  the  period  of  burning. 

The  linear  speed  of  burning  r^,  is  same  for  all  the  solid  particles  and 
is  a  function  of  the  average  pressure  In  the  space  between  the  tube  heal 
end  and  the  initial  position  of  the  piston  Lq.  The  rest  of  the  assumptions 
are  the  same  as  those  for  Case  1. 

The  conservation  equations  in  this  case  are  (See  Appendix  I  for 
derivation) : 
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Solid  continuity: 


3  v 

Jc~  +  vd  '  0 
s 


(32) 


Gas  continuity: 


3p  SII 

_E  +  u— S.  +  p  — 
dt  dX  g  dx 


+  ^-U 


3v 


(p  “P  ) 

s  e  _ _ 

(l-v  )  d  '  (1-v  )  dx 
s  s  s 


(33) 


Momentum: 


$  +  1^ 

<3t  9x 


PSU 


2x 


1_3P _ _ 

p  3x  (l-v  )p  Vd  (l-v  )p  R 
8  s/Hg  s  s^g 


(34) 


Energy: 


DhK  DP  p 

(l-v  )p  — -  (l-v  )—  -  p  (W  +  —  +  -r - h  )v 

ft  f?  Dt  s  IJt  s  n  Z  Q 


s  Dt  s  p  l  go 
s  s 


2h  2x  U 

+  -T- 


(35) 


As  none  of  the  solid  parti'.les  moves  beyond  Lq,  can  be  expressed  as 

s 


(36) 
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The  same  equation  of  state,  i.e.  equation  (6)  is  used  and  by  succes¬ 
sive  use  of  equation  (8)  and  (33)  the  final  form  of  the  energy 
equation  (35)  becomes: 


!U.C  * 

St  Sx  II  Sx  II  d 


BIIU 


Sv 


(1— v  )  Sx 
s 


-  E 


2h 

II  ~R 


2t  U 

1  +  hi  -f- 


(37) 


where 


II  (l-np  ) 

6 

p  u2 

yP(p  -p  )  +  (y-l)p  p  +  p  +  2 
r  S  g _ s_g _ s _ 

II  p  {(1-v  )(l-np  )} 

g  s  g 


(38) 


(39) 


E 


II 


E  =  (Tr1).  .. 

I  (l-v  )(l-np  ) 

s  g 


(40) 


The  initial  and  boundary  conditions  are  the  same  as  those  for  Case  I. 

The  characteristic  equations  are  also  required  to  calculate  the  ballistic 
properties  at  the  two  ends. 

P-U  Characteristic.  The  procedure  is  exactly  same  as  Case  I.  Multiplying 
the  momentum  equation  (34)  by  an  arbitrary  con^ant  A,  and  adding  to 
the  energy  equation  (37)  one  obtains: 
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[If  +  <U+X)f  ]  +  K  If  +  <BII+X‘>8®lf] 


BttU  8v  p  U 

„  .  ,  II  s  ,  s 

=  C„v,  +  -7- - - - A-; 


II  d  (1-v  ,  3x  (1-v  )  d 
s  s)  S  s 


2h , 


2T  r 


-  E  -^(T-T  ,)  +  -=* 


'll  R  w,i/  '  R  L  II  (1-v  ) 

s 


E.JJ  - 


■]  (41) 


The  characteristic  directions  are  such  that: 


Ap, 


dt  1  _____ 

dx  U+A  BTT+Ap  U 

II  8 


X  =  +  /!ii  =  +  / _ IL 

1.2  -Jp„  -/  pgd-npg) 


-  ±a 


(42) 


From  solid  continuity,  i.e.  equation  (32) 


3v 

s 

•  -  -  S»  — V 

at  d 


Bnu 


s 

3v 


(43) 


By  adding  and  substracting  ^  ^  ^  (U+A)  3"t~"  on  rl^lt  hand  side  of 

s 

equation  (41) 


[f +  <u+x>!f]  +  \  [If +  <u+x>f. 


Bn“ 


r  3v  3  ;  -t  B__U  3v 

[  jt  +  -  11  8 


(1-Vg)(u+A)  L  3t  ■  '“"''ax  J  (l-Vg)(U+A)  at 


P  U  2h  2t  r 

+  C  v  -  A— - .i  -  E_ ~(T-T  .)  +  EttU-t~ -J 

11  ds  (1"Vs)  ds  11  K  W»1  R  11  d-'-aO 


s 

(44) 
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Dividing  equation  (2.44)  by  /l+(U+X) ^  and  using 


1 _  3_  (U+X)  3_  _  d_ 


A+(u+x)2  St  /i+(u+\) 2  3x  dn,? 


yields  the  relation: 


BttU  dv 

UP _  .  dU _  _  II _ s_ 

dn,€  l,2pg  dn.C  (l-vs)(U+X)  dn,C 


p  ux 


B„U 


l  fr „  Ks“"  ,  “ir  i  .. 

/T  V+  /  1  \  /lT  .  \  \  •»  V, 

mTvTZ  - 


/iW^11  (1'V  U-vs)(u«)'  ds 


2h  2t  ,  -i 

-  hiT^w.P  +  T  lEnu  *  o^TT1]  W5) 


d  t  X 

Therefore,  along  q-characteristic,  i.e.  —  ■  : 


bii“ 

4P  +  °S  a  4U  ‘  (l-v8) (u+a)  4vs 


r  psu  Bnu  1  . 

+  LCII  "  3  (l-v  )  +  (l-v  ) (U+a) J  Vdo  At 


-  E, 


2h 

II  ~R 


1  <T-Tw,i)st  +  qr  [EIIU  -  O^r]  4t 


(46) 


dt  1 

and  along  ^-characteristic,  i.e.  ^  : 
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BIIU 

4P  -  pg  a  AU  ■  (i-"')(u-7)  Avs 

s 


+  r  +  PSU  +  BIIU  1  .  , 

LCII  a(l-0'  (l-v  )(U-a)J  Vd„ 


2h .  2t  r 

-  EII  R  (T'Tw,i)At  +  IT  LEIIU  +  o^r] At 


P-Pg  Characteristic.  Using  equation  (8)  in  equation  (35)  an  alter¬ 
native  form  of  energy  equation  is: 


3 1  3x  YP  L  3 1  3s  j 


(y-Dp  p  (h  -  w  - 

s  g  g  Pg  2 

+  (l-v  )YP  Vds 


(Y-DP  f"  2h  2x  U I 

- 6_  — i/T-T  ) _ — 

(1-J  )YP  L  R  w,i'  R  J 


Proceeding  in  exactly  the  same  manner  as  for  Case  I,  along  a  particle 

,  dt  1 
path,  i.e.  ^  -  tj  : 


.  n  4  (y-1)p  p  (h  -  w  -  ~ —  y-) 

c>  U-np  )  s  g  g  p  2 

Ap  .  fi  -  &  AP  +  - TT - r— - 2 -  V.  At 

g  YP  (1-v  )yP  d 

S  9 


(Y-l)p_  r2h 


(l-v  )yP 
s 


r2h  2x  U-j 

— i  (X-T  ) - —  At 

L  R  '  v,r  R  J 
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vg  at  the  piston  base  is  always  zero  after  the  piston  starts  moving  and 
vg  at  the  tube  head  end  can  be  obtained  from  equation  (32)  alone. 

Boundary  Layer  Analysis 

The  boundary  layer  part  of  tue  entire  analysis  did  not  receive 

much  attention  in  the  past  because  of  the  nonsteady  and  nonuniform 

nature  of  the  free  stream  flow.  The  flow  is  generally  in  the  turbulent 

region  with  pressure  gradient  in  the  direction  of  the  flow  and  a  large 

temperature  difference  across  the  boundary  layer.  Also,  in  Case  I  a 

gas-solid  mixture  flows  down  the  tube;  hence  the  analysis  is  more 

complicated.  A  number  of  attempts  [8,  10,  11,  12]  have  been  made  in 

the  past  to  model  the  mechanism  of  heat  transfer  in  a  gas-solid  mixture 

with  various  solid  particle  sizes  and  loading  ratios  (w  /w  ) .  It  has 

s  g 

been  found  that  the  effect  of  the  solids  on  convective  heat  transfer  is 
prominent  for  micron-size  particles  whereas  for  millimeter  size  the 
effect  is  not  appreciable.  The  present  problem  deals  with  the  solid 
prope’lant  of  millimeter  size  and  most  of  the  time  it  burns  out  completely 
long  before  the  piston  reaches  the  end  of  the  tube.  It  has  also  been 
found  from  the  study  of  Hicks  and  Thornhill  [2]  that  the  boundary  layer 
'■h/.ckness  is  small  compared  to  the  tube  radius.  Therefore,  to  simplify 
the  model,  it  is  assumed  that  the  solids  always  stay  in  the  core  of  the 
flow  and  never  enter  into  the  thin  boundary  layer  at  the  wall. 

In  the  present  study,  an  integral  approach  is  preferred  to  a 
differential  approach  to  keep  the  model  relatively  simple  and  tractable. 
The  boundary  layer  momentum  integral  for  the  nonsteady  and  nonuniform 
compressible  flow  inside  a  tube  as  derived  in  Appendix  II  is: 


p  (U  ~u) r  dr  + 


r  R  1  rR 

/  pu(Um-u)r  dr  +  I  p(UM-u)r  dr 

R-0  -*  II- G 


/  r  2  r\n  3U  3U 

,.0  3?  00  a 

=  R6-  —  —  +  p  ~  +  p  U  -r— 

2  3x  f  3t  £  00  3x 


+  i  R 
w 


Defining  5  =  Displacement  thickness 


and  ')  =  Momentum  thickness 


such  that, 


rR 

/  J  - 


)  -U  I  ^Ti  I 

£  *  JL 


i2nr(U  -u)dr 


-  1r) 


Pr (U  -u)dr 


pfu!  f  2n 


2nr  dr  ■  f  p  2nr  u(U  -u)dr 

I  uo 

R-0  R-5 


_ 

PfU>(l  "  |r)  =  J  P  u  *(Uw-u)dr 


and  using  the  definition  of  the  profile  shape  factor  H  =  —  ,  and  for 

*  0 

a  thin  boundary  layer  <<  1,  <<  1,  <<  1,  the  momentum  integral 

equation  (50)  becomes. 
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+  |j[c£u^r6]  +  (Pf'J.RHeh 


=  R6  —  +  p  -T—  +  p  U  - -  +  T  R 

3x  f  3t  f  30  3x  w 


AA  ATI  W.  a  O  A  A 

p-U  RH~  +  p.U  RF7  +  p,RH(k~  +  U  RH&~-  +  pcU  IF 
£  ®  3t  f  °°  3t  f  3 1  00  3t  f  00  3x 


8Uoo  2  3pf  3U« 

+  2p,U  RSr--^  +  U  R(h— -  +  pcU  RH0t— 

f  »  3::  «  3X  f  ®  ax 


Tad  SU  3U  -I 

R6  —  +  p  -r— -r  p  U  t— >  +  t  R 

3x  f  3t  f  “ox  w 


Dividing  equation  (53)  by  p  U  RH0  : 

t  CO 


a  3t  II  Dt  U  3t  pr  3t  HO  3x  H  3x  prH  3x  3x 
00  i  f 


,  r,_  3U  3U  1  i 

o _  ap  “*  ,  „  °°  .  w 

,,U  KO  Ux  +  Pf  at  +  PfU“  3x  +  p.U  P 

f  CO  L  J  r  f  nn 


From  the  study  of  steady  compressible  turbulent  boundary  layers  by 

Reshotko  and  Tucker  [l3],  it  is  likely  that  for  moderate  Mach  number 

flow  encountered  in  this  problem  (M  <  1.5),  the  percentage  change  in 

1  3H 

the  shape  factor,  i.e.  —  —  ,  is  small  compared  to  the  percentage  change 

H  at 

1  30 

in  momentum  thickness,  —  —  .  As  a  first  appro Lmation,  therefore,  the 

6  d  t 

shape  factor,  H,  is  assumed  to  be  a  constant.  A  more  rigorous  approach 
would  be  to  derive  another  auxiliary  equation,  say  moment  of  momentum 
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[_  3H 

14 J  to  obtain  an  expression  for  —  .  However,  derivation  of 

such  an  equation  for  the  nonsteady  case  is  extremely  complicated  and 

therefore  neglected  in  the  present  work. 

For  thin  boundary  layers,  =  U  ;  1^  -  T  and  the  film 

temperature, 


(55) 


The  gas  density  at  the  film  temperature,  p  ,  can  be  evaluated  from  the 
equation  of  state  (6),  and  the  final  form  is: 


Equation  (54)  finally  becomes: 


The  initial  condition  is:  6(x,0)  =  0. 

The  boundary  condition  at  the.  piston  end  is  0(L^,t)  a  0,  which  is 
obvious  from  the  fact  that  all  the  particles  at  the  piston  base  are  at 
the  full  piston  velocity  all  the  time.  The  condition  at  the  tube  head 
end  shall  be  established  later. 

It  is  assumed  t  iat  the  entire  flow  is  in  the  turbulent  region 
and  the  wall  shear  s  cess  can  be  obtained  from  the  Ludwieg-Tillmann 
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friction  factor  [6],  which  was  developed  fun  a  series  of  experiments 
with  ail  types  of  pressure  gradients.  The  original  expression  which 
holds  good  for  incompressible  flow  with  small  tetper.iture  differences 
across  the  boundary  layer  is : 


cf  hp  v 


^  .  0.246  x  IQ-0'678"  (!**!)-°-268  («) 


In  the  present  work,  the  expression  is  slightly  modified  by  using 

the  fluid  properties  (p,  p)  at  the  film  temperature,  T.,  instead  of  the 

free  stream  temperature,  T^,  to  take  into  account  the  effect  of  property 

variation  across  the  boundary  layer.  The  expressions  ror  local  friction 

coefficient,  C-,  and  local  shear  stress  at  the  tube  wa] 3 ,  x  ,  used  in 
’  f  w 

the  present  work  are: 


C 


f 


2A 


1 

(Re/ 


T 

W 


a  pfir 


i 

(Re/ 


(39) 


where 


0.123 


10 


0.678H 


Pfue 

B  =  0.268;  Re„  =  — — 
0  uf 


Using  the  above  expression  for  wall  shear  stress  in  equation  (57)  ,  and 

g 

multiplying  equation  (57)  by  (1+B)9  : 
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u  30  .  n+R)A^f_  „(1-B) 
H3x  +  (1+B)H  ~  U 
Pf 


-  (1+8)0 


1  3pf  U  3pf  1  au  (H+2)  3U 
pf  3t  pfH  3x  U  3t  H  3x 


+  (1+B)  C,  +  ^  +  f 

1  prU  dx  U  3t  3x 

L  f 


ahere 


0=0 


(1+B) 


or,  0=0 


1 

(1+B) 


Cl  "7 


At  the  tube  head  end,  U  ■  0  and  the  equation  (60)  becomes: 


30  /,.„x  rTl  8pf  U  9pf  ,  1  3U  ,  (H+2)  3U 

at  "  "  (J  *  0  pf  at  +  pfH  3x  +  U  3t  +  H  3x 


,  .  1  3P  1  3U  3U 

+  (1+B)  G  C  — r  —  +  —  ~  +  r— 

1  (_PfU  3X  U  dt  dX 


and  at  t  ■  0,  0  •  0,  This  implies  that  at  the  tube  head  end,  0,  i.e. 
momentum  thickness  or  boundary  layer  thickness  is  zero  at  all  times. 

The  equation  (60)  is  applicable  to  both  Case  I  and  II  for 
computing  momentum  thickness  0,  and  thence  the  friction  coefficient, 
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C^,  at  each  station  in  the  axial  direction  at  each  time  step. 


Heat  Transfer  Analysis 

As  the  flow  is  in  the  turbulent  region,  the  analogy  between  the 
momentum  transfer  and  the  heat  transfer  provides  the  easiest  way  to 
determine  the  heat  transfer  coefficient,  h.,  at  the  tube  wall.  Because 
of  its  simplicity,  Colburn's  analogy  [7]  has  been  used  for  Prandtl 
numbers  other  than  unity  as  follows: 


St  Pr 


2/3 


C„ 

1 

2 


or. 


hi  "  pfUcp  <2 ^)/Pr2/3  <63> 

u  c 

where  Pr  =  (—%—£)  and  all  properties  are  evaluated  at  T^. 

Kg 

The  heat  transfer  in  the  tube  wall  is  considered  as  a  one¬ 
dimensional  (radial)  unsteady  heat  conduction  problem  in  a  hollow 
cylinder.  Longitudinal  heat  conduction  is  neglected  because  the  tem¬ 
perature  gradient  in  the  radial  direction  is  expected  to  be  steeper 
by  several  order  of  magnitude  than  in  the  axial  direction.  The  differ¬ 
ential  equation  can  be  written  as  [15]: 


5T 

w 

(1 

a2T  . 

w  1 

9T 

w 

3t 

3r 

— 1 

w 

^9r2  r 

(64) 


The  material  properties  will  be  assumed  constant  with  respect  to  tem¬ 
perature.  The  boundary  conditions  are: 
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at  r  =  R,  q .  =  h  (T  -  T  . ) 

’  1  i  *  w,i 


5T 


-  K 


W 


w  9r 


r=R 


3T 

at  r  =  R  ,  q  =  h(T  -T  ,  )  =  -  k  . 

o  o  o  w,o  amb  w  9r 


w 


(65) 


r=R 


at  t  =  0,  T(r)  =  T 


amb 


It  is  possible  to  solve  equation  (64)  numerically  and  obtain  the 

temperature  at  the  inner  surface  of  the  tube  T  .  at  each  station  along 

•  w ,  x 

tne  length  of  the  Lube  at  each  time  step.  The  iocal  heat  transfer  rata 

to  the  wall  per  unit  surface  area  is  given  by  h  (T-T  , ),  and  integrating 
*  i  w,  i 

over  the  entire  surface  and  the  time,  the  total  heat  transfer  to  the 


tube  wall  can  be  determined.  The  values  of  local  wall  shear  stress  x  , 

w 

heat  transfer  coefficient  h. ,  and  inner  surface  temperature  T  as 

x  W  ,  1 

calculated  from  (59) ,  (63)  and  (64)  are  used  in  the  one-dimensional 

analysis  for  the  subsequent  time  step. 


Non-dimensionalization 

Before  proceeding  to  the  solution  technique  that  can  be  applied 
to  solve  che  equations  derived  so  far,  it  is  advantageous  to  non- 
dimensionalize  the  equations  to  obtain  a  general  solution  for  the  geo¬ 
metrically  similar  devices  with  the  same  initial  conditions.  The  non- 
dimensionalized  parameters  are: 


Axial  distance, 

x' 

X 

ae  ■  -  - 

Lt 

Pressure, 

P' 

=  P/P 

♦ 

Temperature, 

T' 

H 

H 

n 
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Density , 
Velocity, 

Time , 

Linear  speed  of 
burning , 

It  can  be  noted  that: 


P*  =  P/Pg 


u*  =  u/u 

„  _o 


U 

o 


•RJi. 


Vu 


o 


(66) 


v  r/ 
s  b 


L1 

l  dx,J 
s 


U 


o 


h 

g 


P'(Y_nPg  Pg} 

p'(y-I) 

g 


(67) 


(68) 


W 


R  T 


-JLJL 

(Y-l) 


R  T 


_S^_ 

(Y-l) 


A  PfU  (Refl) 


1^'po 


'Q' 


(69) 


(70) 


and,  np  is  a  constant  non-dimensional  quantity. 

^o 

Finally,  non-dimensional  forms  of  the  conservation  equations  are: 
Case  I 

Continuity  of  solids: 


3v 

s 

3t* 


+  IT 


3v 

3x 


T  +  v 


fr  +  v  ’  -  0 

s  3x  d 

s 


(71) 


Continuity  of  gases: 
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3P  '  3pe‘  Pk'  9U 1  (ps'“pkT) 

at*  3x '  (l-v  )  3x '  (1-v  )  Vd 

s  s  s 


Momentum: 


3U*  ■  w.  3U’ 
at*  3x' 


1  IF1  ,  °£'U'  1 

■>„’ 3x' '  1  '  %'  ^7 


Energy : 


ip;+u.|pi  +  Bt,|ui  =  c  .  * 

3t  dx'  I  3x'  Id 


where, 


P'U'3 

'  D2,E,hi(T-Tw,i>  + 


&i  ~  "(l-v  )(l-np  p  ’) 

S  g0  g 


Yp'(p  '-P  ,)+(Y-1)p  'p  '(W  +  P  ' 

p  M(i-v  yor-vp  P ')> 


E 1  = 


<l-vs)(l-nPg  Pg  ) 


V  "  “<r> 


V  '  20  pTT 

o  o 


A3 


Equation  of  state: 


U-np  p  ') 

P'  _ 8 

p  '  <1-rlp  ) 

8  go 


(80) 


Characteristic  equations: 

Along  n',  £'  characteristic,  i.e.  “y  =  ^ 


1  m 


AP 


*  ±  p  • /bt * /p  '  AU'  «  CT'v.  '  At  -  D  'E’h, (T-T  .)At' 
m  1  m  Id  2  I  w,i 


+  V 


E*U'  +  /bt7p  ' 
I  m 


Pf'u'2 

(Re0) 


At' 


(81) 


dt'  1 


Along  a  particle  path,  i.e.  =  -jjr  , 


Ap  '  =  G'  AP'  +  HT'v,  'At'  +  D„ ' 
K  Id  2. 


(Y-1)Pk’  ■ 

(1-v  )yP' 


h . (T-T  ,)At' 

i  w,i 


Di’ 


(Y-l)p. 


»  1  «  »n«3 


(1“V  )yP' 


a'1*  J  (Re0) 


p  'U' 

-2 -  At' 


,  B 


(82) 


where, 


G' 


pg’<1-nV*,) 

yP' 


V 


/  n  ,  ,  (h  '  -  W'  -  — ~t) 

(y-DPs  pg  s  Ps 

(l-v^yP1 


(83) 


(84) 
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,  u  .  dt'  1_ 
Again  along  a  particle  path,  i.e.  , 


v  (1-v  ) 


s  s 

A  V  =  - “i — 

s  P„ 


Ap  *  -  ~7  V  'At' 

B  P  a 

b  g  S 


Case  II 


Continuity  of  solids: 


Continuity  of  gases: 


Dv 

S  ,  •  f 

3r  +  Vd 

s 


3P  '  3P  '  .Mi* 


p  'U'  9v 


a  l  ,3U'  >rs  '  ...  ,  .  18 _ s 

^  +  Uj^+P857°  (1-v)  "  Vdq  +  (1-v  )  3x' 


Momentum: 


3U1  +  U'i^ 

at'  u  ax' 


Pf‘u’2 


r  +  u  aT5"  “  "  77  ax' 


(1-v  )p 


Vd_  "  D1  a-vs)Pp,(Ho0;B' 


s'  g  s 


Energy : 


B  1 U 1  3v 

*  +  U1^  +  Bu*ST  "*  cn’Y  +  ^ 

pf'u'3 

-  v*1  w.d  *  h'*’-—* 


where, 
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(90) 


yP 1 

bii'  =  (i-np_  p  ') 

8o  g 


r  T  =  — 

CII 


2 

”<P  '-p  ')  +  (y-x)p  *p  ' (W'  +^-r  + V"  Y} 


p*{(l-v  )(l-np  p  ')> 

g  s  gQ  g 


E',  D^',  and  D^’  are  given  by  equations  (77),  (78),  and  (79),  respec¬ 
tively.  The  equation  of  state  for  the  gas  is  also  same  as  equation 
(80). 

Characteristic  equations: 


Along  n'i  5'  characteristic,  i.e.  r  = 


U’±/V7V 


AP'+p  '  /B  '  / p'""*  AU 


g  II  g 


b  'u' 

_ n_ _ 

(l-v  )(U'±/B'/p  ') 
s  II  g 


r  _ _  p  'u'  B  'u'  1 

+  C'+/BTT/p  *  - - r  + - V,  'At' 

II  II  g  (1-Vg)  (1_v  )(u»+/b  */p  ')  ds 

s  II  g 


-  D  'E'h (T-T  . ) At ' 
2  i  w,i 


+  v[I'»,taTj'Wj^7lt'  <92) 


Along  a  pai  ..icle  path,  i.e.  ~r  =  jyr  » 
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"  (Y-I)p  ' 

Ap/  =  G*  AP'  +  Hn'  vd  '  At'  +  D0’  h,  (T-T  .  ,)At' 


2  L(l-v8)yP‘J  “iVi  w,i; 


(Y-1)P  ’  1  P  ’U'3 


T)  1 - U - 

1  (1-v  )  yP ' 

L  S 


where 


(T-l)0s'p  ,(hg'  "  ~l  '  "  U2  } 

HXl'  _  (1-v  )  y'p ' 

s 


Boundary  Layer  Equation 

0 

Using  the  non-dimensional  momentum  thickness,  0'  =  —  and  all 

K 

the  non-dimenslonalized  parameters  listed  in  (66),  the  boundary  layer 


equation  (60)  becomes; 


U'  30'  A  ^t  Uf  B  U • 

H~  +  (1+B)H(r}  (‘p  U  '  p  “rB~ 

R  0  t 

°0 


n+B,Q.  +  JLL!!ll  +  L.iul  +  H±2iill 

(1+B;0  [_Pf'  3t’  pf'H  3x'  U'  3t'  H  3x' 


+  (i«)C  o'  r~i-  +  L.  SSL  +  i!tL 

U  p  'U*  3x'  U'  3 1 '  3x' 


where 


n  -  ai<1+B> 


O'  =  0' 


O’  =  o’ 
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Now  U',  p  '  and  P'  are  taken  from  the  one-dimensional  solution 
g 

and  p^'  is  obtained  from  the  non-dimensional  form  of  equation  (56) 


( _ _ v  , 

''T+T  /Pg 
w,i 


l+np„ 
gD  g 


*  ( — — - i)] 

VT+T  ,  ’ 

w,l  J 
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SECTION  IV 


SOLUTION  PROCEDURE 


Solution  of  Interior  Points 


As  there  is  no  analytical  solution  to  the  set  of  coupled  non¬ 
linear  partial  differential  equations  derived  in  the  previous  chapter, 
numerical  techniques  have  been  used  to  solve  the  conservation  equations 
along  with  the  boundary  layer  momentum  integral  equation  and  the  equation 
of  state  of  the  combustion  gas.  The  differential  equations  are  written 
in  finite  difference  form  and  MacCormack’s  version  [l6]  of  Lax-Wendroff 
two  step  method  [l7]  is  followed.  The  procedure  is  shown  by  an  example 
below: 


where  c  is  a  constant. 


Equation  (3.1)  can  be  written  as, 


n  At  ,  n  n, 

uj  -  c  s  (uj+i  -  V 


n  At  /-n-H  -n+1. 

uj ' c  3  <uj  '  Vi* 


and  finally. 


,  -n+1  ,  =n+l 

4  u.  +Uj 
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where  un+d  and  u?+d  are  the  first  and  second  estimated  value  of  u?+\ 
j  J  1 

The  integer  n,  and  j  denote  the  time  and  axial  position  of  a  nodal  point 

^  v  n  |  ^ 

shown  in  Figure  2.  If  c  is  a  variable,  i.e.  c(u),  c(u)^  is  used  in 

^  |  i 

the  second  »  ^imation  of  u^.  .  It  can  be  recognized  that  both  forward 

and  backward  space  derivatives  have  been  taken  into  account. 

Case  I 

Using  the  MacCormack  scheme,  the  non-dimens ionalized  conservation 
equations,  i.e.  (71)  through  (74)  can  be  written  as: 


-  n+1  n  tI,n  At'  ,  n  nx  n  At'  /lt,n  Tt,n. 
v  =  v  -  U'  .  — r  (v  -  v  )  -  v  — r  (U -  U.) 

s.  s,  j  AxX  s. ..  s.  s.  Ax^  j+1  j 


j+1  j 


pd  r 

L  sj. 


(100) 


•n+1  =,  n  'n  -  U'n^-  (P  'n  -  n  -n 


j+1  b  j; 


J  ■  ~r  (U’“  U'J)  + - fv,  ']n  At' 

-V  n)  Ax0  j+1  j  (1-v  n)  ds  J 


(p  '-p  '  ) 
S  8  . 


(101) 


U'"+l  =  U'"  -  U'J  t^t  (U'“  -  U'") 

j  j  J  AXq  j+1  1 


—Tit  7^-  (P’“  -  P’n)  -  ULF  n  At’ 

p»  Axi  j+1  J  J 


(102) 


50 


Piston  Path 


Figure  2.  Numerical  Scheme  for  Interior  Points 


characteristic 


Figure  3.  Scheme  for  End  Points 


Thus,  for  al]  the  nodal  points  (except  the  last  point  adjacent  to  the 
piston  base)  the  first  estimation  regarding  the  ballistic  properties 
after  a  time  increment  At'  is  made  by  knowing  the  present  values  of  the 
properties  at  the  point  of  interest  and  at  its  forward  nodal  point.  The 
second  estimation  is  done  as  follows: 
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“  n+1  _  n  .~,n+l  At1  n+1  ~  n+1. 
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(109) 


U'n+1  =  u'n  -  U'n+1  Ail  (£},n+l  _  u'n+1) 

J  J  j  Ax'  v  j  j-l' 
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and  finally, 


n+1 


~  n+1  s  n+1 
v  +  v 

u  SJ  SJ  . 
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p’n+i  -  is  p'n+1  +  p“'n+1 

gj  L  8j  gj 


u'n+1  =  iJu'n+1  +  fj'n+1 

J  j  j 

p,n+l  ,  , +  j,n+lj 


(113) 


(114) 


(115) 


From  the  equation  of  state  (80):' 


,n+l  (1-np  P'n+1) 

,,n+l  a  ij__  8o  gi 

j  *  ,n+l  (1-np  ) 

gO 

gj 


(116) 


Similarly,  to  calculate  the  momentum  thickness  at  a  nodal  point  aftei 
increment  At',  equation  (95)  is  written  as: 
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where 
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taking 


T,.  m 

yf  =  Mg  (F> 

°0  O 


(119) 


m  being  a  suitable  constant. 

The  densities  p^n  and  p^n  are  calculated  from  equation  (97) 

j  j 

with  the  assumption. 
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finally. 


The  same  procedure  is  followed  for  Case  II  starting  with  appro¬ 
priate  conservation  equations,  namely  equations  (86)  through  (89), 
same  equation  of  state  (80)  and  boundary  layer  momentum  equation  (95). 
The  only  points  of  differences  are:  (1)  no  solid  particles  beyond  1/ 
and  (2)  the  burning  rate  r^  is  chosen  corresponding  to  the  average 
pressure  in  the  space  between  the  tube  head  end  and  Lq. 

Solution  for  Boundary  Points 

It  has  been  stated  earlier  that  to  calculate  the  ballistic  proper¬ 
ties  at  the  tube  head  end  and  at  the  piston  base  end,  one  needs  the 
characteristic  equations.  Typical  characteristic  directions  are  shown 
in  Figure  3.  Let,  at  any  time  t',  the  piston  be  at  position  1  with 
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velocity  Its  position  and  velocity  after  time  At'  can  be  calcu¬ 

lated  by  using  equations  (18)  and  (19)  as  follows: 


U'  0  =  U'  +  (a  At) 
p,2  p,l  Uq  p 


and 


U  At  +  h  a  (At) * 
P,1  P 


(123) 


where, 

P  A 

a  =  — EaS — 2.  (p  mean  0f  p  and  P„) 

p  M  p  12 

P 

Then  the  n '-characteristic  is  traced  back  using  the  appropriate  expres¬ 
sion: 

for  Case  I:  Ax'  »  (U’WbT/p7!  At' 

I  m 

(124) 

for  Case  II:  Ax'  -  (U'+ZB^/p^)  At' 


The  point  X  is  thus  determined  and  all  the  properties  are  interpolated 
between  the  nodal  points  in  each  side.  Pressure  at  point  2  is  calculated 
by  applying  equation  (81)  for  Case  I,  and  equation  (92)  for  Case  II, 

fir 
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for  Case  IT,  Pi  =  P'  -  p»  ✓fT/p*  (U *  0-U') 
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The  gas  density  and  volume  fraction  of  solids  at  the  new  base  point  2 
are  determined  from  characteristic  equations  along  a  particle  path  as 
follows : 

for  Case  I  and  II, 


p'  =  p-  +  g'(p;-p!)  +  h;  v  *  At' 

g,2  g,l  2  1  1,11  dg 


(127) 


for  Case  I, 


v  (l-v  )  p* 

v  0  =  v  +  - -  (p'  -p*  )  -  — j-  v  ’  At1 

S>2  S|1  Dg  S>2g,l  Pg  d^ 


(128) 


for  Case  II, 


v  o  =  0 
s,2 


(129) 
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For  both  points  1  and  2,  the  momentum  thickness  is  zero,  which  implies 
that  both  friction  factor  and  film  heat  transfer  coefficient  at  point 

1  and  2  are  infinitely  large.  Therefore,  the  last  two  terms  of 

equations  (82)  and  (93)  have  been  deleted  while  writing  the 

equation  (127)  For  the  same  reason,  in  equations  (125)  and  (126) 

the  values  for  h , (T-T 

i  w,i 

nodal  point  adjacent  to  the  first  base  point  1.  All  the  coefficients 
used  in  equations  (124)  through  (128)  are  mean  values  between  point 

2  and  X  or  point  2  and  1  depending  on  the  characteristic  used.  The 
properties  at  point  2  are  first  assumed  to  be  the  same  as  point  1  and 
then  iteration  is  carried  on  until  the  values  converge  within  the 
specified  limit. 

For  the  tube  head  end, 

a  UJ*  =  0  (130) 

and  the  momentum  thickness  and  heat  transfer  coefficient  are  also  zero. 
The  ^'-characteristic  is  traced  by  using, 

for  Case  I:  Ax'  «  (U'-Zo'/p')  At' 

1  “  (131) 

for  Case  II:  Ax'  “  (U’-Zb^/p^)  At' 

By  knowing  the  properties  at  X’  and  1',  properties  at  point  2'  are 
obtained  in  the  following  manner: 


and 


(Ree)B 


are  taken  corresponding  to  the 
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Case  I 
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and  from  (2.86), 
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The  expression  for  p^2,  is  obtained  by  replacing  Hj  by  in  equation 
(133). 
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The  iteration  procedure  for  the  tube  head  end  is  the  same  as 
that  for  the  piston  base  end  stated  earlier. 

The  properties  at  the  nodal  point (s)  adjacent  to  the  base  point 
(shown  by  *  in  Figure  3),  which  cannot  be  calculated  from  the  Lax-Wendroff 
method, are  determined  by  linear  interpolation  between  the  piston  base 
point  and  the  nearest  point  where  properties  have  been  calculated  from 
the  Lax-Wendroff  method. 

It  is  noted  that  the  spatial  interval  Ax^  is  fixed  for  the  entire 

solution  and  can  be  chosen  arbitrarily  depending  upon  the  desired 

accuracy.  But,  for  the  stability  of  the  Law-Wendroff  solution,  the  time 

At' 

interval  At'  must  be  chosen  such  that  nowhere  exceeds  the  slope  of 

o 

any  characteristic  [l7] .  This  implies  that  at  every  time  step, 


Ax' 

At'  £ - •°"  ■  -  -  for  Case  I 

|u* |+/b;/p' 
l  m 

Ax' 

At'  <_ - °  ■■  ■  for  Case  II 

,u' ,+/®Ti7pI 


(137) 


Therefore,  before  selecting  a  new  time  interval,  the  right  hand  side  of 
(137)  is  calculated  at  each  nodal  point  (including  the  end  points)  and 
then  the  lowest  value  is  chosen  as  the  next  time  step. 

Determination  of  Wall  Temperature 
The  differential  equation  (64)  in  finite  difference  form  can 
be  written  as  (see  Figure  4  for  notations) : 
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t 


Figure  4.  Numerical  Scheme  for  Determination  of  Tube  Wall  Temperature. 


Figure  5.  Heat  Balance  for  a  Thin  Circular  Element  at  the  Inside 
Surface  of  the  Tube 
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Therefore,  the  temperature,  at  any  interior  point  in  the  tube  wall  after 
a  time  interval  At  can  be  calculated  from  the  knowledge  of  present  tem¬ 
peratures  at  and  around  the  point  of  interest.  For  boundary  points, 
however,  a  heat  balance  as  described  below  is  required: 

Inner  Surface:  With  reference  to  the  Figure  5: 
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In  finite  difference  form, 
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Similarly,  for  the  outer  surface, 
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The  stability  conditions  [l8,  19]  for  the  equations  (138)  through  (140) 
are,  respectively: 
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Very  small  values  of  Ar  (0.05  millimeter)  are  taken,  and  the  selected 
At  is  the  least  of  the  values  calculated  from  the  right-hand  3ide  of 
(141),  (142),  and  (143). 

For  a  thick  wall  and  initially  cold  tube,  the  temperature  wave 
in  a  single  shot  does  not  generally  reach  the  outer  surface  and,  there¬ 
fore,  equation  (140)  can  be  disregarded. 

Summary  of  the  Procedure 

Once  the  piston-cylinder  arrangement,  the  initial  conditions  and 
all  other  input  parameters  are  chosen,  the  solution  proceeds  according 
to  the  following  steps: 

(L)  The  time  interval  At*  is  determined  in  accordance  with 
expression  (3.40)  and  the  burning  rate  is  taken  corresponding  to  the 
average  burning  pressure. 

(2)  The  new  piston  position  and  its  velocity  are  calculated, 
and  using  the  appropriate  characteristic  equations  as  indicated  earlier 
the  new  ballistic  properties  at  both  the  piston  base  end  and  the  tube 
head  end  are  determined. 

(3)  The  interior  points  are  solved  either  by  the  Lax-Wendroff 
method  or  by  linear  interpolation  as  discussed  earlier. 

(4)  The  new  heat  transfer  coefficient  is  determined  from  the  new 
ballistic  properties  and  the  momentum  thickness  at  all  the  nodal  points. 
The  new  wall  temperature  is  also  calculated  using  the  mean  heat  transfer 
coefficient. 

(5)  All  the  calculated  values  are  stored  as  the  present  values 
and  reused  for  the  next  time  step. 
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Thus  the  solution  proceeds  until  the  piston  reaches  the  desired 
position.  A  computer  program  for  the  entire  solution  procedure  was 
written  in  FORTRAN  V  and  was  run  to  obtain  all  of  the  results  presented 
in  the  following  chapter.  The  flow  chart  for  the  program  has  been  shown 
in  Appendix  III.  The  computation  time  is  approximately  four  minutes  for 
the  typical  cases  run  on  the  Georgia  Institute  of  Technology's  UNIVAC  1108 


machine . 


SECTION  V 
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Covolume,  n 

Specific  heat  at  constant 
pressure,  c 

P 

Viscosity  (at  3000°K) ,  p 

8o 

Thermal  conductivity  (at 
3000°K) ,  kt 
go 


0.00095  m3/kg. 

0.412  kcal/kg-°K 
0.00007  kg/m-sec 

0.000034  kcal/m-sec-°K 


Tube  material  properties; 


Thermal  diffusivity,  a 


0.126  cm  /sec 


Thermal  conductivity,  * 


0.0138  kcal/m-sec-  K 


Initial  tube  temperature, 
Tamb 


300°K 


The  initial  gas  density  p  as  calculated  from  equ:  *■  (2.15)  is 

3  8° 

19.14  kg/m  and  the  potential  of  the  propellant  is: 


W  "  "  985  • 66  kcal/kg 


The  burning  rate  versus  pressure  data  for  the  propellant  has 
been  taken  from  reference  [20]  and  is  presented  in  Table  I.  Linear 
interpolation  is  used  to  determine  the  burning  rate  at  the  desired 
pressure.  To  ensure  the  convergence  of  the  solution,  a  single  iteration 
on  the  burning  rate  is  performed  in  each  time  step  as  shown  in  the  flow 
diagram  in  Appendix  III. 

For  Case  1,  the  solid  particles  are  initially  assumed  to  be 
evenly  distributed  in  the  chamber.  But,  in  Case  II,  a  specific  initial 
distribution,  namely  a  constant  value  up  to  the  second  nodal  point  from 
the  piston  and  then  linearly  to  zero  at  the  piston  base,  is  chosen  to 
avoid  the  discontinuity  at  x equal  to  Lq.  This  has  been  shown  in  Figure 
6. 


68 


TABLE  I.  PRESSURE  VERSUS  BURNING  RATE  DATA  FOR  THE  PROPELLANT  [20] 


•■5 

Pressure  X  10 

2 

Newton/m 

Burning  rate, 

m/sec 

20.68 

0.00330 

34.46 

0.00508 

48.25 

0.00711 

68.93 

0.00965 

103.39 

0.01320 

137.86 

0.01727 

172.32 

0.02057 

2C6.78 

0.02438 

275.71 

0.03048 

344.64 

0.03683 

413.57 

0.04369 

551.43 

0.05588 

689.28 

0.06858 

1378.57 

0.11684 

2067. 86 

0.17018 

2757.14 

0.21082 

3446.43 

0.24384 

4825.00 

0.3098 f, 

6892.86 

0.40132 

13735.71 

0.63500 

*9 
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The  one-seventh  profile  has  been  used  quite  extensively  in  the  past  to 
compute  the  turbulent  boundary  layers  with  favorable  pressure  gradients 
[22,  23].  The  same  profile  is  assumed  under  the  "standard  conditions" 
and  the  corresponding  value  for  the  shape  factor,  i.e.  1.2857,  is  taken 
for  the  boundary  layei  computation. 

The  viscosity  of  the  combustion  gas  is  assumed  to  be  proportional 
to  the  square  root  of  the  absolute  temperature  which  implies  that  the 
value  of  m  in  (119)  is  0,5.  The  same  relation  is  assumed  between  the 
gas  conductivity  and  its  absolute  temperature.  These  yield  a  constant 
value  of  0.8482  for  the  Prandtl  number  of  the  gas. 

One-Dimensional  Solution 

Case  I.  The  results  of  the  one-dimensional  analysis  have  been 
presented  in  Figures  7  through  14.  Comparison  with  the  solution  neglecting 
the  heat  transfer  and  skin  friction  shows  insignificant  effect  of  these 
phenomena  on  the  ballistic  properties  of  the  piston-cylinder  arrangement. 
But  the  following  observations  can  be  made  from  these  results: 

1)  The  Lagrange  approximation  of  lJnear  velocity  distribution 
and  constant  gas  density  is  not  a  good  approximation  of  the  real  situation. 
It  can  be  noted  from  Figure  10  that  a  considerable  amount  of  time 
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Figure  7.  Piston  Path  for  Case  I. 
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Figure  10.  Spacewise  Distribution  of  Velocity  at  Various 
Times  (Case  1) 
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Figure  12.  Spacewise  Distribution  of  Gas  Temperature  at  Various 
Times  (Case  I) 
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(approximately  one-third  of  the  total  time)  is  elapsed  before  the 
velocity  distribution  along  the  length  of  the  tube  becomes  linear. 

The  gas  density,  on  Lhe  other  hand,  has  always  a  drooping  character¬ 
istic  from  the  breech  end  to  the  piston  base  end  (Figure  9  and  13). 

The  same  characteristic  is  observed  for  the  volume  fraction  of  solids, 
(Figure  14). 

(2)  The  gas  pressure  varies  considerably  along  the  length  of 

the  tube,  and  at  the  peak  pressure,  the  difference  between  the  pressures 
at  the  two  ends  is  as  high  as  20  per  cent  of  the  breech  pressure 
(Figures  8  and  11).  The  heat  transfer  to  the  tube  wall  reduces  the 
pressure  at  all  points,  whereas  the  skin  friction  reduces  the  piston 
base  pressure  as  the  piston  reaches  the  end  of  the  tube.  These  two 
effects  together  reduce  the  final  piston  velocity  to  some  extent. 

(3)  The  gas  temperatures  at  the  breech  and  the  piston  base  are 

of  the  same  value  all  the  time  except  for  a  short  period  in  the  beginning 
(Figure  9).  There  is,  however,  a  sag  in  between  the  two  end  points  due 
to  the  heat  loss  to  the  tube  wall  (Figure  12). 

Case  II.  Similar  results  for  Case  II  have  been  presented  in 
Figures  15  through  21.  A  comparison  with  Case  I  reveals  that  the  final 
values  of  piston  velocity  and  total  time  of  travel  do  not  differ  much 
from  those  in  Case  I.  But  the  peak  breech  p-.essure  can  be  10  to  15  per¬ 
cent  higher  than  the  corresponding  pressuie  in  Case  I.  This  causes  the 
propellant  to  burn  faster.  The  piston  base  pressure,  however,  remains 
very  close  to  the  corresponding  pressure  in  Case  I  except  for  a  short 
period  towards  the  end.  This  accounts  for  the  slight  variation  in  the 
final  piston  velocity  between  the  two  cases. 
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Figure  17.  Spacewise  Distribution  of  Gas  Velocity  at  Various 
Tiues  (Case  II) 
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Figure  20,  Spacewise  Distribution  of  Gas  Density  at  Various 
Timas  (Case  II) 


A  couple  of  interesting  phenomena  are  observed  in  this  case  of 
stationary  solids: 

(1)  Just  after  the  propellant  is  completely  burnt,  the  breech 
pressure  falls  rapidly  and  finally  becomes  less  than  the  piston  base 
pressure.  The  same  phenomenon  was  also  observed  by  Carriere  [5].  This 
happens  because  all  the  solid  particles  are  assumed  to  stay  near  the 
breech  end  all  the  time,  whereas  there  are  no  solids  at  the  piston  base. 
Therefore,  when  the  solids  are  completely  burnt,  no  sudden  change  in  the 
pressure  slope  occurs  at  the  piston  base  as  it  is  observed  at  the  breech. 

(2)  Due  to  the  rapid  change  in  the  volume  fraction  of  solids  v 

s 

near  x  equal  to  LQ,  a  pressure  difference  sufficient  to  produce  a  local 
gas  velocity  higher  than  the  piston  velocity  is  created  (Figure  17). 

The  gas  having  higher  velocity  slams  at  the  back  of  the  piston  and  thus 
increases  the  temperature  (Figure  19) .  This  large  temperature  rise  is 
confined  within  a  thin  layer  at  the  piston  base  and  does  not  affect  the 
rest  of  the  gas.  The  oscillations  observed  in  Figures  19  and  20  are  not 
due  to  the  numerical  instability,  but  are  most  probably  due  to  the  sudden 
area  change  near  x  equal  to  Lq. 

Boundary  Layer  and  Heat  Transfer  Solution 

The  results  showing  the  boundary  layer  thickness,  heat  transfer 
coefficient,  and  the  wall  temperature  for  moving  solids  (i.e.,  Case  I) 
are  presented  in  Figures  22  through  27.  As  the  time  increases,  the 
boundary  layer  thickness  increases  to  a  maximum  value  of  approximately 
20  percent  of  the  tube  radius  when  the  piston  reaches  the  end  of  the 
tube.  This  implies  that  the  maximum  value  of  the  displacement  thickness 
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Figure  24.  Spacewise  Distribution  of  Wall  Surface  Temperature  at 
Various  Times  (Case  i) 
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Figure  25.  Variation  of  Wall  Temperature  at  a  Particular  Position 
with  Depth  at  Various  Times  (Case  I) 
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Flgura  26. 


Variation  of  Wall  Surface  Temperatures  at  Certain  Fixed  Locations 
with  Time  (Case  I) 


is  less  than  three  percent  of  the  tube  radius  and,  therefore,  the 
assumption  of  a  thin  boundary  layer  is  valid. 

The  order  of  magnitude  of  the  heat  transfer  coefficient  is 
extremely  high  due  to  high  gas  density  and  velocity  (Figures  23  and 

26) .  The  surface  temperature  of  the  tube  wall  reaches  as  high  as 
1100°K,  and  it  occurs  at  the  initial  piston  position  (Figures  24  and 

27) .  The  total  time  is  so  short  that  in  spite  of  a  very  steep  radial 
temperature  gradient  and  high  thermal  diffusivity  of  the  tube  material, 
the  temperature  T  ave  cannot  penetrate  more  than  one  millimeter  into 
the  tube  wall  (Figure  25).  This  justifies  the  exclusion  of  equation 
(140)  from  the  computer  program. 

The  heat  flux,  h^(T-Tw  ^ )  at  certain  fixed  positions  along  the 

length  of  the  tube  are  shown  in  Figure  28.  Although  the  maximum  value 

2 

of  heat  flux  could  be  as  high  as  350,000  kcal/m  -sec,  these  types  of 
fantastically  high  values  last  only  for  one  or  two  microseconds.  The 
average  value  of  heat  flux  would  be  around  50,000  kcal/m  -sec. 

The  same  type  of  results  was  also  obtained  for  Case  II,  and  the 
total  heat  losses  for  both  the  cases  are  compared  in  Figure  29.  It 
shows  that  the  heat  loss  in  Case  II  is  about  ten  percent  higher  than 
that  in  Case  I.  This  is  mainly  due  to  the  higher  gas  velocity  in  the 
initial  period  of  Case  II.  After  this  Initial  period,  the  heat  transfer 
coefficients  in  the  two  cases  are  quite  close. 

lhe  important  results  of  both  the  cases  are  tabulated  in  Table 
II.  A  mass  and  energy  balance  detailed  in  Table  III  shows  that  the  compu- 
tati  i  error  is  less  than  0.5  per  cent.  The  ballistic  efficiency  (i.e.  , 
the  ratio  of  final  kinetic  energy  of  the  piston  and  the  propellant  energy) 
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Variation  of  Heat  Fluxes  at  Certain  Fixed  Locations  with  Time  (Case  I) 


Comparison  of  Total  Heat  Losses  to  the  Tube  Wall  in  Case  I  and  Case  II 


TABLE  II.  COMPARISON  OF  RESULTS  FOR  TOO  LIMITING  CASES  OF  SOLIDS 
VELOCITY 


Case  I 

(solids  moving) 

Case  II 

(solids  stationary) 

Time  of  travel 
(second) 

0.002465 

0.002448 

Final  piston  velocity 
(m/sec) 

1242.3 

1271.3 

Peak  breech  pressure 
(atmosphere) 

6450 

7200 

Peak  surface  temperature 
(°K) 

1067 

1270 

Ballistic  efficiency 
(%) 

35.47 

37.14 

Total  heat  loss 
(kcal) 

8.64 

9.60 

Keat  loss  in  percentage 
of  input  energy  (%) 

5.10 

5.66 
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TABLE  III.  MASS  AND  ENERGY  BALANCE  FOR  CASE  I  and  CASE  II 


Case  I 

Case  II 

Initial  Conditions: 

Propellant  charge 
(kg) 

0.172 

0.17? 

Propellant  energy 
(kcal) 

165.53 

169.53 

Final  Conditions: 

Total  gas  mass  (kg) 

0.1716 

0.1713 

Gas  internal  energy 
(kcal) 

90.10 

87.63 

Gas  kinetic  energy 
(kcal) 

9.90 

3.56 

Piston  kinetic  energy 
(kcal) 

60.13 

62.96 

Heat  loss  (kcal) 

8.64 

9.60 

Total  energy  (kcal) 

168.77 

168.75 

Error  in  Mass  Balance  (%) 

-0.23? 

-0.407 

Error  in  Energy  Balance  (%) 

-0.450 

-0.460 
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has  also  been  presented.  The  total  heat  loss  to  the  tube  wall  is  found 
to  be  five  to  six  percent  of  the  input  energy  and  about  15  percent 
of  the  final  piston  kinetic  energy. 

Parameter  Variation 

Because  of  the  large  number  of  independent  design  parameters, 
no  general  correlation  is  attempted  here.  Only  a  few  important 
parameters  are  varied  for  Case  I  to  study  their  effect  on  the  ballis¬ 
tic  as  well  as  heat  transfer  solution,  and  the  important  results  are 
presented  in  Table  IV. 

Initial  Chamber  Pressure,  P 

The  selection  of  initial  chamber  pressure,  i.e.,  piston  start 
pressure  is  quite  arbitrary  as  it  is  very  difficult  in  practice  to 
determine  the  exact  pressure  at  which  the  piston  starts  to  move. 

Therefore,  two  different  initial  chamber  pressures  (100  atmosphere 
and  300  atmosphere)  other  than  the  standard  200  atmosphere  are 
considered, and  the  solution  neglecting  the  heat  transfer  and  skin 
friction  is  presented  in  Figure  30.  The  peak  breech  pressure  and  the 
final  piston  velocity  are  very  much  the  same  for  all  the  three  cases 
of  different  piston  start  pressures.  Only  the  time  of  travel  is  prolonged 
as  the  initial  pressure  decreases.  It  is  obvious  that  the  heat  transfer 
solution  would  be  very  close  for  all  the  three  initial  pressures  because 
of  close  ballistic  properties.  This  implies  that  the  piston  start  pres¬ 
sure  has  insignificant  effect  on  the  overall  performance  of  the  device. 
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Profile  Shape  Factor,  H 


In  this  effort  the  profile  shape  factor,  H,  is  assumed  to 
be  a  constant  throughout  the  entire  solution;  however,  this  may  not 
be  true  in  the  real  situation.  For  non-uniform  steady  flow,  a  favor¬ 
able  pressure  gradient  lowers  the  value  of  H  [?4]  whereas  in  uniform 
steady  flow  an  increase  in  flow  Mach  number  increases  H  [l3j.  In  the 
present  case,  however,  it  is  difficult  to  predict  its  probable  variation. 
Therefore,  two  different  values  of  shape  factor  (1.4  and  1.2222), 
corresponding  to  the  one--fifth  and  the  one-ninth  velocity  profile,  are 
taken  and  the  results  are  compared  with  those  for  H  equal  to  1.2857. 

It  is  clear  from  Figure  31  that  lower  values  of  H  cause  higher  values 
of  boundary  layer  thickness  and  heat  transfer  coefficient,  i.e.  the 
viscous  effect  of  the  fluid  is  higher.  Because  of  this,  the  total 
heat  loss  to  the  tub  wall  and  the  peak  inner  surface  temperature  go  up 
as  the  shape  factor  decreases.  However,  at  H  equal  to  1.2222,  these 
values  do  not  exceed  the  corresponding  values  for  the  standard  conditions 
by  more  than  ten  percent.  The  total  heat  losses  to  the  tube  wall  for 
the  three  different  values  of  shape  factor  are  compared  in  Figure  32. 

Tube  Inside  Diameter,  D 

The  tube  inside  diameter  is  varied  keeping  the  loading  density 

(m  /V  ) ,  and  the  piston  mass  per  unit  area  (M  /A  )  constant.  Two  tube 
'  si  o  ’  r  p  p 

diameters  (2cm  and  4cm)  are  selected  with  appropriate  propellant 
charge  mg  and  piston  mass  M^.  The  results  are  shown  in  Figure  33.  It 
can  be  noted  from  Table  IV  that,  although  an  increase  in  tube  diameter 
means  an  increase  in  the  total  heat  loss,  the  heat  loss  per  unit  input 
energy  decreases  with  larger  tube  diameter.  This  is  due  to  decrease  in 
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Figure  31.  Comparison  of  Boundary  Layer  Thickness  and  Heat 
Transfer  Coefficient  for  Various  Profile  Shape 
Factors  at  0.002337  sec  (Case  I) 


(kcal/m  -sec-  K) 


Figure  32.  Comparison  of  Total  Heat  Loss  to  the  Tube  Wall  for  Various  Profile  Shape 
Factors  (Case  I) 
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Figure  33.  Comparison  of  Breech  Pressure,  Piston  Velocity, and  Total  Heat  Loss  to  the  Tube  Wal 
for  Various  Tube  Diameters  wl.th  Same  Loading  Density  and  Same  Piston  Mass  per  Unit 

Area  (Case  I) 
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the  surface  Lu  volume  ratio.  In  other  words ,  for  ballistically  similar 
devices,  the  increase  in  tube  diameter  reduces  the  heat  loss  per  unit 
mass  of  gas  and  this  accounts  for  the  slightly  better  ballistic  results 
obtained  for  the  4-cm-diameter  tube. 

Propellant  Charge,  mg . 

Two  different  values  of  the  propellant  charge  (0.15  kg  and 
0.19  kg)  are  chosen  apart  from  the  standard  value  of  0.172  kg,  and  the 
results  are  presented  in  Figure  34.  It  is  obvious  that  an  increase 
in  propellant  charge  improves  the  ballastic  efficiency  of  the  device. 

At  the  same  time  this  increases  the  peak  pressure,  heat  transfer 
coefficient,  and  Hie  peak  surface  temperature  which  put  a  limit  on  the 
propellant  charge. 

Piston  Mass,  M 

- ; - 1 — p 

Piston  mass  plays  an  important  role  in  the  problem  of  internal 
ballistics.  Therefore,  besides  the  standard  mass  of  0.326  kg,  two  other 
pistons  having  masses  equal  to  0.2  kg  and  0.5  kg  are  considered.  The 
results  are  shown  in  Figure  35.  The  heavier  the  piston,  the  slower  it 
moves,  thereby  leaving  less  room  for  the  combustion  gas  to  expand  which 
causes  an  increase  in  the  peak  pressure.  Although  the  heavier  piston 
moves  slower,  the  ballistic  efficiency  of  the  device  is  improved  and 
therefore  suitable  for  the  application  where  energy  conversion  is  of 
prime  interest.  If  higher  velocity  is  desired,  a  lighter  piston  would 
be  chosen.  The  surface  temperature  is  also  lower  in  the  case  of  a 
lighter  piston  due  to  lower  pressures  which  lead  to  lower  heat  transfer 
coefficients. 
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Figure  34.  Comparison  of  Breech  Pressure,  Piston  Velocity  and  Total  Heat  Loss  to  the  Tube  Wall 
for  Various  Propellant  Charges  (Case  I) 


Figure  35.  Comparison  of  Breech  Pressure,  Piston  Velocity  and  Total  Heat  Loss  to  the  Tube  Wall 
for  Various  Piston  Masses  (Case  I) 
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Web  Thickness,  w 

Although  the  web  thickness  is  only  a  geometrical  property  of  the 
solid  particles,  it  is  important  because  it  determines  the  total 
burning  surface  for  a  given  propellant  charge.  For  thinner  webs  more 
surface  i.  available  for  burning  and  consequently  the  pressure  rise  is 
more  rapid.  This  aspect  is  clear  from  Figure  36  where  results  of  three 
different  web  thicknesses  (0.5mm,  0711  mm  (standard),  and  0.9  mm; 
are  presented.  A  rapid  pressure  rise  naturally  accelerates  the  piston 
faster  and  thus  improves  the  ballistic  efficiency.  This  gain  is 
neutralized  by  a  much  higher  peak  pressure  and  wall  surface  temperature. 

Comparison  with  Other  Work 

No  analytical  work  in  the  past  considered  the  movement  of  the 
solid  particles  in  the  one-dimensional  ballistic  solution.  Although 
Carriere  [5]  studied  the  problem  of  internal  ballistics  assuming  the 
solids  to  be  stationary,  it  was  not  possible  to  determine  the  input  data 
and  final  results  from  his  publication.  Therefore,  a  quantitative 
comparison  could  not  be  made.  However,  an  excellent  qualitative  agree¬ 
ment  is  observed  between  his  results  showing  the  piston  path,  piston 
velocity  and  end  pressures,  and  the  results  obtained  from  the  present 
analysis  for  the  case  of  stationary  solids.  Unfortunately,  no  work 
presently  exists  that  the  rpacewise  distribution  of  the  ballistic  properties; 
therefore, no  comparison  can  be  made. 

The  boundary  layer  and  heat  transfer  analysis  of  Hicks  and 
Thornhill  [2]  assumed  the  Lagrange  approximation  and  omitted  the  zero 
boundary  layer  thickness  condition  at  the  piston  base.  The  values  of 
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Figure  36.  Comparison  of  Breech  Pressure,  Piston  Velocity  and  Total  Heat  Loss  to  the 
Tube  Wall  for  Various  Web  Thicknesses  (Case  I) 


TABLE  IV.  RESULTS  FOR  VARIOUS  INPUT  PARAMETERS 


P 


! 

I 


III 


heat  transfer  coefficient  for  a  typical  case  (input  data  not  indicated) , 

as  presented  in  reference  2,  are  lower  than  the  values  obtained  for  the 

typical  case  in  the  present  study,  by  a  factor  of  approximately  two. 

The  reason  could  be  due  to  entirely  different  input  data  and  different 

boundary  layer  thickness  condition  at  the  piston  base.  However,  the 

value  of  peak  surface  temperature  and  the  location  when  it  occurs  are 

in  good  agreement  with  the  study  of  Hicks  and  Thornhill. 

A  very  recent  analysis  on  convective  heat  transfer  in  gun  barrels 

[25],  which  is  also  based  on  the  Lagrange  approximation,  indicates  that 

the  heat  flux  at  the  inside  surface  of  the  barrel  can  be  as  high  as 
5  2 

2.7x10  kcal/m  -sec.  This  value  is  quite  close  co  the  expected  maxima 
shown  in  Figure  28  for  the  typical  case. 
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Efelin  Condition 

In  order  to  compare  results  with  experimental  temperature  data 
gathered  at  Eglin  Air  Force  Base,  Florida,  conditions  other  than  the 
Standard  conditions  were  run  with  the  following  input  conditions: 


Tube  length,  Lt 

Tube  inside  diameter,  D 

PLston  Mass,  M 
P 

Initial  conditions: 

PLston  position,  Lq 

Chamber  pressure  (Piston 
start  pressure) ,  Pq 

Gas  temperature  (explosion 
temperature,  tq 

Change  of  propellant,  mg^ 

Propellant  properties: 

Density,  P 

s 

Initial  web  thickness,  w 


Gas  Properties: 

Molecular  weight,  M 
Ratio  of  specific  heats,  y 
covolume ,  n 

Specific  heat  at  constant 
Pressure,  c 

P 

Viscosity  (at  3000°K) ,  y 

go 

Thermal  conductivity  (at 
3000°K) ,  k 

go 


2m 
3  cm 
.326  kg 


25  cm 

200  atm 

3000°K 
0.140  kg 


1670  kg/m3 
0.7366  mm 

M-10,  single 
perforated 


24 

1,252 

0.00095  m3/kg 

0.412  kcal/kg°K 
0.00007  kg/rc/sec 

0.000034  kcal./m/sec°K 
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Tube  material  properties: 

2 

Thermal  diffusivity,  a  0.062  cm  /sec 

w 

Thermal  conductivity,^  0.0138  xcal/m-bec°K 

Initial  tube  temperature,  Tam^  300°K 

The  burning  rate  was  again  taken  from  Table  I,  and  only  Case  I 
was  taken,  i.e.,  the  soil  ids  move  with  the  gas. 

The  results  of  this  calculation  are  shown  in  Figures  37  through 
43,  and  agree  extremely  well  with  experimental  data  for  the  peak  wall 
temperature  taken  by  Eglin  Air  Force  Base.  The  calculated  maximum  wall 
temperature  shown  in  Figure  43  is  1020°K  at  x  =  25  cm.  The  experimental 
maximum  measured  at  Eglin  was  approximately  1000°K.  No  measured  parameters 
such  as  chamber  pressures  were  input  to  the  program  to  calculate  this  wall 
temperature.  The  only  necessary  input  conditions  are  propellant,  gun, 
and  loading  parameters  listed  at  the  beginning  of  this  section.  These 
Eglin  conditions  are  obviously  more  realistic  than  the  original  standard 
conditions  due  to  the  smaller  propellant  loading  (140  grams  instead  of 
172  grams).  The  peak  head  and  chamber  pressure  is  approximately  58,000  psi, 
as  seen  from  Figure  38. 

This  excellent  agreement  verifies  the  accuracy  of  the  detailed 
ballistic  and  heat  transfer  model  developed  here  and  used  without  the 
aid  of  experimental  ballistic  data. 
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Figure  38.  End  Pressures,  Piston  Velocity,  and  Mass  Fractions  of  Solids  (Case 
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Figure  40.  Gas  Temperature  Versus  Axial  Location 
Eglin  Conditions) 


at  Various  Times  tCase  I, 
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Surface 


Wall  temperature 


SECTION  VI 


CONCLUSIONS 

The  following  conclusions  can  be  drawn  in  the  context  of  the 
i exults  presented  in  Part  I. 

1.  The  Lagrange  approximation  of  linear  velocity  distribution 
and  uniform  gas  density  along  the  length  of  the  tube  is  not  a  good 
representation  of  the  real  case.  It  takes  a  considerable  amount  of 
total  tine  before  the  velocity  distribution  can  be  linear.  Further¬ 
more,  the  gas  density  cannot  be  called  uniform  at  any  time. 

2.  There  is  a  large  pressure  gradient  along  the  length  of  the 
tube  and  at  the  peak  condition  the  difference  between  the  pressures 

at  the  two  ends  can  be  as  high  as  20  to  30  percent  of  the  maximum  pressure. 

3.  As  the  piston  moves,  the  gas  temperature  continuously  decreases 
with  almost  a  uniform  spacewise  distribution.  For  the  case  of  stationary 
solids,  however,  a  steep  spacewise  temperature  rise  is  observed  at  the 
back  of  the  piston. 

4.  While  the  final  ballistic  results  are  essentially  the  same 
for  the  fo  extreme  cases  of  solids  velocity,  the  peak  pressure  in  the 
case  of  stationary  solids  is  about  10  to  15  percent  higher  than  in  the 
case  of  moving  solids.. 

5.  The  maximum  boundary  layer  thickness  can  be  on  the  order  of 
20  percent  of  the  tube  radius  in  typical  cases,  ^or  ballistically 
similar  devices  the  ratio  of  the  maximum  boundary  layer  thickness  to  the 
tube  radius  increases  as  the  tube  diameter  is  reduced.  Therefore, 
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the  assumption  of  a  thin  boundary  layer  may  not  hold  good  for  very  small 
diameter  tubes. 

6.  The  order  of  magnitude  of  the  heat  transfer  coefficient  and 

the  heat  flux  at  the  inner  surface  of  the  tube  is  extremely  high. 

2  o  2 

Average  values  of  50  kcal/m  -sec-  K  and  50,000  kcal/m  -sec  for  the  heat 
transfer  coefficient  and  the  heat  flux,  respectively,  can  be  expected 
for  typical  cases.  The  maximum  values  can  be  five  to  six  times  higher 
than  the  average  values  but  the  maxima  do  not  last  for  more  than  a 
few  microseconds. 

7.  The  tube  inner  surface  temperature  can  reach  1000° k  for  typical 
cases,  and  it  occurs  near  the  initial  piston  position.  The  time  of  travel 
is  so  short  that, even  with  the  extremely  high  values  of  heat  fluxes  and 
high  thermal  properties  of  tube  material,  the  temperature  wave  cannot 
penetrate  more  than  one  millimeter  into  the  tube  wall. 

8.  The  total  heat  loss  to  the  tube  wall  is  five  to  six  percent 
of  the  input  energy  for  typical  cases  and  has  only  a  minor  effect  on 
the  final  ballistic  resalts.  The  same  conclusion  is  valid  for  the  skin 
friction. 

9.  The  piston  start  pressure,  although  difficult  to  determine 
in  practice,  does  not  pose  any  real  problem  due  to  its  insignificant 
effect  on  the  ballistic  solution. 

10.  Improvement  in  ballistic  efficiency  car  be  brought  about,  by 
in-  ►‘easing  the  propellan-  charge  and  the  piston  mass  or  by  reducing 
trie  v.;b  thickness.  In  each  of  these  cases,  there  exists  an  adverse 
eiL<*et  jf  higher  peak  pressure  and  higher  wall  temperature.  Therefore, 


great  care  should  be  taken  to  obtain  the  optimum  design  conditions. 

11.  The  heat  transfer  model  is  verified  by  the  limited  wall 
temperature  data  available  at  this  time.  This  agreement  is  gained  with¬ 
out  the  use  of  experimentally  measured  interior  ballistic  data  such  as 
chamber  pressure  versus  time,  as  used  in  other  analysis  [25]. 
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PART  II  -  PARTICLE  SEEDING  EFFECTS  ON  GUN  BARREL 
HEAT  TRANSFER  AND  EROSION 
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SECTION  VII 


INTRODUCTION 

The  seeding  of  gun  propellants  with  additives,  such  as  Ti02  in  a 
wax  matrix,  appears  to  be  an  effective  method  of  reducing  gun  barrel 
erosion  [26,  27,  28,  29,  30].  The  mechanism  by  which  this  is  accomplished 
has  not,  however,  been  satisfactorily  explained.  This  makes  it  impossible 
to  optimize  the  :ffect  or  properly  predict  it  without  expensive  and  time 
consuming  trial-and-error  experiments. 

Since  gun  barrel  erosion  has  been  shown  to  be  primarily  tube 
temperature  [31]  and  thereby  heat  transfer,  sensitive,  the  quantitative 
heat  transfer  model  developed  in  the  previous  section  will  be  used  to  test 
various  possible  mechanisms  by  which  additives  might  reduce  the  temperature 
of  the  gun  tube. 

The  fact  that  the  reduced  erosion  is  primarily  a  temperature 
phenomena  rather  than  a  chemical  phenomena  is  substantiated  not  only  by 
reference  31  but  also  by  reference  29  which  shows  an  increase  in  CO  due 
to  a  wax  additive.  Reference  31  shows  the  CO  increase  should  increase 
rather  than  reduce  chemical  erosion.  The  fact  that  erosion  is  actually 
reduced  by  the  TiOj  -  wax  additive  causes  one  to  search  for  a  reduced 
heat  transfer  mechanism  created  by  the  additive. 

This  is  carried  out  by  looking  at  (1)  temperatures  reached  under 
repetitive  fire  conditions,  (2)  experimental  results  from  dusty  gas  heat 
transfer  experiments  usually  related  to  gas-cooled  nuclear  reactors,  and 
(3)  insulation  effects  of  particles  imbedded  in  valleys  of  the  barrel 
surface  roughness. 
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SECTION  VIII 


REPETITIVE  FIRE  ANALYSIS 
Introduction 

Since  rapid  fire  aircraft  cannons  are  of  primary  interest  in  this 
effort,  the  wall  temperatures  predicted  in  the  foregoing  single-shot  analy¬ 
sis  are  net  directly  applicable.  Since  this  wall  temperature  is  the 
parameter  of  most  interest,  a  repetitive  fire  analysis  is  carried  out 
here  to  determine  what  peak  temperatures  actually  occur  in  aircraft 
cannons.  This  allows  one  to  determine  if  the  temperatures  are  sufficiently 
high  to  cause  melting  of  a  thin  layer  of  metal  on  the  inside  of  the  barrel. 

Model 

In  carrying  out  a  repetitive  fire  analysis,  one  is  immediately  con¬ 
fronted  with  the  gun  emptying  problem,  i.e.,  calculating  the  gas  flow 
properties  and  gas  to  metal  heat  transfer  rates  after  the  projectile  leaves 
the  tube  but  before  the  next  projectile  is  fired.  The  previous  single¬ 
shot  analysis  is  not  applicable  because  it  assumes  a  solid  projectile 
base  closes  the  tube  at  all  times. 

Rather  than  attempt  to  solve  this  complex  problem,  it  was  decided 
to  carry  out  the  calculations  using  an  upper  and  lower  limit  on  the  heat 
transfer  during  this  emptying  or  soaking  period.  Since  the  total  heat 
transfer  during  this  time  was  expected  to  be  small,  it  was  hoped  that 
these  two  limits  would  put  a  rather  small  bracket  on  the  actual  results. 
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An  upper  limit  is  obviously  taken  if  one  assumes  that  the  gas  tempera¬ 
ture  and  heat  transfer  coefficient  remain  at  their  values  which  existed 
at  the  time  the  projectile  left  the  tube.  After  the  projectile  leaves, 
the  gas  temperature  will  be  reduced,  the  gas  pressure  will  be  reduced,  and 
the  gas  velocity  will  decay  to  zero.  Each  of  these  factors  will  tend  to 
reduce  tne  heat  transfer  to  the  barrel  below  this  assumec  upper  limit. 

A  lower  limit  is  certainly  given  by  assuming  that  the  inside  of  the 
barrel  is  suddenly  insulated  from  the  gases  at  the  instant  the  projectile 
leaver  the  barrel.  This  assumes  that  no  heat  is  transferred  from  the  hot 
gases  to  the  barrel  during  the  soak  period  between  shots. 

The  gas  temperature  and  heat  transfer  coefficient  as  a  function  of 
axial  distance  and  time  were  taken  from  the  single  shot  analysis  and 
repetitively  applied  during  each  and  every  shot  while  the  projectile  was 
moving  down  the  barrel.  The  conduction  analysis  was  started  at  the  time  of 
the  initial  shot  and  successively  calculated  each  At  during  the  firing 
period  and  the  soaking  period.  Since  the  metal  temperature,  T^,  at  any 
t  +  At  depends  on  the  metal  temperature  profile  T  (r)  at  t,  this  conduction 
calculation  must  be  carried  on  continuously  with  the  repetitive  heat  trans¬ 
fer  coefficient,  h, (t),  and  gas  temperature,  T  (t),  being  repeated  for  each 
i  00 

firing  cycle,  as  detailed  in  Chapter  IV,  except  that  the  wall 

temperatures  are  not  reset  to  at  the  beginning  of  each  firing. 

One  should  note  that  this  model  does  not  assume  that  the  total  heat 

transferred  during  each  cycle  is  the  same  from  one  cycle  to  the  next  because 

the  wall  temperature,  T  .  ,is  increasing  from  cycle  to  cycle, thereby  re- 

w ,  l 

ducing  (T  -  T  .  )and  the  heat  transfer. 

*  w,i 
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Conduction  Series  Solution  During  Soaking 


For  short  cyclic  firing  times  (i.e.,  0.010  second  or  less)  the  con¬ 
duction  solution  during  soaking  was  carried  out  by  the  finite  difference 
technique  detailed  in  Chapter  IV.  However,  at  larger  cyclic  firing 
times  (i-s*>  up  to  0.060  second),  the  number  or  time  steps  required 
during  the  soak  period  was  so  large  that  the  computer  computation  time  was 
excessive.  Therefore,  a  closed  form  series  solution  was  developed  and 
used  during  the  time  period  from  the  instant  the  projectile  left  the  barrel 
to  the  time  the  next  shot  was  initiated.  This  series  solution  was  made 
possible  during  soaking  due  to  the  constant  heat  transfer  coefficient,  h^, 
with  respect  to  time.  This  parameter  is  not  constant  with  time  during  the 
firing  itself. 

Taking  the  conduction  equation  and  applying  it  at  some  specific  axial 
location  during  the  nth  soaking  period: 

(145) 

and  separating  variables,  one  finds  that  the  solution  may  be  expressed  as 
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where  J  and  Y  are  zero  order  Bessel  functions  of  the  1st  and  2nd  kind, 
respectively,  and  and  B^  must  be  found  from  the  boundary  and  initial 
conditions  which,  for  the  lower  limit  case  of  zero  heat  transfer  during 
soaking,  are: 

3T 

—  =  0  at  r  =  r,  (147) 

or  i 
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K  ,  t-SL  +  h  (T  -  T  v)  =  o  at  r=r 
w  5r  o  ^  amb  c 


(148) 


Vr>  "  Tamb  =  f(r)  at  t  -  tn 


(149) 


where  t  is  the  time  when  the  nth  projectile  leaves  the  tube  and  the  nth  soak 
period  begins  and  f(r)  is  the  temperature  distribution  with  respect  to  r 
existing  at  that  time. 

Applying  these  boundary  conditions  to  (146),  one  finds: 


(  r) 
ml  m 


1amb  ^ 


Cm  * j7 

m  m  w  1 


J-.  (  r  )  -  h  J  (  _r  ) 

0  0  0  04  o 


VW> 

2 

\n  wYl  ^  mro)  “  ^oYo  ^  mro^ 


e'Vw*  (150) 


where  U,  and  Y,  are  first  order  Bessel  functions  of  the  1st  and  2nd  kind. 


respectively,  ^  is  the  mth  root  of  the  equation: 
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=  0  (151) 


and  C  is  calculated  by  use  of  the  function  in  equation  (149) 


C  =  N  r  R  (X  ,r)  f  (r)  dr 
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Using  an  of 
equation  (151)  are: 


3  cm  and  an  tq  of  6  cm  the  first  five  roots  of 


X±  -  3.1965 
X,  =  6.3275 
Xi  =  9.4522 
Xx  =  12.5844 
X1  =  15.7230 


and  only  five  terms  were  retained  in  the  series.  Accuracy  was  checked 
by  comparing  the  results  from  the  series  solution  containing  five  and  ten 
terms  with  the  finite  difference  solution  results.  This  comparison  showed 
that  five  terms  yielded  sufficient  accuracy. 


Results 


Results  of  the  repetitive  fire  calculations  were  carried  out  for  two 

firing  rates:  0.010  second  between  firings  and  0.060  second  between  firings. 

Since  the  projectile  is  in  the  tube  for  less  than  0.003  second,  the  soaking 

times  are  about:  0.007  and  0.057  second  for  the  two  firing  rates,  respectively. 

The  results  for  three  axial  locations  are  shown  in  Figures  44  through 

46  for  the  0.010  second  firing  time  with  the  heat  transfer  coefficient 

h^  during  soaking  at  its  upper  limit  which  is  equal  to  its  value  at  the 

time  the  projectile  leaves  the  tube,  devoted  by  h^  .  These  are  for 

last 

ten  repetitive  firings. 
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Figure  44,  Repetitive  Fire  Wall  Temperature.  Versus  Time  at 
4xial  Location  of  25  cm  (Upper  Limit) 
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Figure  45.  Repetitive  Fire  Wr LI  Temperature  Versus  Time  at 
Axial  Location  of  12  cm  (Upper  Limit) 
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Figure  46,  Repetitive  Fire  Wall  Temperature  Versus  Time  ?t 
Axial  Location  of  83  cm.  (Upper  Limit) 
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Figures  47  through  49  show  the  same  results  assuming  zero  convective 
heat  transfer  during  soaking,  i.e.,  h^  equal  zero  whicl:  is  the  lower  limit. 
These  sets  of  curves  then  bracket  the  upper  and  lower  limits  for  the  inside 
wall  temperature.  At  an  axial  location  of  25  cm,  which  sees  the  highest 
temperature  (Figure  24),  the  upper  limit  case  after  ten  firings  show  a 
maximum  temperature  of  about  1600°K  compared  with  the  lower  limit  case  of 
about  1500°K.  This  is  a  small  bracket  (8  percent)  and  shows  that  assumption 
concerning  the  heat  transfer  during  the  soaking  period  has  little  effect 
on  the  maximum  wall  temperature  of  interest.  Therefore,  at  a  given  axial 
position,  rubsequent  calculations  assumed  zero  heat  transfer  from  the 
combustion  gas  to  the  tube  walls  from  the  time  projectile  n  leaves  the  tube 
to  the  time  projectile  n+1  passes  that  axial  position. 

Results  using  this  lower  limit  assumption  for  the  lower  cyclic  firing 
time  of  0.060  second/round  are  shown  in  Figures  50  through  55  for  50  repetitive 
firings  at  the  same  three  axial  positions  as  shown  for  the  previous  case  in 
Figures  44  through  49. 

Figure  50  shows  the  wall  temperature  as  a  function  of  time  at  the 
high  temperature  axial  position  of  25  cm.  Figure  51  shows  the  temperature 
profile  inside  the  tube  wall  at  this  same  axial  position  at  four  different 
times.  These  same  type  of  curves  are  repeated  at  the  two  other  axial 
positions  of  12  and  34  cm  in  Figures  52-53  and  Figures  54-55,  respectively. 

With  the  longer  soaking  periods  between  firings,  the  maximum  temperature 
•if tor  50  cycles  is  about  1400°K.  It  is  seen  from  the  wall  temperature 
ver  ,us  time  graphs  that  the  peak  temperature  during  ench  cycle  rises  very 
slowly  niter  the  initial  10  rounds.  Also,  it  is  seen  from  the  radial 
temperature  distribution  inside  the  tube  wall  that  the  outside  tube 
temperature  has  not  risen  significantly  after  the  50  rounds. 
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Figure  47.  Repetitive  Fire  Wall  Temperature  Versus  Time  at 
Axial  Location  of  25  cm  (Lower  Limit) 
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Figure  49.  Repetitive  Fire  Wall  Temperature  Versus  Time  at 
Axial  Location  of  83  cm  (Lower  Limit) 
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Figure  50.  Repetitive  Fire  Wall  Temperature  Versus  Time  at  Axial  Position 
of  25  cm  (0.060  second  firing) 
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Figure  51.  Variation  of  Repetitive  Fire  Wall  Temperature  with 
Depth  at  Axial  Position  of  25  cm  at  Various  Tines 
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of  84  cm  (O.OfaO  si 


[The  reverse 


I 

2  ,n:i  j 

SF.C  J 

Figure  54.  Repetitive  Fire  Wall  Temperature  Versus  Time  at  Axial  Position 
of  84  cm  (0.060  second  firing) 
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Figure  55.  Variation  of  Repetitive  Fire  Wall  Temperature  with  Depth 
at  Axial  Position  of  84  cm  at  Various  Times 


Of  course,  with  soaking  time  between  series  of 
outside  wall  temperature  would  rise  and  result 
t^mnpratures  as  well. 
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SECTION  IX 


PARTICLE  SEEDING  EFFECT  ON  TURBULENT  BOUNDARY 
LAYER  HEAT  TRANSFER 


Introduction 

The  section  will  be  devoted  to  a  discussion  of  the  physical  phenomena  where¬ 
by  certain  additives  (e.g.,  small  titanium  dioxide  particles  in  a  wax  matrix 
dramatically  reduce  erosion  in  gun  barrels.  The  limited  open  literature 
documents  the  effectiveness  of  several  additives;  however,  no  consistent 
physical  explanation  is  offered,  and  no  satisfactory  mathematical  model 
has  been  proposed. 

Since  hot  gas  erosion  rates  appear  to  be  strongly  dependent  on 
temperature  level,  an  effort  has  been  made  to  relate  erosion  reduction  to 
a  reduction  of  the  barrel  inner  surface  temperatures.  Some  mechanism 
must  be  presented  then  which  attenuates  the  heat  transfer  from  com¬ 
bustion  products  to  the  gun  barrel.  Three  distinct  possibilities  have 
been  considered: 

1.  The  additive  significantly  alters  combustion  gas  properties 

so  as  to  reduce  combustion  gas  temperatures  and/or  heat  transfer 
coefficients, 

2.  The  additive  particles  reduce  the  turbulence  level  near  the 
wall,  thereby  reducing  the  heat  transfer  coefficient. 

3.  The  particles  and/or  wax  matrix  deposit  on  the  wall,  forming 
a  protective  coating  which  insulates  the  metal  from  the  hot 
gas. 
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Effec<s  of  Additive  on  Combustion  Products 


The  work  by  Lenchitz  et  al  [29]  indicates  that  TiC^  is  not  in¬ 
volved  in  the  combustion  process,  as  the  heat  of  reaction  remains  essentially 
the  same  when  relatively  large  amounts  of  the  micronized  additive  are  in¬ 
cluded  with  the  propellant.  On  the  other  hand,  the  addition  of  about  10 
percent  by  weight  of  wax  decreased  the  heat  of  reaction  by  about  10  percent. 
This  is  apparently  due  to  a  combination  of  reaction  of  wax  with  the  explosion 
products  and  prevention  of  complete  combustion  of  the  propellant.  Less 
reduction  in  heat  of  reaction  is  observed  when  TiO^  is  added  to  the  wax, 
suggesting  that  the  TiO^  inhibits  the  interaction  between  the  wax  and  com¬ 
bustion  products  of  the  propellant.  There  seems  to  be  a  trade-off  between 
heat  of  reaction  and  quantity  of  gas  such  that  the  pressure  is  unaffected 
by  the  additives;  hence,  the  interior  ballistic  properties  are  observed  to 
be  nearly  the  same. 

From  the  heat  transfer  standpoint,  the  effects  of  interest  are  a 
probable  lowering  of  the  adiabatic  flame  temperature  and  a  change  in  com¬ 
position  of  the  combustion  gases.  In  addition,  the  gas  temperature  could 
be  reduced  due  to  sensible  heating  of  the  TiC^  particles  and  wax  [34]. 

The  magnitude  of  this  cooling  due  to  the  inert  particles  having  to 
be  heated  to  the  flame  temperature  can  be  approximated  by  an  energy 
balance  on  the  gas-particle  mixture.  Since  the  heat  transfer  to  the 
particles  is  equal  to  the  heat  transfer  from  the  gas: 
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Qgas  =  cp  M  AT 

gas  gas  gas 


in  Equation  (155) 
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Since  the  combustion  gas  temperatures  are  around  2000  K  and  the 


particles  start  out  at  ambient  temperature,  AT  is  about  1700°K. 

Using  cD  of  0.412  k  cal/kg  -°K,  °p  of  0.23  k  cal/kg  -°K,  and 
^gas  part 

M  /M  of  0.05: 
part  gas 


AT  =  -28°K 

gas 


This  is  obviously  much  too  small  of  a  gas  temperature  drop  to  result  in 
a  significant  tube  wall  temperature  reduction. 

Effects  of  Particle  Seeding  on  Turbulent  Heat  Transfer 

There  has  been  considerable  interest  in  recent  years  in  the  fluid 
mechanical  and  heat  transfer  characteristics  of  suspensions,  both  liquid- 
solid  and  gas-solid.  Gas-solid  flows,  for  example,  have  been  extensively 
studied  during  the  past  15  years  in  connection  with  pneumatic  conveying 
of  solids,  coolants  for  gas-cooled  nuclear  reactors,  etc. 

Liquid-Solid  Suspensions 

The  technology  of  drag  reduction  phenomena  dates  back  to  1945  when 
it  was  discovered  that  under  turbulent  flow  conditions,  gasoline-aluminum 
soap  gels  exhibited  much  less  hydraulic  resistance  than  untreated  gasoline  [35]. 

Since  then,  hierarchies  of  drag  reduction  effects  have  been  found 
in  liquids  containing  low  molecular  weight  micelle-type  additives,  high 
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molecular  weight  polymers,  and  particulates.  The  decrease  in  drag  or 
pressure  drop  is  usually  attributed  to  the  non -Newtonian  nature  of  the 
suspension.  The  solids  absorb  energy  from  the  stream,  thereby  de¬ 
creasing  the  turbulence  level  possibly  to  the  point  of  laminarizing 
the  flow. 

In  spite  of  the  more  than  400  references  in  the  open  literature,  it 
seems  that  practical  exploitation  of  this  technology  area  has  not  yet 
been  realized.  Furthermore,  there  is  relatively  little  information 
available  regarding  the  heat  transfer  characteristics  of  these  suspensions. 
The  natural  temptation,  however,  is  to  assume  that  heat  transfer  will  be 
reduced  in  accordance  with  the  analogy  between  heat  and  momentum  transfer. 
Alteration  of  Transport  Properties  of  Particle  Seeding . 

It  is  appropriate  at  first  to  turn  to  experience  on  the  dynamics  of 
dense  clouds  of  particles  or  dusty  gases.  The  classical  Einstein  relation 
for  an  incompressible  fluid  containing  solid  spheres  gives  the  viscosity 
as  [36] 


*  u  (1  +  2. Sot) 
m  g 


(159) 


where  um  and  u  are  the  viscosity  of  the  mixture  and  pure  gas  respectively 
and  a  is  the  volume  fraction  of  solids.  In  order  to  estimate  the  increase 
iu  viscosity  due  to  solids  addition,  assume  that  TiOj  is  added  in  the 
large  proportion  of  20  percent  by  weight  of  propellant.  Presume  that  the 
propellant  undergoes  complete  combustion.  The  volume  fraction  of  TiO^, 
for  no  relative  velocity  between  phases,  is  then 
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where  M  is  the  mass  of  the  solids,  (  )  ,  or  gas.  (  )  .  In  other  words, 

s  g 

the  loading  ratio  is  too  low  to  cause  any  significant  increase  in 
viscosity.  Even  if  the  local  concentration  at  the  wall  were  much  higher, 
it  is  probable  that  the  density  would  be  sufficiently  elevated  so  that 
the  Reynolds  number  would  be  increased  rather  than  decreased.  The  usual 
description  of  dusty  gas  flows  is  that  the  particle  addition  results  in 
a  flow  of  higher  effective  Reynolds  number  [37,  38].  The  heat  transfer 
coefficient  would  then  be  expected  to  increase. 

This  is  consistent  with  theoretical  and  experimental  work  which 
suggests  that  particle-particle  interactions,  which  could  increase 
viscous  drag,  are  negligible  for  particle  volume  fractions  less  than 
0.05  [39]. 

Studies  of  Gas-Solids  Suspensions  (>10  yi  Particles) 

During  the  past  15  years  several  studies  have  been  reported 
of  heat  transfer  to  (or  from)  gas-solid  suspensions,  with  particle  sizes 
generally  in  the  30  to  200p  range  and  with  relatively  large  loading  ratios. 
The  results  have  been  applied  to  numerous  industrial  applications, in¬ 
cluding  coal-fired  boilers,  heat  recovery  equipment,  solid  propellant 
rockets,  and  gas-cooled  nuclear  reactors.  Numerous  investigators  have 
reported  a  decrease  in  local  heat  transfer  coefficients  at  low  solids 
loading  ratios.  This  observation  is  of  direct  interest  to  the  present 
study. 

Farber  and  Morley  t^O]  made  average  measurements  of  heat  transfer 
from  a  steam-heated  tube  to  air-alumina  catalyst  mixtures.  No  con¬ 
sistent  trend  of  decreasing  coefficient  can  be  discerned  in  the  data, 
which  exhibit  consJderable  scatter  at  low  loading  ratios.  It  was  noted 
that  bulk  temperature  rises  were  lower  (and  AT^'s  higher)  with  the 
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suspensions  than  for  pure  air;  however,  this  effect  was  more  than  com¬ 
pensated  for  by  the  increased  heat  transfer  coefficient,  i.e., 
the  heat  transfer  rate  was  always  greater  with  solids  addition. 

One  of  the  few  analytical  attempts  was  reported  by  Tien  [4] 
who  considered  the  case  of  mass  loading  ratios  lower  than  about  1.0. 

A  linear  increase  in  local  coefficient  with  loading  ratio  was  predicted. 

In  addition,  the  solid  particles  were  expected  to  have  negligible 
effect  on  heat  transfer  far  downstream  from  the  thermal  entrance. 

Average  predictions  compared  reasonably  well  with  the  data  of  Farber 
and  Morley. 

Farber  and  Depew  [42]  took  average  data  which  demonstrated  a 
distinct  decrease  in  coefficient  at  low  loading  ratios;  however,  no 
particular  notice  was  taken  of  this  behavior.  The  overall  technique 
thus  precludes  a  clear  look  at  what  is  happening  to  the  local  heat 
transfer  coefficient. 

Depew  and  Farber  [43]  subsequently  reported  local  data  which  clearly 
demonstrated  decreases  in  heat  transfer  coefficients  at  low  loading 
ratios,  as  shown  in  Figure  57.  A  theoretical  analysis  qualitatively  pre¬ 
dicted  this  decrease;  however,  the  predicted  decrease  became  negligible 
as  the  particle  size  was  reduced  to  30  microns,  a  trend  which  was  con¬ 
trary  to  the  data.  Tien  and  Quan  [44]  also  took  local  data  (Figure  57) 
which  demonstrated  the  decrease  in  coefficient  at  low  loading  ratios. 

It  was  speculated  that  the  particles  altered  the  flow  field. 

Briller  and  Peskin  [45]  reported  average  measurements  of  heat 
transfer  to  gas-solids  mixtures  and  suggested  that  there  was  no  observable 
increase  in  heat  transfer.  This  conclusion  is  perhaps  justified  consider¬ 
ing  the  inordinately  large  scatter  in  data;  however,  a  closer  examination 
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Transfer  Coefficient  Variation  with  Figure  57.  Keat  Transfer  Coefficient 

cle  Loading  (Briller  and  Peskin)  Variation  with  Particles 


of  the  cooling  data  shown  in  Figure  56  suggests  a  distinct  decrease  in  the 
coefficient  at  low  loading  ratio  since  the  heating  data  exhibited  no 
particular  dependence  on  M  /M.  This  could  suggest  that  decreases  are  more 
pronounced  with  cooling. 

Soo  [46]  has  developed  an  interesting  analogy  explanation  for  the 
effect  of  solids  on  heat  transfer  at  low  loading  ratios.  The  decrease 
in  friction  factor [44]  in  the  presence  of  solid  particles  (Figure  5S) 
is  attributed  to  reduced  mixing  length  due  to  dissipation  by  the  solid 
particles.  Upon  introduction  of  the  reduced  shear  stress  into  the 
Reynolds  analogy,  a  Nusselt  number  -  loading  ratio  expression  is 
determined  which  is  similar  to  the  experimental  behavior,  as  shown  in 
Figure  59.  While  this  approach  is  interesting,  it  serves  only  to  express 
one  dependent  variable  in  terms  of  another;  hence,  little  insight  is 
given  into  the  mechanism  of  transport  process  reduction. 

Rossetti  and  Pfeffer  [47]  have  recently  reported  friction  data 
which  are  of  interest  in  the  application  of  the  analogy  concept.  As 
illustrated  in  Figure  60, when  glass  beads  of  10  to  60u  were  added  to  air  at 
Reynolds  numbers  of  10,000  to  25,000,  friction  factor  reduction  up  to 
75  percent  was  observed  in  horizontal  test  sections.*  In  order  to 
clarify  this  effect,  turbulence  intensities  were  measured  at  the  pipe 
centerline.  The  intensity  of  turbulence  was  always  greater  with  soiids 
addition,  an  effect  which  is  in  direct  conflict  with  the  drag  reduction. 

It  was  postulated  that  the  interaction  of  particles  wich  the  turbulence 
eddies  actually  occurs  very  close  to  the  wall  and  that  at  the  wall  the 
turbulence  decreases. 


*  It  is  noted  that  the  reduction  is  real  since  the  suspension  friction 
factor  is  evaluated  on  the  basis  of  the  gas-phase  density.  Use  of  a 
suspension  density  will  always  yield  a  reduced  friction  factor. 


Brandon  [48]  has  also  pointed  out  the  need  for  wall  turbulence  measurements 
as  a  necessary  step  toward  understanding  friction  and  heat  transfer  be¬ 
havior  of  gas-solids  suspensions. 

Survey  Conclusions. 

A  detailed  examination  of  the  data  suggests  that  the  following  con¬ 
clusions  can  be  drawn  from  the  experimental  studies  cited  above. 

1.  The  decrease  in  heat  transfer  coefficient  is  generally  larger  with 
smaller  particles  [44]. 

2.  The  percentage  decrease  in  h  is  smaller  with  higher  Ifeynolds 
numbers  [44]. 

3.  The  loading  ratio  at  which  the  minimum  h  occurs  seems  to  decrease 
as  particle  size  is  reduced;  however,  increases  in  Ifeynolds  number  shift 
the  minimum  _o  higher  values  of  loading  ratio  [44]. 

4.  The  decrease  in  h  is  greater  for  lead  particles  than  for  glass 
particles  [44]. 

5.  Thermal  entrance  lengths  are  larger  for  glass  particles  than  for 
lead  particles.,  and  are  considerably  larger  for  smaller  particles.  In  all 
cases,  however,  the  entrance  length  was  increased  over  the  pure  gas 
value  [44] . 

6.  The  decrease  in  heat  transfer  coefficient  at  low  loading  ratios 
appears  to  be  more  pronounced  with  cooling  than  with  heating  [45]. 

7.  Several  available  analytical  and  semi-analytical  studies  of 

heat  transfer  and  friction  do  not  satisfactorily  account  for  the  observed 
reduction  in  heat  transfer  and  friction  at  low  loading  ratios  [41,46,49]. 

Observation  1  is  encouraging  as  it  suggests  that  still  smaller 
particles  (<30p)  would  exhibit  greater  decreases  in  h.  However, 
observations  2  through  5  are  either  nebulous  or  suggest  a  trend  In  the  opposite 


161 


direction  to  what  is  desired  to  explain  reduced  heat  transfer  in  gun 
tubes  with  seeding.  Observation  6  is  favorable;  however,  the  data  are 
too  limited  to  draw  firm  conclusions.  Finally,  it  is  noted  from 
observation  7  that  the  theory  is  of  little  help  in  suggesting  what  will 
happen  with  micron-size  additives. 

Studies  of  Gas-Solids  Suspensions  (~ly  Particles) 

Particles  in  the  size  range  of  more  direct  interest  have  been 
employed  in  studies  relating  to  augmentation  of  heat  transfer  in  gas- 
cooled  nuclear  reactors.  The  general  idea  has  been  to  add  graphite  to 
the  helium  coolant  in  such  proportions  as  to  appreciably  increase  the 
coolant  specific  heat  and  the  heat  transfer  coefficient.  Graphite  is 
the  normal  additive  since  it  will  pose  few  problems  for  the  circulating 
equipment;  regardless  of  the  size  of  addition,  the  particles  shortly 
pulverize  to  an  average  size  less  than  lp. 

Schuderberg  et  al  [50]  summarized  an  extensive  program  at  Babcock  and 
Wilcox  to  study  the  characteristics  of  dispersions  of  micron-size 
graphite  particles.  Friction  factors  at  higher  Reynolds  numbers 
dropped  below  the  nominal  value.  It  was  suggested  that  flow  turbulence 
was  suppressed  by  the  suspended  particles.  By  analogy,  heat  transfer 
should  be  affected  to  a  comparable  extent.  Unfortunately,  the  heat 
transfer  tests  were  done  only  with  loading  ratios  considerably  in  excess 
of  unity. 

Woodcock  and  Worley  [51]  summarized  extensive  Ehglish  studies  of 
boiler  performance  m  a  gas-cooled  nuclear  reactor  system,  where  the 
CC>2  coolant  was  loaded  with  graphite.  Only  loading  ratios  greater 
than  4  were  considered,  and  the  average  heat  transfer  coefficient  for 
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the  tube  bank  always  increased.  One  interesting  aspect  of  their  work 
was  that  the  effective  particle  sizes  appeared  to  be  much  larger  than 
those  of  the  basic  particle.  It  is  doubtful,  however,  that  the  load¬ 
ing  ratios  encountered  with  barrel  wear-reducing  additives  are  high 
enough  for  the  mixture  to  exhibit  this  effect. 

One  of  the  more  serious  difficulties  encountered  with  gas- 
graphite  suspensions  is  the  deposition  of  the  solids  on  cold  surfaces  [52]. 
This  is  apparently  due  to  a  combination  of  Brownion  motion  and 
electrostatic  effects.  For  the  present  application,  this  observation 
suggests  that  particle  concentration  near  the  wall  will  be  relatively 
high.  It  would  seem,  however,  that  the  high  velocities  and  short  firing 
times  would  tend  to  limit  the  buildup  of  as  heavy  particle  layer  as  is 
experienced  in  gas-cooled  reactor  system  tests.  The  effects  of  any 
particle  buildup  will  be  examined  in  the  next  section  of  this  report. 

The  IITRI  Analysis 

In  1966  the  Illinois  Institute  of  Technology  Research  Institute 
was  awarded  a  contract  by  Picatinny  Arsenal  to  explain  the  mechanism  by 
which  an  additive  reduces  gun  barrel  erosion.  While  the  primary 
emphasis  in  this  investigation  was  devoted  to  tests  with  a  vented 
chamber,  an  analytical  model  was  also  developed  as  described  in  the 
summary  report  by  Remaly  and  Stanley  [52], 

The  basic  premise  of  this  model  is  that  reduced  erosion  is  related 
to  a  reduction  in  heat  transfer  from  the  gas  to  the  metal.  Following 
the  qualitative  explanation  offered  for  reduction  of  transport  rate  with 
solid-liquid  and  gas-solid  suspensions,  turbulent  damping  was  hypothesized. 
The  gas  particles  were  assumed  to  be  in  simple  harmonic  motion;  the 
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particles  then  interact  with  t gas  through  a  dasnpot  that  represent*  *. he 
drag.  The  equations  of  motion  vere  then  solved  numerically  for  various 
spectra  of  turbulence..  The  output  parameter  of  interest  i«=  the  loss  of 
turbulent  kinetic  energy  due  to  injection  of  the  particles.  The  rate  of 
energy  loss  was  then  computed  for  various  particle  diameters  (0.1  to  10  microns), 
densities  (2  to  8  g/cc),  and  concentration  (2  to  9  percent  by  weight). 
Computations  were  limited  by  computer  running  time  to  the  first  three  cycles 
of  gas  oscillation  (about  100  ysec).  Typical  increases  of  5  percent  were 
found  in  the  boundary  layer  thickness. 

The  computations  indicated  that  turbulent  damping  increases  with  particle 
diameter,  density,  and  concentration.  Typical  results  are  showr  in  Figure 
61.  The  experimental  program  was  then  designed  to  provide  at  least  a  quali¬ 
tative  verification  of  the  analytical  model. 

Controlled  Firing  Tests 

IITRI  Tests 

Remaly  and  Stanley  [52]  tested  the  effects  of  various  propellant 
additives  in  a  37  mm  vented  chamber.  Additive  effectiveness  was  by 
measuring  wear  of  an  orifice-type  plug. 

A  variety  of  additives  was  tested  for  multiple  firings  since  weight 
loss  data  from  single  shots  were  inconclusive.  No  wax  was  utilized  in 
tests  presented  in  this  report,  and  the  additives  were  introduced  in  such 
a  fashion  that  good  mixing  with  the  explosion  products  was  obtained. 

All  of  the  additives  (TiC^,  talc,  WO^,  and  microballoons)  were 
found  to  be  effective  in  reducing  erosion  of  the  orifice  insert.  The 
data  shown  in  Figure  62  are  typical  of  the  results  obtained.  Low  density 
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Figure  62.  Cumulative  Weight  Loss  Versus  Shot  Number  (Remaly  and 
Stanley) 
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materials  ware  found  to  be  the  most  effective  in  contrast  to  the  results 
of  the  model  described  in  the  previous  section.  Furthermore,  small  particles 
gave  greater  erosion  reduction,  which  was  again  an  effect  opposite  to  that 
predicted  (Figure  61).  A  re-examination  of  the  model  was  thus  suggested. 

This  leads  one  to  suspect  that  the  basic  premise  of  turbulence  alteration  may 
not  be  correct. 

One  of  the  incidental  observations  of  these  tests  was  concerned  with  a 
residue  in  the  orifice  insert.  This  propellant  residue  was  apparently 
substantial  since  soaking  in  a  chlorinated  hydrocarbon  solvent  was  necessary 
to  remove  it  before  the  weight  was  recorded.  While  no  indications  of  the 
thickness,  composition,  or  influence  of  the  residue  on  erosion  were  given, 
it  seems  reasonable  to  expect  that  the  residue  could  play  a  role  in  reducing 
the  erosion. 

Picatinny  Arsenal  Tests 

The  Swedish  Addillve,  titanium  dioxide  or  lungs tun  trioxid*  in  a  wax 
matrix,  was  adopted  for  U.S.  military  application  in  1961.  This  additive 
was  highly  effective  in  reducing  (up  to  90  percent)  bore  erosion  in  a  variety 
of  gun  barrels.  The  empirical  results  of  previous  tests  reported  by  Swedish 
consultants  indicated  that  many  metal  compounds  have  wear-reducing  proper¬ 
ties.  The  percentage  composition  of  the  metal  compound  and  wax  and  the 
location  of  the  coated  cloth  liner  were  found  to  be  important. 

Several  reasons  were  suggested  for  the  effectiveness  of  the  additive. 

The  wax  component  vaporizes  and  forms  a  coul  insulating  gas  layer  along  the 
wall,  the  metal  oxide  particles  in  the  flew  have  a  high  reflectivity  and 
reduce  radiant  heat  transfer  from  the  propellant  gas  to  the  tube,  and 
the  metal  oxide  particles  form  a  protective  coating  on  the  tube  wall 
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so  as  to  reflect  heat,  away  from  the  surface. 

Extensive  tests  were  undertaken  at  Picatinny  Arsenal  to  evaluate 
various  additives  and  suggest  a  mechanism  for  reduced  erosion  which  is 
consistent  with  the  test  firing  data.  The  first  major  report  issued 
is  by  Wolff  [26]  who  described  extensive  tests  with  90  and  i05mm  tank 
ammunition.  The  results  of  a  preliminary  series  of  tests  with  W0^  wax 
additive  and  a  full  series  of  test  with  TiO£  wax  additive  are  summarized 

in  Figure  63.  Control  tests  with  no  additive  and  with  wax  additive  only 
are  also  indicated.  The  general  conclusion  is  that  any  additive  is 
better  than  ordinary  untreated  ammunition.  The  Swedish-type  additives 
are  definitely  superior  to  the  other  laminar  coolant.  The  TiO additive 
seems  to  be  better  than  the  W0^  additive;  however,  a  larger  number  of 
firings  for  the  W03  would  seem  to  be  necessary  to  confirm  this  finding.  Most 
surprising  perhaps  is  the  observation  that  the  plain  wax  is  also 
reasonably  effective. 

This  report  is  strictly  empirical  and  does  not  pursue  the  mechanism 
of  barrel  wear  reduction.  The  mechanism  advanced  by  the  Swedish  con¬ 
sultants  concerning  reflection  of  radiation  is  not  a  reasonable  explana¬ 
tion  because  the  radiative  heat  transfer  from  the  gas  to  the  tube  is 
negligible  compared  to  the  convective  heat,  transfer.  The  maximum 
radiative  heat  transfer  is  given  by: 


4  4 

q  .  =  o  (T  -  T  ) 

^rad  gas  w 


(161) 


which  for  T  of  500°K  and  T  of  2000°K  yields  a  q  ,  of  20.6  kcai/rn  -sec 
w  gas  J  rad 

compared  with  the  average  convective  heat  flux  calculated  in  Part  I 

2 

ol  50,000  kcal/m  -sec. 
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Comparison  of  Tube  Wear  Versus  Number  o 
Gun  Tubes  (Wolff) 


On  the  other  hand,  the  existence  of  a  residue  on  the  barrel  wall  is  signi¬ 
ficant.  The  coating  was  readily  removed  before  measurements  were  made.  For 
the  main  series  of  tests  a  chemical  analysis  indicated  that  the  coating  was 
mainly  TiO^. 

Subsequent  tests  at  the  Picatinny  Propellants  Laboratory  were  reported 
by  Lenchitz  et  al  [36].  Heir  conclusions  relating  to  the  effect  of  addi¬ 
tives  on  the  combustion  products  is  given  early  in  this  section  of  the 
report.  While  this  effect  tends  in  the  right  direction,  the  coating  effect 
of  the  additive  appears  to  be  of  greater  significance.  Lenchitz  and  his 
co-workers  suggested  that  the  wax  is  the  prime  factor  in  building  up  an 
effective  insulating  coating.  The  role  of  TiC^  was  considered  to  be  that 
of  strengthening  and  permitting  proper  dispersion  of  the  wax.  The  test 
results  with  a  3/8-inch  vented  chamber  (see  Table  V)  abstracted  from  their 
work  was  used  to  justify  this  conclusion: 

TABLE  V.  EXPERIMENTAL  EFFECTS  OF  ADDITIVES  ON  EROSIVITY 

Sample  Proportions  Erosivity  mg  at 

16,000  psi  (5-shot  averaged 


M-2  Propellant 

35g 

24.1 

Propellant  +  TiO^ 

35g  v  2g 

19.0 

Propellant  +  wax 

35g  +  2-4g 

2.8 

Propellant  +  wax/TiO 

35g  +  4g 

0.9 

(552/45%)  2 

The  wax  is  clearly  indicated  as 

the  major  erosion  reducer. 

In  addition 

small  chemical  effects,  the  thermal  insulation  of  the  carbonized  wax  was 
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considered  to  be  important. 

Taking  an  overall  look  at  the  Picatinny  Arsenal  test,  it  is  reason¬ 
able  to  conclude  that  barrel  wear  reducing  additives  produce  a  deposit 
on  the  tube  wall  which  thermally  insulates  the  metal  from  the  hot  com¬ 
bustion  gases.  This  is  consistent  with  the  observation  that  either 
wax  alone,  TiO^  (WO^)  alone,  or  TiO^/wax  (WO^/wax)  reduces  erosion.  In 
all  cases  a  deposit  was  noticed  on  the  barrel.  The  IITRI  tests  mentioned 
previously  are  in  substantial  agreement  for  the  effects  of  TiO^  WO3 
as  well  as  for  talc  and  microballoons.  It  now  remains  to  demonstrate 
that  thin  layers  of  any  of  these  materials  can  produce  sufficient  re¬ 
duction  in  metal  temperature  to  inhibit  erosion. 
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SECTION  IX 


MECHANISMS  FOR  REDUCED  EROSION  WITH  TiO 
PARTICLES  1 

Introduction 

In  light  of  the  repetitive  firing  analysis  and  the  survey  of  the  dusty 
gas  heat  transfer  literature,  the  possible  mechanisms  by  which  reduced 
erosion  may  result  from  TiO^  seeding  of  the  propellant  was  studied  for 
their  plausibility.  The  fact  that  the  repetitive  firing  analysis  showed 
that  temperatures  on  the  inside  surface  of  the  barrel  may  reach  the  metal 
melting  temperature  plus  the  findings  in  the  literature  led  to  concentrating 
efforts  on  two  wall  temperature  reducing  mechanisms.  The  first  to  be  dis¬ 
cussed  is  the  reduced  heat  transfer  coefficient  resulting  from  the  particles 
reducing  the  turbulence  level  and  the  second  is  the  insulation  of  the  tube 
wall  by  particles  forming  a  packed  bed  in  the  crevices  of  the  tube  surface 
roughness. 

Reduced  Heat  Transfer  Coefficient 

As  shown  in  the  previous  section,  the  literature  in  the  dusty  gas 
area  does  not  include  any  experimental  data  for  the  Reynolds  numbers  range 

g 

(~10  )  and  particle  size  (<1  micron)  applicable  to  the  present  gun  problem. 
Also,  all  experimental  work  has  been  carried  out  at  steady  conditions, 
while  the  gun  flow  is  highlv  unsteady. 

However,  it  can  be  seen  from  Figures  56  through  60  that  a  reduc¬ 
tion  in  the  heat  transfer  coefficient  of  around  20  percent  has  been 
demonstrated  in  several  cases.  There  are  also  indications  that  the  higher 
Reynolds  numbers  and  smaller  particle  sizes  may  lead  to  an  even  higher 
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reductions  in  the  heat  transfer  coefficient,  h,  particularly  at  light  load¬ 
ings.  It  may  thereby  be  argued  that  a  20  percent  reduction  in  h  idght  occur 
with  the  addition  of  to  gun  propellants.  There  is  evidence  that  it 

might  possibly  occur  but  no  experimental  proof  that  it  does.  The  plausibility 
of  a  20  percent  reduction  in  h  causing  a  sufficient  reduction  in  the  wall 
temperature  to  produce  a  reduced  erosion  rate  was  studied. 

The  model  and  program  for  calculating  heat  transfer  and  temperatures 
in  the  tube  presented  in  Section  II  was  used  for  this  calculation.  The  h 
which  the  program  calculated  was  reduced  by  20  percent  and  this  reduced  h 
was  used  in  the  heat  conduction  analysis  to  find  the  tube  wall  temperature. 

The  results  are  shown  in  Figure  64  where  the  wall  surface  temperature 
is  displaced  as  a  function  of  time.  The  peak  temperature  at  the  hottest 
location  (x=25  cm)  is  shown  to  be  about  950°K  compared  to  1050°K  from  the 
normal  heat  transfer  coefficient  results  shown  in  Figure  27. 

This  100°K  drop  which  would  result  for  the  20  percent  lower  h  value 
would  appear  to  be  large  enough  to  substantially  reduce  the  erosion. 

This  is  particularly  true  when  we  consider  that  the  temperature-related 
erosion  mechanism  would  be  expected  to  have  somewhat  of  a  threshold  temperature. 
Above  this  threshold  temperature  the  metal  would  possibly  be  liquid  so  as  to 
erode  away  easily,  whereas  below  this  temperature  the  metal  would  be  solid 
with  strength  to  resist  erosion.  This  100°K  reduction  could  bring  the  metal 
surface  below  this  threshold  temperature  and  thereby  substantially  reduce  the 
temperature-related  erosion.  The  remaining  erosion  would  be  primarily  related 
to  other  mechanisms  such  as  chemical  or  friction.  This  could  explain  the 
factor  of  10  reduction  in  erosion  rate  experienced  in  tests  using  Ti02  addi¬ 
tives,  however,  many  assumptions  would  have  to  be  verified  to  inclusively 
show  this. 
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rature  (°K) 


Insulation  Effect 


In  order  to  assess  what  the  effect  an  insulating  layer  of  Ti0o 
particles  forming  on  the  tube  surface  would  have  on  the  wall  temperature, 
the  model  and  computer  program  were  modified.  The  increased  thermal 
resistance  was  included  in  a  new  effective  convective  heat  transfer 
coefficient  h  ^  given  by: 


h 


eff 


h 

1  +  h(AL/K) 

ins 


(162) 


where.  AL  is  the  thickness  of  the  insulating  layer  and  K,  is  the  thermal 
conductivity  of  the  insulation. 

The  order  of  magnitude  of  AL  can  be  found  by  considering  the  surface 
roughness  of  the  tube  wall  and  assuming  the  interstices  are  filled  with 
TiC^particles .  These  interstices  are  on  the  order  of  5  to  10  microns 
(0.0002  to  0.0004  inches)  into  which  the  micron  size  TiO ^  particles  can 
be  packed  forming  a  packed  bed  5  to  ±0  microns,  thick.  The  order  of 

-4 

magnitude  of  tne  thermal  conductivity  for  packed  beds  is  around  4x10 

kcal/sec-tn-°K.  This  low  value  for  packed  beds  is  due  to  the  voids 

existing  between  the  individual  packed  particles. 

The  results  of  the  calculations  for  a  single  shot  firing  are  shown 

in  Figures  65  and  66  assuming  5  and  10  microns  for  AL, respectively .  This 

calculation  should  be  compared  with  Figure  27  in  Section  II  which  is  the 

identical  case  without  insulation  It  can  be  seen  that  the  effect  of 

TiO^  particles  embedded  in  the  surface  roughness  crevices  is  to  lower  the 

peak  metal  surface  temperature  from  about  1050°K  shown  in  Figure  27  to 
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original 


Figure  66.  Variation  of  Wail  Surface  Temperature  at  Certain  Fixed 
Locations  with  Time  (Case  I,  5-Micron  Insulation) 


825°K  for  5  nicrons  of  insluation  or  650°K  for  10  microns  of  insulation. 

This  is  a  very  substantial  decrease  of  225°K  and  400°K  in  the  metal 
temperature  and  should  account  for  the  decreased  erosion  rate. 

Tn  order  to  further  assess  the  plausibility  of  this  hypothesis,  a 
calculation  has  been  carried  out  to  determine  the  mass  of  1102  particles 
necessary  to  coat  the  tube  with  a  5-micron-thick  layer.  Taking  the  barrel 
dimensions  used  throughout  this  work  of  200  cm  in  length  and  3cm  inside 
diameter,  along  with  a  packed  bed  Ti02  density  of  4  gm/cm  ,  yields  a  mass  of 
3.8  gm  of  1i02.  For  the  prt pellant  loading  of  172  gm,  this  is  about  a 
2  percent  loading  fraction  compared  to  typical  loading  fractions  of  5  percent. 
Therefore,  *~here  are  adequate  TL02  particles  in  the  propellant  so  that  after 
only  a  few  firings  the  interstices  of  the  surface  roughness  can  be  filled 
with  an  insulating  layer  ot  Ti02  particles. 
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SECTION  XI 


CONCLUSIONS 

The  following  conclusions  can  be  drawn  from  the  findings 
presented  in  Part  II  of  this  report. 

1.  Repeated  firings  every  60  milliseconds  produced  increasing 
gun  tube  wall  temperatures  to  1500°K  after  fifty  firings.  At  this 
time  (3  seconds)  the  thermal  wave  does  not  reach  the  outside  of 

the  1.5-cm-thick  wall.  Repeated  firing  bursts  could  produce 
temperatures  high  enough  to  cause  melting  of  a  very  thin  layer 
of  metal  inside  the  tube.  Even  after  30  firings,  the  temperature 
drops  very  rapidly  with  increasing  distance  from  the  inside  wall 
surface,  i.e.,with  a  surface  temperature  of  1500°K,  at  a  depth  of 
0,1mm  below  the  surface  the  temperature  is  typically  below  800°K. 

2.  A  review  of  the  existing  literature  concerning  experimental 
measurements  on  dusty  gases  reveal  that  a  reduction  in  wall 
friction  and  convective  heat  transfer  is  possible  in  many  cases 

of  lightly  loaded  dusty  gases.  Reductions  up  to  30  percent  have 
been  measured  in  some  cases.  However,  no  experimental  measures  of 
any  kind  are  available  which  are  applicable  to  the  gun  tube  gas 
flow  problem. 

3.  Available  experimental  data  on  gun  propellant  additives 
containing  micron  size  particles  all  produce  evidence  of  a  deposit 
on  the  inside  surface  of  the  gun  tube. 

4.  If  a  reduction  of  the  heat  transfer  coefficient  of  20 
percent  were  caused  by  the  presence  of  micron-size  particles  In  the 
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combustion  gases,  the  inside  wall  temperature  would  be  reduced  by  about 
100°K.  This  20  percent  reduction  is  considered  possible,  but  experimental 
evidence  of  this  reduction  applicable  to  this  specific  problem  is  not 
available. 

5.  Insulation  of  the  tube  wall  from  the  hot  combustion  gases  can 
result  from  the  micron-size  particles  becoming  packed  into  the  metal 
surface  roughness  crevices.  These  crevices  have  dimensions  on  the  order 
of  five  to  ten  microns;  therefore,  the  packed  bed  of  particles  in  these 
crevices  form  an  insulating  layer  of  five  to  ten  microns  thickness  between 
the  wall  and  the  gases.  Taking  typical  values,  this  insulation  layer  can 
reduce  the  tube  metal  temperature  by  up  to  300°K.  This  hypothesis  is 
supported  in  many  ways,  and  it  is  felt  that  it  represents  the  mechanism 
by  which  gun  tube  erosion  is  reduced  when  TiQ-j  particles  are  added  to 
propellants.  Theiefore,  in  order  to  minimize  erosion  in  gun  tubes,  this 
deposition  of  any  insulating  particles  into  the  surface  roughness  crevices 
el.ould  be  optimized. 
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APPENDIX  I 


*9 


DERIVATION  OF  CONSERVATION  EQUATIONS 
A  control  volume  approach  has  been  taken  for  the  derivation  of 
the  conservation  equations.  In  Figure  37  a  control  volume  having  a 
cross-sectional  area  A^  and  length  Ax  has  been  shown.  The  volume 
fraction  of  solids  per  unit  length  is  v  ,  or  in  the  other  words,  is 
the  fraction  of  the  total  cross-sectional  area  A^  occupied  by  the  solids. 
Therefore,  (1-v  )  is  the  fractional  area  occupied  by  the  gases  at  any 
position  and  time. 

Due  to  the  assumptions  regarding  the  burning  rate  of  the  solids 

(same  for  all  the  particles  at  a  particular  instant)  and  the  constant 

total  burning  surface  Sb  ,  it  is  easy  to  estimate  the  burning  surface 

t 

available  in  the  chosen  control  volume, 


(1-1) 


Therefore,  the  rate  of  gas  produced  (by  mass)  from  the  solids  (or,  rate 
of  decrease  of  solids  by  mass)  within  the  control  volume  is  given  by, 


p  S.  r 
s  L  b 


»  Pc  SK 
s  b 


v  Ax  r, 
s  b 


c  (  P 

v  dx 


(1-2) 
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Case  I  (Solid  particles  moving  at  gas  velocity) 
Continuity  of  solids. 


Rate  of  increase  of 
solid  mass  in  c.v. 


Rate  of  solid 
flowing  in 


Rate  of  solid 
flowing  out 


Rate  of  gas 
produced  in  c.v. 


t —  (p  A  v  Ax)  =  (p  U  A  v  ) 

dt  s  p  s  s  s  p  s 


(p  U  A  v  ) 
s  s  p  s 
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with  U  =  U  =  U, 
s  g 


3v  3v  b.  v  r 

3t  3x  s  3x  rtp  '  ~jx 
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Continuity  of  Gases, 


Rate  of  increase  of  _  Rate  of  gas 
mass  of  gas  in  c.v.  flowing  in 


Rate  of  gas  Rate  of  gas 

flowing  out  produced  in  c.v. 


3_  n 

3t 


(1-v  )p  A  Ax 

(1-v  )p  A  U  1 
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It  Rl-v  )o  1  +  I—  f(l-v  )p  U  ]  =  p  (-7—^)  — - — — 

dt  L  s  gj  dx  L.  s  "g  gj  S  A  (1 

P  I  Pvgdx 


Using  equation  (1-3)  and  U  =  U  =  U: 

s  g 


3p  3p  ... 

(1"Vs)3Fa  +  (1"Vs)  U  3^  +  Pg  3x  =  (ps  “pg>  *d 


8pe  3pg  pe  3U  (ps“pe) 

- 6.  4.  n  - 6.  4.  - fi - -  - 2 - S_  y, 
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A  general  equation  of  continuity  can  be  obtained  by  considering  the  gas- 
solid  mixture  as  a  whole  which  gives, 

Rate  of  increase  of  Rate  of  gas-solid  Rate  of  gas-solid 
mass  of  gas-solid  -  mixture  flowing  -  mixture  flowing 
mixture  in  c.v.  in  out 
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Momentum  Equation. 


Rate  of  increase  of  Momentum  flux 

momentum  in  c.v.  out 


Momentum  flux 
in 


T,  External 
forces 
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Using  equation  (1-5) : 
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Now,  mixture  density  p^  = 
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Energy  Equation. 

Rate  of  increase 
of  energy  in  c.v. 


Energy  flux 
flowing  in 


Energy  flux  Rate  of  increase  of 
flowing  out  energy  due  to  con¬ 
version  of  solids 
into  gases  in  c.v. 


Rate  of  work  done 
by  the  gas-solid  mixture 


Rate  of  heat  loss 
to  the  tube  wall 


or, 
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where,  AE  =  Additional  energy  release  per  unit  mass  due  to  conversion 
of  solids  into  gases, 
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Using  the  definition  of  enthalpy  for  solids  and  gas,  i.e. 
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Using  U  *  U  =  U , 
s  g 
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Using  the  general  equation  of  continuity  (1-4)  and  the  momentum  equation 
(1-5): 
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Ou  differentiation  and  by  use  of  continuity  of  solids  (164)  and  continuity 
of  gases  (1-3):' 
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Case  II  (Solid  particles  stationary  at  their  initial  positions) 
Continuity  of  Solids. 

Rate  of  increase  of  _  Ra^e  of  solids  Rate  of  solids  Rate  of  gas 

solid  mass  in  c.v.  flowing  in  flowing  out  produced  in 
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Continuity  of  Gases, 


Rate  of  increase  of 
mass  of  gas  in  c.v. 


Rate  of  gas  Rate  of  gas  Rate  of  gas 

flowing  in  flowing  out  produced  in  c.v. 
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Putting  U  =  U  and  using  equation  (1-8): 
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Momentum  Equation.  As  the  solids  are  at  rest,  the  free  as  volume  in 
the  control  volume  shown  in  Figure  37  is  taken  as  the  new  control  volume 
in  the  following  derivation.  It  is  assumed  that  the  solids  are  at  the 
core  of  the  flow  and  the  skin  friction  at  the  surface  of  the  solid 
particles  is  negligible. 

Rate  of  increase  Momentum  _  Momentum  _  ^External 

of  momentum  in  c.v.  flux  out  flux  in  forces 
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Using  equation  (1-9)  and  putting  l  =  U, 
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luation.  The  same  control  volume  as  used  in  the  derivation  of 


momentum  equation  is  taken. 

Rate  of  energy  _  Energy  flux  _  Energy  flux  +  Rate  of  energy  increase 
increaoe  in  c.v.  flowing  in  flowing  out  due  to  gas  coming  into 

the  c.v.  from  the  solids 

Rate  of  work  Rate  of  heat 
done  by  the  -  loss  to  the 
flowing  gas  tube  wall 
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mm*** *■ 


f<Y’£T^1!«"»A"  "AW  W  u*r  ' 


3_ 

at. 


u  'I 

a-Vs)ApPg(eg  +-|-)AxJ  = 


L 


1 

(1-v  )A  p  U  (e  +  -S-) 
s  P  g  8  g  2  J 


-  (1-v  )A  p  U  (e  +  -4-) 


s  P  g  g  g 


x+dx 


v  Ax  r, 

+ps  sb  ^ ^  W 


t  P 


v  dx 
s 


f(l-v  )A  d  U  (~)l 
t  s  P’  g  8  P  j 


x+dx 


j(l-v  )A  p  U  (— )  ) 

l  s  P  8  8  V  1 


-  2ttR  Ax  h  (T-T  .  i 
i  w,i' 


xJ 


or. 


3_ 

3t 


U 


(1-v  )p  (e  +  -*-)  ,  + 


1  ,A_ 


s  g 


8  T  2  'J  *  * 


U 


(1-v  )p  It  (e  +  —  +  -1 ) 
s  g  g  g  P  2  ' 
8 


b  v  r  2h , 

P..(t^>  TT^ -  W  -  ~±  (T-T.  ) 

IV  w  |  1 


s  A 


P  /  P  ^dx 


Put  U  “  U  and  li  ■  e  +  —  , 
8  8  go 

g 


JL 

at 


(l-vs)pg(hg-I-) 


+  «  L(1"VV'gLJ  *  4  7F  [(1’vs)p/r?j 


+  >.  JL 

1  3x 


{(1-v  )p  U}U‘ 
s  8 


p  v  W 
s  u 

s 


2hi 

—  (T-Tw,i) 
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Using  equations  (1-9)  and  (1-11) 

|~(l-v  )p  h  "I  +|-  r(l-v)phu!  -  |- 
3t  I  s/Kg  g  3x  s"g  g  j  31 


(i-vs)p  +  Ps  vd  2- 

J  s 


+  U  -p  v .  U 
s  a 

L  s 


/1  %  3P 

(1-v  )  — 
S  3x 


—  -  p  V ,  W 

R  J  s  d 
J  s 


2h 

-—(T-T  .) 

R  w,i 


Using  equation  (1-8) 


—  j  (1-v  )p  h  1  +  ~~  [~(l-v  )p  h  U  |  -  (1-v  )  ("H  +  U 

3t  [_  s/Fg  gj  3x  s  Kg  g  J  s  [_3t  3x 


s  3t  3x 


P  u2  2hi  2l»u 

*s  <“ -  r +  f>  -  -T  (T-T»,i)  +  T- 

s  Ks 


f  3h  3h 

(1-V  )p  Ir^+U^ 
s  gl_3t  3x 


’h*[tt  {<1-V’g)  +  fe{‘1-VV’ 

-  (J  -v  )  j"~  +  U  1^1  -  p  v  (W  +  --  +  --) 
s  [_3t  ox  J  s  dg  pg  2 


2h  2t 

_r<T'Tw,i>  +  "T u 


D  3  $ 

Using  equation  (1-9)  and  notation  —  =  +  U  ,  finally  : 
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(1~vs)pg5^  "  (1“V  m  =  PS  *d .  (w  +  r  +  T  ~  V 

w5  S 


2h.  2x 

— -  (T-T  .  )  + 

R  w,i  R 


(1-12) 


Computation  of  burning  Surface 

A  typical  solid  particle,  a  single  perforated  cylinder  in  shape, 
is  shown  in  Figure  38. 

Let, 

=  initial  inner  radius  of  the  particle 
rQ  =  initial  outer  radius  of  the  particle 
1  =  length  of  the  particle 

n  =  total  number  of  the  particles  in  the  chamber 
Therefore,  total  initial  burning  surface  »  2TT(r^+rQ)  In.  It  is 
assumed  that  combustion  gas  is  produced  from  both  inner  and  outer  cylin¬ 
drical  surfaces  of  a  particle  but  not  from  two  ends.  If  r^  is  the  linear 
speed  of  burning,  which  is  assumed  to  be  same  for  all  the  particles  at 
a  particular  instant,  the  total  burning  surface  after  time  At  is 

2n  ^(r^+r^  At)  +  (rQ-rb  At)J  1  n 

=  2tt  (r.+r  )1  n 
i  o 

It  is  clear  that  for  hollow  cylindrical  particles  the  total 
burning  surface  is  constant  and  can  be  given  by  •' 
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2  2 

2tt  <r  -r  )1  n  p 
o  l _ _2_ 

p  (r  -r . ) 
s  o  i 


where 


m  =  initial  mass  of  solids 
s . 

l 

p  —  solid  mass  density 
s 

w  =  initial  web  thickness  of  a  solid  particle 
s . 

l 

Expressions  for  Enthalpies  of  Solids  and  Gases 
lor  any  pure  substance,  h  =  h(P,T) 


(1-13) 


dh  -  (§>  OP  +  (f)  01 

T  p 


(1-14) 


From  thermodynamics,  dh  =  Tds  +  vdP 


,  J**)  +  v 
,3P  ,’T  wp;t 


(1-15) 


Again  from  Gibb's  form  of  first  law  of  thermodynamics, 


dG  a  -  s  (JT  +  v  dP 


(1-16) 


As  Gi 


Gibb’s  function  G  is  a  property  of  the  system, dG  must  be  an  exact 


differential 
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How,  coefficient  of  thermal  expansion, 
..  Equation  (1-15)  becomes, 


c  -  i/9v) 

S  '  vl3T/_ 


-  T8v  +  v  =  v(l-T$) 


Equation  (1-14)  becomes, 


dh  =  v(l-T8)  dP  +  c  dT 

P 


(1-17) 


(1-18) 


Solids 

For  solids,  c  is  equal  to  c  and  it  has  been  assumed  that  the 
P  s 

temperature  of  the  solids  T  remains  constant  throughout  the  period  of 

s 

burning 

•  dh  -  --  (1-T  3  )  dP 

•  •  s  p  s  s 


It  has  also  been  assumed  that  the  coefficient  of  expansion  for  solids 
Bg  is  negligible. 


(1-19) 


and  h 

s 


c  T 
s  s 


(1-20) 
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Gases 


The  equation  of  state  for  the  combustion  gas  at  high  pressure 
can  be  taken  as , 


P(v  -n)  =  R  T 
g  g 

where,  the  covclume  n  is  a  constant. 
Now, 


From  equation  (1-18) 


dh 

g 


r  R  T  1 

v  1 - £-  dP  +  c  dT 

g  L  yj  p 


d-21) 


or, 


dh  =■  1  d?  +  c  dT 
g  P 


(1-22) 


Now,  the  specific  heat  at  constant  pressure,  c  =  -~r  R  yields, 

P  Y-l  g 


dh  «  n  dP  +  — '—r  R  dT 
g  Y-l  g 


(1-23) 


Differentiating  equation  (1-21)  and  using  v  *  1/p  : 

8  8 

R  dT  =  (— - n)  dP  -  — — r  dp 

g  Pg  Pg2  g 
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Equation  (1-23)  becomes, 


(Y-nP„) 

>S~-  dP 


dhg  ~  Pg(Y-l) 


-I* 

(y-Dp7 


dp. 


From  equation  (1-22) 


h  =  n  p  +  c  T 
g  P 


=  n  P  + 


.1— 

(Y-l) 


P 


(Y-nPg) 

(Y-l)pg 


P 


(1-24) 


(1-25) 
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APPENDIX  II 


DERIVATION  OF  BOUNDARY  LAYER  MOfffiNTUM  EQUATION 
With  reference  to  Figure  39,  the  boundary  layer  momentum  integral 
equation  is  derived  as  follows: 

Equation  of  Continuity 


Rate  of  increase 
of  mass  in  c.v. 


Rate  of  mass  flow 
=  into  the  c.v. 

through  surface  AB 


a__ 

at 


J 


r 

I 


2-rrr  Ax  p  dr 
R-6 


Rate  of  mass  flow 
out  through  surface 
CD 

Rate  of  mass  flow 
into  the  c.v. 
through  surface  BC 


x+dx 


+  m. 


BC 


CII-1) 


Momentum  Equation  (x-directional) 

Rate  of  increase  Momentum  Momentum  '^T  External 
of  momentum  in  -  flux  in  +  flux  out  =  forces 
c.v. 


9_ 

at 


rR 

"  rR  2 

"  rR 

I  2irr  Ax  p  u  dr  -  mD_  U  - 

BC  *° 

1 2nr  p  u  dr 

1 

+ 

|  2rrr  p  u4r 

R-6 

-  'R-6  -J 

X 

] 

<0 

1 

_ 1 

x+ctx 


r  f 

|  j  2 nr  P  dr 

L  4-6 


x+dx 


+  2tt(R-6)(— ■  Ax)P  -  2ttR  Axt 

OX  w 


( 1 1  -  2 ) 
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Using  the  expression  for  m  from  equation  (II-l),  equation  ( 1 1 - 2 ) 


becomes , 


(pu  r)dr  "  UK  |  (pr)dr  (pu  r)dr  +  (pu2r)dr 


-  -  I  (P  r)dr  +  iR-6)P  — ■  -  t  R 
3x  dx  w 


( 1 1  -  3 ) 


Now,  for  thin  boundary  layer,  radial  component  of  velocity  is  very 
small  and  consequently, 


Also,  as  U  i*  f  (r) , 


(II -4  D 


it  I  u«pr  dr  =  u»it*  pr  dr  +  pr  dr  ir 


(II-5) 


r 2  f R  rR  su« 

h  I  u»(pur)dr  =  u~  h  pur  dr  +  j pur  dr  9T 

'R-6  4-6  L/r-5 

Using  relations  (1 1—4)  ,  (II-5),  and  (1 1-6) , finally 


( 1 1  ~  6 ) 


rR  rR  r R  9ua 

P(Uoo“u)r  dr  +  ~  P“<Vu)r  dr  -  pr  dr  ^ 


'R  >.  [f  1  9P 

pur  dr  - —  =  r  dr  —  +  x  R 
^  3x  I  3x  w 

1-6  J  L4-6  J 


(H*7) 
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now. 


-R 


r  dr  =  R6  - 


R-6 


and  using  p^,  i.e.  gas  density  corresponding  to  the  film  temperature, 
as  the  average  density  in  the  boundary  layer,  the  equation  (II-7)  becomes, 


3t 


p(U  -u)r  dr  +  —  pu(U  -u)r  dr  + 


'R-6 


3x 

'R-6 

,2 

(R6  -  y-) 


R 

p(Uro-u)r  dr 
Mi-6 


3P 


3U 


3U 


3x  +  Pf  3t  +  pfU°°  3x 


3U 

a 

Bx 


+T  R 
w 


(II-8) 


The  equation  (II-8)  is  the  required  momentum  integral  equation  for  a 
nonsteady,  nonuniform,  and  developing  compressible  flow  in  a  tube. 
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APPENDIX  III 


FLOW  CHART  FOR  THE  SINGLE  SHOT  COMPUTER  PROGRAM 
ram 
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’  F  * 


‘-No 
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Subroutine  BP 

The  purpose  of  the  subroutine  is  to  calculate  the  ballistic 
properties  at  the  piston  base  at  time  t’+At’  knowing  the  properties 
at  all  the  points  at  time  t',  and  the  velocity  and  the  position  of 
the  piston  at  time  t'+At*.  The  following  flow  diagram  should  be 
read  along  with  Figure  3. 
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Subroutine  WP 


This  subroutine  is  called  to  calculate  the  properties  at  the 
tube  head  end  at  time  t'+At'.  The  logic  is  same  as  that  for  subroutine 
BP  and  the  flow  diagram  should  also  be  read  along  with  Figure  3. 


Calculate  the  temperature  at 
point  2’  using  the  equation  ot 

state  (6) 


Return 
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Subroutine  HTW 


The  purpose  of  the  s ubroutine  is  to  calculate  the  wall 

temperature  at  a  particular  station  along  the  length  of  the  tube 

at  time  t+At,  knowing  the  temperature  distribution  at  the  present 

cime  t,  and  the  mean  heat  transfer  coefficient  h  over  the  time 

m 

interval  At.  In  the  following  flow  diagram,  At  is  the  time  step 
selected  in  the  main  program  whereas  At  is  the  time  step  selected 
in  accordance  with  the  stability  condition  for  the  temperature  calcu¬ 
lation. 
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APPENDIX  IV 


PR 

BUR 

TUBEL 

XNI 

PI 

DIAT 

HASSP 

P0 

T0 

RSP 

WM 

ETA 

DENS 

CHAR0 

BETAS 

WEBT 

MN 

CPC, 

C0NG 

H 


SINGLE  SHOT  COMPUTER  PROGRAM  INPUT 

Input  Parameters 

Pressure  at  which  experimental  burning 
rates  are  input  (nt/m^) 

Burning  rate  at  PR  pressure  (m/sec) 

Tube  Length,  also  reference  length  (m) 

Nondimensional  projectile  starting 
position 

Nondimensional  axial  tube  step  size 

Tube  diameter  (m) 

Projectile,  mass  (kg) 

Shot  start  pressure  (nt/m^) 

Propellant  flame  temperature  (°K) 

Ratio  of  gas  specific  heats 

Propellant  molecular  weight 

3 

Co-volume  of  gases  (m  /kg) 

Propellant  solid  density  (kg/nr*) 

Initial  propellant  mass  (kg) 

Zero 

Propellant  web  thickness  (m) 

Gas  viscosity  at  T0 

Specific  heat  of  gas  (kcal/kg  -°K) 

Thermal  conductivity  of  gas  at  T0 
(kcal/m  -sec-  ’K) 

Boundary  layer  velocity  profile  parameter, 
for  one  seventh  profile  =  1.2857 


212 


B 

Ludwig  -  Tillman  friction  factor  parameter 
=  0.286 

R0UT 

- 

Tube  outside  radius  (m) 

RIN 

- 

Tube  inside  radius  (m) 

DELR 

- 

Tube  wall  radial  step  size  for  numerical 
conduction  solution  (m) 

DIFFU 

- 

Thermal  diffusivity  of  metal  (m^/sec) 

C0NDUC 

- 

Metal  thermal  conductivity  (kcal/m  -  sec  -  °K) 

TAMB 

- 

Ambient  air  temperature  (°K) 

H0 

- 

Heat  transfer  coefficient  for  tube  outside 
surface  (kcal/m^  -  sec) 

UNI 

- 

Zero 

UNBI 

- 

Zero 

PRNI 

- 

One 

PRNBI 

- 

One 

DENNI 

- 

One 

DENNBI 

- 

One 

TEMNI 

- 

One 

TEMNI 

__ 

One 

All  data  Is  to  be.  input  in  format  form  to  conform  with  read  and 
format  statements  as  follows.  An  example  is  also  shown. 
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READ  199,  (PR(IJ) ) t  BUR(IJ),  IJ  =  1,20) 

199  FORMAT  (4(E12.2,  F8.j,) 

READ  20C,  TUBE2,  XNj ,  PI,  DIAT,  MASSP 

200  FORMAT  (5F16.4) 

READ  201,  P0,  T0,  RSP,  WM,  ETA 

201  F0RMAT  (E16.3,  F16.1,  F16.4,  F16.2,  F16.6) 

READ  202,  DENS,  CHAR0,  BETAS,  WEBT 

202  FORMAT  (F10.1,  2F10.3,  F10.2) 

READ  203,  MU,  CPG,  C0NG,  H,  B 

203  F0RMAT  (5F16.7) 

READ  204,  R0UT,  DELR,  DIFFU,  C0NDIJC ,  TAMB.H0 

204  F0RMAT  (6F13.7) 

READ  101,  UNi(l) ,  UNBI,  PRNI(l),  PRNB1,  DENNl(l), 
DENNBI,  TEMNl(l) ,  TEMNB1 
101  F0RMAT  (8F9.2) 
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APPENDIX  V 


SINGLE  SHOT  COMPUTER  PROGRAM  LISTING 

DIMENSION  UNH2W1  »  .PRN1  ( 2*>  1 )  » DENN 1  (  2 ^  1 )  »  T£ '-'N 1  ( 20  ]  )  ,  VpS  \  (  2^1 )  • 
lPRi^0X»ItURt2'-  1.*  J 1 12 i-i  1 )  *C  1  (2U14-.CM1  1 2  Oil. UN CAR  (20 1L.UN2  L2014-.- 
2PRN2I2U1 ) »DtKN2(201 1 ,1EVN2(201 ) »VFS2 ( 201 ) *U?E (201 ) »PR2E<  201 ) . 
3uEN2E(2£-14»Vf2E(2L14-.L)2£(2ul4  »-C2L(201ijOM2lll2CU-*.L!2.E2.(2QU-» 

4PR2E2  (  201 )  »DLa2£2(201  1  .VF2E2(  201 1  * ; i T <  1  (  20 1  )  ,MT2E  ( 201 )  .'-1T2E2  ( 201 1  . 


-Mil  i  .A?-T>:t?ul  <  till  (  2lU  1 « H  ?  I  /  0  1  1  .  TEM.M12.1  L*S11.»,  REU-L20LL* 


6CH2ul)»E2E(2i-l>  »i>CNF  1(201)  .DENF2 (201 1  .OcN2F(2Gl 1  »DEn2F2<  201 )  * 
7RIMF.1 2(114. 

DIMENSION  HP(4.50J) .TG(4»500) 

-XQMIdON-  yUCOLF /OENNS-tJ  bETAj - ELLA.V. 

COMMON  /  M  B  P  W  P  C  /  R  S  P  .  E  T  A  D  0 

COMMON  /  V.nP-.P  /DTiM  «  P  I  .  PR.Ml  .11N.Ej.D2MU.j-V  FS 1  i  QSV  .  2UMft  iilal » 
1H1»REMT  *B*DENF1 


COMMON  /MftP  /PRNB1..D&NNB  l-tUNEl  ».UNa2->XN-LtXC12. - V-E-Sai 

COMMON  /MRAUT /T  IMEP  (  500  )  ,TWXl  (600  »TWX2(600)  ,TVX3(600)  .T'.VX4(600>  » 
!Illl£5EX£a-»JKiAXtX^LSLl_l.»2mi£L*A^ 

COMMON  /MBPR/KT 

COMMON  /MHTf/NN.i  .  R I  N  .  ROUT  tJ3f4_H-mlll.F  F  Llt£HNIlLlf  «  H I  » HO  I  T/VM11  »iIL 


coMMor;  /mhbpkp/to.temni  *temm 

_  RE  AL.M  AS  Sp-.K  TK1  »MT2E-tHJ 
199  FORMAT (4(E12»2*F8»5) 1 
2A10.  EQRMAJLl  AE16..-4J 

201  FORMAT ( E 16. 3. FI  6. 1. F 16. 4, F 16. 2. F 16. 6) 
2112  f  ORMAT  (F10.1  .?F  11). 3  »F  1  0 .7  I 


203  FORMAT ( 5F16 • 7 ) 
2(14 _ FORMAT  16  F 13. 74- 
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211L 


FORMAT (17H 


MOLECULAR 


WT«»F6. 1 

..  1  t)  OHM 


2?  H 


.Xhl 


RATIO  OF 
=  >F5»4jSH 


1 


MASSP*  .F5.3.13H  ICG 


DIATs,F5.4,llHM 
20, OH  PRANDIL-.~«F6»4t/ L 
211  FORMAT ( 1 1H  PO-.E8.4.23HNEWTON 

1FH  CU.M  TO  a  «F6« 1 .  15HDEG — X - UN.LI.2: 


SP.HEATS«.F6.4./1 
P  I  =  .  F5._4_»4_lld — 
COVOL*  *F8 •  » * 7HCU.M/A 


PER  SQM  DENGO*  »F6 . 2  * 17HXG  P 
« Ff»  .  1  « 1  9HM  PER  &EC _ TIMEC*.. 


2F7.5.3HSEC/1 
21?  FORMAT  1  12H— 


CHARGE-=-*E6.»  3 12 IH-KjG- 


*Eri  Thi  LCJLNESSjljE-Sj^S  1 2 .1 H  M 


P0TENTIAL*»F1C,2.3H  KCAL/KG/ 


1  SOLID  DENSITY*. F7. 1 .22H  XG/CU.M 

It  1 

218  FORMAT ( 24H  G. V ISCOS I T Y ( REF .TEMP  1 « »F8 .6 . 17HKG/M-SEC-  GAS  CP»iF6.4* 

219  FORMAT ( 13H  KETL.DIFFU*.F10.8.20HSQM/SEC  METAL  COND»»F8.6.23HXCAL 
l/M-SFf~X  AMB.TEMPa.F6 . 1 1 3npF 0  1 1 \ 

213  FORMAT  (//10H  T IME ( ND 1 «  »F9 • 6 » 15H  SOL  3Y  MASS* . F8 .6 , 14H  AV. PRESS 

HIBFa.FlU.A.lHH  PISTON  on  IND)a.F9.6a?il - HE.  AT  LQS.WE  'WvaXl - 

214  FORMAT  ( 115H  D 1ST ( ND )  VEL(ND)  PRtSS(ND)  SOL  BY  VOL  G4$  D(ND) 


1  mxjaiADJ  Tf  MP(ND) _ THICK  ( ND  )  REYNOLD  NO  H.JL.CO££ - IV.1  ALLA 


215  FORMAT  (7F11.6) 

>16  FORMAT  17fll.fe.Fll. 7 »2E1  U.AtF-111.21. 

101  FORMAT ( 8F9.2 ) 

77  FORMAT  I 3GI 6 . 9 1 

250  FORMAT ( /7H  TIME*.F8.6*13HSEC  PR0.VEL*.F8.1 » 15HM/SEC  PRO.X.E*. 
1FB.2.20HXCAL  INPUT  ENERGY*  ,F8. 2 .22H<CAL _ SAL  L  I  ST.UL-LF-F.I  *  .  Ffl  t  U. 
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FORMAT (5F15.7) 

READ  199, (PR .( JLJJj  BUR  (  1  J  ) .  I  J«  1 1 20  ) 

READ  200.  TUoEL.XNl.Pl .D1AT.MASSP 

READ  201.  PO.TO.RSP.WM.ETA 

READ  202.  DENS. CHARO. BETAS. WEST 

READ  203*  MU.CPG.CONG.H.B 

READ  204.  ROUT. DELR.DIFFU.CONDUC. IAMB. HO 

READ  101 .UNI ( 1 ) .UN61 »PKN 1 ( 1 ) .PRNBl .DENN1 ( 1 » .DENNB1 . TFMN1 I 1 ) .TEMNB1 
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AREAd=2*CHAR0/  <0E«S«W£*iT4- 
N=XN1/P1+0.0001 
AREAP=3.14159*(01AT**2)/A« 
RU=03!4./WM 


TJM£C=TUBEL/UN1  TV 

G£N60= 1 ./  ( tT  A.tRU*m/PO  J _ 

VT=ARE AP*XMl*TUd£L 
VFS1=JCha.po/vt-OLNGQ>/(OENS-D£NGO> 
00  500  1=1, N 

< ! 1 sUPA 


VFSU1=VFSI 

0  EN  N  S= D£_S1  S/.  DiLMGjQ 

ETADO=ETA*DENGO 

I£aS=3Q0-. _ 

T8ETA=TEMS*BETAS 


HVS=RU*TO/(  (RSP-1>*4184.) 


CALV=HVS 

Em=DJLAIZ2-.- 

NN= (ROUT -RIN)/DELR+u. 000001 


00  251  1=1,201 

i  i=i  . 


251  TEMM( I ,J)=TAMB 

=  i  . 


MTK1( I  1=0. 


■■till 


REMT ( I ) =0, 

111  t  T  1  = 


QT  =  0. 


A»U#123/( 1U*C**(0*678*H) ) 

/  i  .  /  1 1  " 


BUM1*<  (1+B)*A/H)*(TUBEL/RIN)*BUM 


IKii  iii.mm  u«nM3kM(i«n  nniiiiirniiui  a« 


WRITE (6, 21 2 )  CHARO.WEBT ,DENS,CALV 

uDtTci  a  .  noivn.ror.  .rnwr.  .  lj 


WRITE (6.21 9 )0!FFU,C0NDUC»TAMB 
jao.  1 . JX«2,H 

UNinX)«UNl(l)  +  (UNBl-UNl(in*UX-l)/N 


DENNK  IX  )*DcNNl ( l )  +  ( DENNiil-DENNl  (1))*(IX-1)/N 


M-N 


ua.1  ■■iE3ifl3ui;iiicau^:iiMBi.iii  dimwit 


TlMEW-0. 

X  T«  1 

TIMES-0.2 

1K»0 

CALCULATION  OF  AVERAGE  CHAMBER  PRESSURE  AND  BURNING  RATE 
5  PDUM1«PRN1(1) 
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107  VDUM1=VDUM1+VFS1 C I ) 

_rsvH=v£um*p  l±vf.siuu«  Ix:a--U‘--Xjjscj4 
V0UH2=VFS1C2) 

_ DO-  108  _1  =  3-»H 

108  VDU«2=VDUM2+VFS1( 1) 


ikvi  ■R'rtii  iz&jliim  tv  4  a»T>iu>i;iPijg  ■  kj’Jib  \ 


TSV=(TSVH+TSVL»/2. 


I [•  1  yj a »T wtCW  J *  1  ‘ijl ipi a U fcvfi ilLVii ■ 


IF  (M.LT.  13 )  GO  TO  73 


GO  TO  72 


i  ■ nw^  «■¥=!■  cult  si»i 


72  HP( 2  »KT )=H1 ( 13 ) 


iriBid 


TG( l.KT I=TEKN1 (26)*TO 
r,n  rn  7m 


700  TG(1,KTI=TEKN31*T0 


DO  623  IJ=lt51 


>  in  i  a  ■  m  m  1*3  ■  a  *  w  ■ 


623  XWMCIJ)*UJ-1)*0.02 

i  i  if  i 


TIMES=TIMES+C«2 


371  WRITE(6t213)  T 1ME1 *PR0PM* AVPRN*UNB1 »QT 


DO  90  IX«ltM.5 


DM1 ( 1X)*VFS1< IX)«DENNS-M  1-VFS1 IIX >>*DENN1  IIX > 


^fR  I  TE  <  6,216)  X.UNK  IX)  .PRNl  (IX).VFSI  <  IX  I  .DENNK  IX)  »DM1  IIXI  » 
1TEMN1 ( I X ) * AMTK ( IX ) »REMT (IX)tH 1(1X1 *TEMM( IX  *1 ) 

GO  TO  90 

89  WRITE(6.215>X,UNl(lX)»PRNlUX)»VFSUIXI»DENNl!IX)tDMlUX)» 
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WEMNltlXI- 
91  CONTINUE 

OMiil=VFSbl*OfNNS+ll-VFSUU*OENN.-U 

WRIT t »  6 . 2 1 S  >  X  N 1 .  1 ,  PRfJd  1 .  VF brt  1  .OENfii) 1 ,0'M  1 ,  TEKN3  1 

-I  Li  Li:  ?J  liiEXtll.  J 2h 

C  SCI. EC  T  1  ON  OF  NEW  TLir.  I  NCWCAMENT 
.  -  37a-DQ„.2_L=a»M 

RTMFII  )=2*70*TfcVU<  I  )  /  (  TO»TEMNl  t  I  l  +  TF.MMI  1*1)) 

_ OtNF.l  1 1 )  sRTKF  ( I )  *  JE1N1U  12  1 1.  ♦•EI4,>0*l>C*4;ilLlJL*  (K.TMF  U.j  L) 

CALC  COEFIFRM  !  1  )  iDENNl  (  I  )  ,VFS1 1  I  )  .HI  (  I  )  ,C1  <  I  ),D-.i!  I  )  .DUM1  ,DUM2, 
IfiLUL 

2  UNCAR  I  I  )=J.\1  (  I  )  +SQRT  I  B1  (  I  J/OM1  II)) 

..CALL  CQEF  iPRNiil  .UCNNiL.YESDl  »Bdl..»CBL^DMtll-»IlUiVlJDUM2»EHL^L 
0NCAKt3=UNI31+SCRT  UiU/DMtil  ) 

.UMAX.*  UNCAR.1  1J 
00  12  1*2. M 

11  UMAX=UNCARU) 

12-COMUNUE 

IF  <UMAX~UNC4?3ll3,l4.14 
12-UMAXsUNCARH. 

14  DTN=P I /UMAX 

MaOINftl.UBELV.ULlIU 

KP=1 

.AVgRl*AVgRN*PQ. 

AVPR2*AV1R1 

120,  AVPR=lAV?Rl*AVPR2)/2. 

00  17  IJ«1»20 
IF  ( AVPR-PR (1J))1S«18*17 
17  CONTINUE 

lfl.  SURD*  BUR  I  I.J-.l  L±i3URLLallsflllRl  IJ.rrll  )»1AV£R-PALI  J-lD/IPP.Ll  H- 
PR1IJ-H) 
flUR&sflURO/UNlIV 
IF1T3V.LE.0.)  GO  TO  360 
HSltafl  AI  IflA.*allRfi2 I S  V 
GO  TO  361 

360  DSVaU.  * 

C  CALCULATION  OF  NEW  PISTON  POSITION, ITS  VELOCITY  AND  GAS  PROPS  AT  BASE 

.161  .ACCl»P0«PRNBlM3&Aez*&S.S£ 

0ELS»UN81*UNlTV*DT+ACCL*(DT»*2)/2. 

QE1  SHyDF-1  S/lUSfi 
6  XN2-XN1+0ELSN 

CALL  BO<PB.DL0,TEU,VSB> 

-PiiEW  \l  RRMLtEaJJ2ji 

ACN<-W«PO»PNEW*AREAP/MASSP 

ttiEaL"  imamua  ly«&UACNEW*-l  0T««2l  Z2a/IUBEt 

IFC (ABS (ONE W-DELSNI/OELSN) -0.000 1)20.20,25 

2i_QELSU-DaEW~ 

ACCL-ACNEW 

ift  TO  .6 

20  VN2-XNUDNEW 

i)Nfl?»tJN3)  +ACh!FM*DT /UN  I  TO 

PANB2-P9 

DF>'Nfl?«DFR 

TEMN32»TEB 

VFSfl2«VSB 

C  CALCULATION  OF  GAS  PROPERTIES  AT  TUBE  HEAD  END 


MTK2 1 1 ) *0. 

CALL  WP IPRN2 ( 1 1 .OENN2 ( 1 ) . TEMN2 ( l ) ,VFS2 (II) 

C  CALCULATION  FOR  INTERIOR  POINTS  START 

U2E( 1 1 *UN1 ( 1)-DTN*UN1 ( 1 )*(UN1 ( 2 1 -UNI  111  I /Pl-DTN* (PRN1 ( 2 l-PRNl (1)1/ 
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-UCKlt  1  >»PIl 

DEN2£  (  1  KDENN1  ( 1 )  -DTN*UN  1  (  1) *  ( DENN1  ( 2  ) -DEN‘11  ( 1 )  J  /P I  -OTN+DENN1  ( 1 )  * 

1, ( UN  1  f  2  )  r.UNU  U  )  /  li  1-ViLSli  X-JJ  »PU+.UJN*  !  OL'NNS-OLNN  1  ( 1  )  )  »VE"1(1 )  *DSV2 
2  C 1-VFS1 *  1  * » 

PK  ?  F  I  1  1  =  P  R  N  1 . 1 . 1  >  —  i  *  T  V  ♦  1 1 N  1  I  1  HIPRM1  1  2_Lr.tL2N.l_L.LL)  /l>I-l)TN«:U  l  I  I  »LUhlI2-L=- 
1UNH1  n/PI+CICl  )*USV*VFSK  1  )»OTN 

_  VF.2E-1.1  J=VFSlL(.U-DTN»UN1  LI )  *  ( Vf  SU2-)  +4LES1  i  1.L1/P  l-OJN*VESUJL)  « 

1 (UNI (2  )-UNl ( 1 I )/Pi-DTN*DSV*VFSl  ( 1 > 

_IF_  t VF2  H  U .»  L L .u  . J-  VF-2f.  (.  1 4  =  Q •- 

CALL  COEF ( PR2t ( 1 >  «DfcN2E< 1 > *VF2fc ( 1»  »B2E(1 J  »C2E ( 1 J »DM2F(1) *DUM1 » 
]pllf.<?.F  7  V  LU-X 
MT2EI 1 )=0. 

-X-L=XN.l  /£  L 
J=Xl+0. GOOOl 

_L£.I  Ai3S(  X±r-  J-L=0. .  QaiUl  15jO-»5jj.»  5J. 

■50  UNK  J  +  l  )*UNB1 

DENN1 ( J  +  l ) *DENNB1 
JLESUJiLUVESai 
MTK1 ( J+l ) *0. 

HI ( J+l )«0> 

RTMF  (  J+l >  =2»T0*TEMNB1 /  (  T0*TF.MNB1+T  AM6 > 

QFNF1  (J  +  l  i  gRT.MF  ( J+1  )  *DENM1J  J+ 1  )  1 1,1  t.t£TADO*DF.NNl..l  J.tl.Lt.i£TMF  t  JiJL ).=-L«-. 
1 1 ) 

51  DQ  30  1  X«2  «  J- . 

MULSPU*$QRT ( ( T0*TEMN1 ( I X l+TEMM ( IX  *  1 « )  / ( 2*T0) ) 

,M9»HULA*B _ 

I F C REMT I IXJ.LE.0.)  GO  TO  6000 

III  nvsFaHiiMfe+DFNFi  ( lx  i in i  t  tx  )+»?> UHUU  t  l.xut-U*£Ml (1  X  i (tmm 
PGA1NF«3UM6*E1(IX)  0ENF1 ( IX) * ( UNI ( IX > **3 . > / ( REMT < I X ) **B ) 

GO  TO  b£QX 

6000  ULOSSF«0« 

PGAIMF«Qj^. 

6001  U2E(lX)«UNlUX)-0TM*UfaUX)»<UNl(lX+l)-UNl(IX))/P!-DTN*(PRNl(  IX+1J 
i-ppmi  i  mwinMimuPii-ui  n£SF»nTN 

IF  (U2E ( IX) *LE • 0. )  U2E ( I X ) -0* 

nFMPF I  lx  I »DFNM1  LLXJ -flTN*lJNl  1 1 X  1  +  ( DENN1UX.+JLUOENN1  Ux)  i/.21dimt 
1DENN1 1 IX  I* (UN1( IX  +  1 1-UNI i IX) ) / ( ( l-VFSK 1X5 i*PI )+DTN*(DENNS- 
2CENN1  UXXl«VFSinx>*r)SV/ll-VFSltUU 
DEN2F(  1X»-RT.MF(  IX  )*0EN2E  <  I X )  / 1 1+ETA00*DEN2E  f  I X  I*  <  RTMF  ( IXI  —  lJ  I 
DB|Fim.pPMi  i  Tyi-nTM»HMi  I  nitiBB.nl  MX  +  1  I-P3N1 I  1 ¥1  I /P 1  -flTNiRLI IXX* 
1  (UNI  I  IX  +  1)  -UNI  (  IX)  5/Pl+CK  IX)  *VFSi  ( I X  )*DSV*D7N-DTN*BUM7*”1  (IX)* 
;nii|yi»iTniTFHM1im.TFvWI  I X  » 1  1 1+dLNAPGALNE. 

VF2E( IX)»VFS1( IX)-0TN«UN1( IX)*(VFJ.1( IX+1 >-VFSl( IX) )/Pl-DTN# 

ji/f$ii  lxitmNiiix+i  i-iiNi  i  mwPirnTN»n.sv»VFM  uxt 
IF  (VF2EIIXI.LE.0.)  VF2E(IX)«0. 

f*i  i  rr>FFiPP?cf  im.n»nm»i  .\zF7FHiM  .ft?F  1 1 x  I  .f  ?F 1 1 X  I  *nM2EX!XU. 
1DUM1 *OUM2»E2E( IX) ) 

U2E24-U-UUM1 >.nTN<iJ?EJJjiiAlii?F  I IX1-U2E  (  TV-Ill  /PI  -HTNi(PR?F!  lXl-~- 
1PR2E( I X— 1 ) )/(0M2t ( IX)*PI)-UL0SSF*0TN 

DEH2EZI  lxi»HENNl  (LX  i“QTN»ll?F(  IXL81DEH2E1 1XI-OFN2E1  UULl  )  /PI-DTNf 
10EN2E ( IX l»CU2Et IXI-U2EJ 1X-1 ) )/K l-VF2E( IX) )»PI )+0TN*(DENNS- 
2P£M2£ <  I X  l  ltVFAl  1 1 X  )  »QSV/ 1 1~VF?F LLXJ-L. 

OEN2F2( IX)"RTMF( 1X)«0EN2E2< IX ) /» l+ETAOO*DEN2E2< IX)*(RTMF( IX)-l) ) 
Pff?F?  I  LX La&BiU 1 1  y  * -r>TMtn?F 1 1 y  > »i PR2F I  lm -PR2F 1 1  x-1  ) )  /Pl-QItii. 

1B2E 1 1 X ) * (U2E ( IX »-U2E ( 1 X— 1 ) )/PI +C2E ( I X )»1>SV*VFSI t  IX ) *f>TN-DTN«BUM7t 
2E2EM»l»HKlX)»(TOtTEMNniX>-T£MM(IX.in+0m»g6AlN£- 
VF2E21 IX)»VFS1( lX)-OTN*U2E( lX)*IVF2E( IX >-VF2E( IX-1 > I/PI-DTN* 

1VF2E(  I X  )  +  (U2E ( IX) -U2E  (  I X~i  ) ) /P 1  -PTN«OSVVF.&m)LL 
UN2I IX ) *(U2E( IX)+U2E2 ( IX ) 1/2* 

!F(UN2(  m.LE.O.I  UN2(IX)>0. 

DEMN2I IXla(DEN2E( IX I+0EN2E2! IX) )/2« 

PRN2I IX )■ (PR2E ( IX)+PR2E2  CIXII/2* 
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V-F52  U  X  » M-VF2E4-I X  >  tV««UX-l  l  /£- 
IF  (VFS2<  IXWLE.O.  1  VFS2UX)=0. 

-3a-l£M.M2  l.JXisPJlIta  ( IX  1*1  AO 0  *  0  L NJ42  UX-U/-1-UEMN2UX4JU  UE-TADOl ) 

XP=XN2/P1 

X£X2 _ 

MM=J+1 

jm=ic+-i 

00  80  IX=MM*MN 

UN2  C 1 X  >  =UN2  (J )  +  ( UN02— UN.2  LJ  LI*  UX^l*£l/iJlh2=J  J=U*I>.U 
PRN2I IX»=PRN2<  J  ) +■  ( PRNo2-PRN2  <  J))*(  I X-J I  *P  I  /  ( XN2- <  J-l  )*PI  ) 


VFS2UX  )=VFS2(  J)-MVFSU2-VFS2(  JJ  )*(  I X- J ) *P  I / 1 XN2- ( J-l  >  *P I  ) 

.  I£_.l.Vf ji.2-UXI-.i-E.  C  .  1 _ Vf-SZULX.  1  =  G.*. 

80  TEMN2I  IX)=PRN2(  IX)  M  1-ETAOO*DENN2 <  I X  ))  /  ( DLNN2  ( I X  )  *  ( 1-ETADO  U 
J4=K±Jl 

IF  I XP»E0*2 )  GO  TO  1UU 
i?0UM3.-PR:i211J,.. 

KM*M-1 

DO.  lit) . I  =2-«KM. 

110  PDUM3=PDUM3+PRN2II > 

POUM4=PRI'42  ( 2 ) 

■00 . 

111  PDUM4=PDUM4+PRN2 ( I ) 

AVPR2*P0* ( AVPRH+AVPRL  Ml* 

JK£=2. 

, GO  TO  120 
1,00  DQ  31  IX-2.J. 

MUl_*MU*SURT ( ( T 0*TEMN 1 ( I  X ) +TEMK ( I  X  *  1 ) )/(2*T0)) 

Ha=»l’lUL*«.a 

IF  (UNlUX).LE.O.oOl  I  GO  TO  1200 
U»  LU2EI 1X1  :D«11 1XJJ7JJHLLLXJ 

UP*8.0*<U+DTN*MPRN1  <  IX  +  1 J-PRNU  IX  >  >  /  (DENF1 II X » *UN1 1  IX )  »PI  1  + 

ir  - 


GO  TO  1201 


] MB*DTN* (UNI ( I X I ** ( 1-B ) ) / IDENF1 ( IX »**B )-( l*fl)*MTKl (I  X } *( DTN* (I H+2 ) * 
7  man  m+1  wini  mn UhlB. U  +uni  i ix i»ldenfi  lix+xl-denfi  t  ixi  LOHiEl* 
3UENF 1 ! I X ) ) ) +U+ 1 0EN2 F I J X ) -DENF 1 1 1 X ) I/DENF 1 1 1 X  M ♦ 1 1+B ) «MT XI ( 1 X ) *UP 

_  If  1MT2ELLIXJ  «JLE  *  .0  *. 1... M.T  2£ LUL 15 0. « _ 

IF  (U2EI IXl.Lt.O.JJl)  GO  TO  1300 

UU«  IU2E2  LIXL-.UNll.lXi  ) /UiEIJXJ - 

UK*8#0* (UU*DTN*C  CPR2EI IXI-PR2EI IX-i) )/(OEN2F( 1XI*U2E{ IXI *P 1 1 ♦ 

1  m?Fi  IXiaU2El LXsl  >  >/pi  U 
GO  TO  1301 
10a  UUiO. 

UK«0. 

101  MT2EJI  IXliMTKJ  l  IX  L-DTN*U2EI  1X1  »tMT2E  ( lXI.-MT2£tlX-l  li/(H»PJ  l>aUML» 
lMa»OTN»IU2EllXI**«l-ini/tDEN2FUX»»*OI-U+B)*MT2E(lX)»lOTN«IIH+2>* 
3  1U3F  ULX1  ~U7E  tlx— 11.1X  t  H*PI  1+U2EI  LX  1*1  DE142-F 1  1 X 1  — DEN2F  I  I  X- 1 1 1  /  IH*Plf 
40EN2F IIX ) 1 I+UU+C0EN2F2I IXI-DENF1 I IX) I/0EN2FI IXI) ♦  (  i*B1*M*2E( I XI*UK 
.  lfJHT2E2  I IX.)  »lL*Ot  1HL2E2JJXUS1U 
31  MTX2I IXUIMT2EI IX»*MT2E2(  lXII/2. 

OO  12  1  X  »MM  *  MW 

61  MTK2,  IX»*MTK2<J)-MT<2<  JJ*I JX-Ul*PI/(XN2-( J-l)*Pl  I 


MUE»MU*SGRT 1 1 T0»TEMN2  < 1 » +TEMM (I , 1 ) ) / 1 2#T0 1 » 

IF (HT<2 ( 1 1 »LT.O. |  WR1TEI6*217)  HTK21 1 1 
ANTK II) -MTK2 1 1 1 **BUM4 
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.LEI  L.GT.MM)  RIll-Il  I )  =  2*TQ»T£.*LN2  (  l  >  /ITFKN2  1 1  HlOtTEMMi  l  » 1 J  1 
OENF2  (  1  )*RT.VF  I  1  )  #l)h  NN2  (  I  )  /  (  1+1.  T  AD0*DLNN2  {  I  1  *  l  KU’F  t  1) -1)  I 
1£JUENFZ(  I )  .LL.O.  )...  v.'Rl  TE(6»  77  1  DENF2  I.IJ  .RT.VFllJ  ,DENN2  tl). 

REMT  ( I  )=OUV.5*DENF2U  )*UN2<  I  )*AMT<<  ll/MUL 

FR I  F  =  2*  A/  (RE-  iT  (  I  >  **t3  ) 

GQ_1 0—4.32 

401  FRI F=0. 

402  112 LU  =aUM2*CEUE2.(  Ll  *UK2.LU-*£RJ  F. 

1FU.GT.MM)  Hi  ( I  - 

HlstHU  1  UH2I  1  1,1/2.- 

C  COMPUT  AT  I (  1  CF  hEAT  LOSS  IN  Tlf'E  DELTA 
_1E__U-.EC  .2  l-GC.  JQ_ LujGOl- 
1F  (I .EQ.MMJ  GO  TO  1001 
..I  E-  LL  . JIT..  MM.1_ G.0— LU- 1.01.2.. 

01=2.*3.14l59*RIN*PI*TUBEL*Hl  (  1  )*<T0*TEKN1<  1  >~TE.XF(  I  ,1 ))*DT 
■Gq... 10- 1003 

1000  Gl=2.#3.l4l59*RIN*I1.5*PI*TU&EL)*Hl(2)*<TO*TEMNl(2J-TEM:!»2,n  ) *DT 
-GO  -1-0-1003 

1001  Q1»2.*3.14159*R1N*TUBFL*<0.5*P1+(XN1-J*PI M*H1(  I )#(T0*TEMN1(I)- 
UEH'lLIULl.lJtqi 

GO  TO  1003 

...1.002-01. =g^ 

1003  QT»QT+01 

IF—LL.GL.y.M ! . GO  TO  4J1U. 

CALL  HTW( I  1 
_ 400  FONT  1  NlIF 

c  reinitialisation 

T  1  MF.l.ml  I MEUQT-N 
KT»KT  +  1 
XNUXN2- 
UN31«UN32 

POMRl »PPMH7 

DENN81-DENNB2 

VESB1,«.VFSQ2— 

TEMNB1*TEMNB2 
00  QA  I X » 1 « M 
UN1( IX)*UN2( IX) 

PRN1  (  IX  )  -'PRN2  t  IX 1 
DENN1 ( IX ) *DENN2 ( IX  I 

VFS1  LLXlalt£&2 1 1X1 _ 

MTX1( IX)«MT<2( IX) 

96  ,  T.F.MN 1 1  I X  >  »T.EMN2 1 IX 1 
DO  53  !«2,M 
53  HU1HH2.U2 

IF  ( XN1-1 • )  5,7.7 

7  UR 1TFIG.2131  TIME*  .PR PPM  .  AVPRN .UNO  LtXU. 

PRINT  214 
DO  91  U»b. K.^5. 

X»(IX-1)*PI 

QHllU(l.«tfFSll  IXUDENNS+II-I/FSUIXI  ItPENNlUXl 
IF  (IX.EQ.l)  GO  TO  88 

*RITF  (4.216)  H.uMl  li^l.PRNimn.l/FSlUXl.DENNIUXl.DMlMXL. 
1TEMNK  IXt.AKtRU/  I ,REMT < lx ) tHl < IX) .TEMMI IX »1> 

GO  TO  91 

8£  WR1TE(6»215)X,UN1UX) ,PRN1 ( IX I »VFS1 ( IXI *0ENN1 ( IX ) ,DMl ( IX) » 

1 TEMN1 1 IX 1 
91  CONTINUE 

DMBl.VFSBl»DLNNS4M-UFSB3)»DENUai 

WRITE  1 6,21 5)  XNl,UNBl,PRN3i.VFSBl.DENN01.DMBl,TEMNBl 

TGI2»KT )*TEMN1 C 13 )*T0 

HP(2«KT)>Hltl3) 

T1MEPUT)«TIME1*T1MEC 
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-  T  Gi  1  »K  I)-T  E  j.;m  1L?  6  )  *  1 0- 
HP!  1.KT  )=H1I26) 

HP  t  3  *  <T  )=Hl(oM 
T  G.  i  A*.  K-LU-T-Ly-’lL  t  .1.A  2_l  *Xa 
HP(4*KT )=Hl ( 142  ) 

JJO-  2a9-  JD-=L.A_ 

DO  710  J<=1*KT 

jftJiJ  .1 1L- 1 3-1 tiR-t  J  K 1 »  LG.L  UD-i  J  k  ).jiTill£PJ  JK  L 

710  COM Ti HUE 
7 OP  c om  T l mu  K 

T'VX  1  (  KT  )  =ThKH (16*1) 

X4X  2-UCT .U  LEMMLZ6. .  Lj_ 

7 WX  3 ( K  T )  =  T  EMM ( 51 »1 ) 
iia4.xja.i-=-i£K:-LULu  iiii 
wR I TE ( 6  *70 )  KT 
m  FORMAT  MS -i. 

END 


SUBROUTINE  COEF ( PRF *DENG , VFS . B *C *D ,G »H.E ) 

I-BETA  *  RTAV _ _ 

COMMON  /MOPwPC/RSP .ETADO 
n=r 


ENS=PRE*( l-TbETA)/D£NNS 

ENG=PRb  » (RSP-ET  AOO.»D£MG  )  /  (  DENG*  (.BSP—  1 )  I 

l- ( 1-VFS ) *( 1-ETADC*DCNG )-( RSP-1 ) *VFS*TBETA 

B=RSP»PRE/2 

C= <RSP*PRE*(DENNS-DENG)-( RSP-1 )  *DENNS*DENG* ( EN6-ENS-R T AV ) )/<DEN5* 
-1ZJ- 

G=DENG/t ( 1-VFS )*B) 

H  ~  I RSP.-.  l.l-»D.£flNSADlN-(L»JL£Xl^iLtOidX  I.AY  LL (. . 

E= ( RSP-1 ) /Z 
E£IURfi 
END 


SUBROUTINE  BP(PR,DEN,TEM,VSM> 

COMMON  /MBPWP/DTN.PI .PRN1I201 ) ,UN1(201 ) .DENN1 ( 20 1 ) . VFS1 1  201 }  * DSV . 


COMMON  /MBP/PRN131  iDENNBl  *UNB1  »UNB2  *XN1  *XN2  *  VFSB1 
,-OMMON  /MRP  B/ XT 

COMMON  /MHBPWP/T0*TEMN1(201»*TEMM< 201*51 ) 


PRZ^RHBI 
DEN2*DENNB1 
VS?  »VFSO 1 

CALL  COEF (PR2 *DEN2 *VS2*B2*C2*D2*G2*H2*E2) 

CALL  COEF_tPRMn»DENNa-l,VFSm.iiBl*-Cai  ,DBl»GBl*Hai  .FBll 

SB1 »VFSB 1* ( 1-VFSb 1 ) /Dfc'NNBl 

T til  «Dill  /DENNB1 

PRX-PRNB1 

D£NX»Q£NiiBl 

UX-UNB1 

ILSXjQlLSlLl- 

DX»DB1 

BX»B01 _ 

CX-CB1 

EX«EB1  _ 

5  DAV«(D2+OX>/2. 

BAVb1B?+BX)/2. 

UAV« (UNB2+UX I /2 • 

IF  tlU.EQ.n  GO  TO  50 
IF  IKT.GT.l)  GO  TO  52 
50  PRX-1. 
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Ut-NX?  L» _ 

UX  =  C  t 
VSX?VFSB1 


52  DXN=IUAV+SuRTtllAV/OAV) )*DTN 
_JU&aXfl2_-UXi'< 

I »XN1/P I 
-j=xxMyp_i — - 
IFI I.EQ.J)  GO  TO  6 


6  PRX=PRN1<  J  +  l  1-MPRNBl-PKNll  J  +  l  )  )*(XXN-J*P1  >/(XNl-J*PI  ) 

.  DLNX=i)LNMl  (  J-UJ.+.IIiE.'WB  L-Q£.VilXJ+  1 )  j  *iXX\-  J  * P  Uxm: LrU *  P 1 1 
UX=UN  1  (  J+l  )  +  (  UNol-UN  1 !  J  +  l )  )*(XXN-J*PI  )  /  !XM-J*Pl  ) 
_VLSX=_Vr_S.li  J+L)  t-LViLSllLri! F S.LLJLtl J  J  *(XXWrJ.*Pi \J  lXMr_J#PJj 
JJ=J+1 

HldsaUjAlJ. _ 

TB=TCMN1< J+l) 
lialsXEMM  (  JJ  f  1 ) 

IF  IRtNTIJJl.LE.O.)  GO  TO  30 
■SX-LJl=£l£M£l.(.JJ-)i  ♦[■UAll.JJ-)  ».*2-«jyjLP t  ■  T  LJLJJ .»*£ ) 

GO  TO  8 
Q.  .SKINsa*. 

GO  TO  8 

7  PRX=PRNi(J+.L)-+(PRNl  LJ.+2  .LrPRMl  J.J.+.l.J.i-t-lXXN^-JXP.l JVPJ 
D£NX  =  0ENN 1 ( J+ 1 ) + ( JENN 1 1 J+2 ) -OENNl ( J+l )  >*fXXN- J*PI  )/P| 

jJX=jJJfmj+XliLl  UN.1.UA2  L-UfLU  J  +.1  J. )  JUXXLLi  J* P-Li  J2A . . 

VSX=VFS1(J+1)+(VFS1 (J+2 l-VFSl (J+l ) )*<XXN-J*PI ) /P I 

■  ).)  =  J  +  ? _ 

H I a=H I  I J+2) 

TflaTEM.Nl  t  J+2) 

TW«TEMM< JJtl) 


SKI  N=DENF1 I JJ ) *  I  UNI ( J J ) **2 • ) / I  RENT ( J J ) **U ) 
ftn  Tn  a 
40  SKINaO. 


DAV* ! D2+DX  )  /2 


UAV*  (UNt)2+UX)/2 


COE AV* I  CX+C2 ) * ( VSX+VFSd 1 ) *0SV/4 


1 ( T0*TB-TW ) +0TN*BUM6# ( E  AV*UAV-SORT ( BAV/OAV ) )*SKIN 

DC0E1«HA1*D$V*VFSB1 
-&CQL2«H2»DSV*VFSdl 
PCOEAV* ( DCOE1+OCOE2 ) /2 . 

-D£MM£W»DEKA>&t*7»A\4txaa^t^^PfiNBl.U'-atN*OCO£AV 

S2«VS2*I1-V$2)/DEN2 

■T2«D2/0ei2 - 

SAV«ISBl+S2)/2. 

TA\/«O.K\/»IT81»\/FSHl+T2»\/F.5An/7. _ 

VFSNEW«VFS81+SAV*(0ENNEW-DENNdl>-TAV*DTN 

IE  lVF5NEWjLE..a».l--VFSflEM«Q.« _ 

IFI  I  At3S (PRNEW-PR2  )/PR2  )-0»U001  )10*1Q»15 
LO  IF  1 1 ABS  IP£MNEiti-DEN2.)  /DJLNLL-0 >0001.1 11  tli-i.L5 
11  IFI ABS(VFSNEW-VS2)*0.wU«01 ) 12  *12  » 15 

15  PR2.PRNEM _ 

DEN2»DENNEW 

VS2-VFSNEW 

CAU  COEF(PR2,DEN2»VS2*B2*C2*D2*G2,H2»E2) 

GO  TO  5 


nwgto 


12-PR=PRNEW - 

DtN=DENN£W 

JLM=PR *  (-l_-tIADQ*Otft  J  JLI  OLN  *  l  1_-E  TADO )  ) 

VSM=VFSNEW 

.Kf-IlJRN. 

END 

SUBPOUT  I  ME  WP(P.V2*DW2*TW2*VFSW2) 

COMMON  V.vePUPC'/RSP  .E.TADO 

COMMON  /M3PWP  N  .  P I  ,  PPM  1  (  201  )  .  UNH  2  O.l)  »  DENN 1  ( 20 1 }  .  VF  S  H  20 1 )  »  DSV  # 

COMMON  /MHUP'aP/  T0»  T uMN  1  (  20 1 )  .  T  EMM {201*51) 

-PK  2_=  PKjNJLIXI 
OtN2=D£NNl  ( 1 )  '* 

_Vi>2.= V£.SilIJ  _ 

CALL  COEF(PR2,OEN2*VS2tl32*C2.D2»G2.H2»E2  ) 
f  Al  I  rnFF  (  PB-M  (  1  1  M  )  .  VFS  1  (  1 )  , B  1  . C  1  .01  .0)  «H.l  . E_Ll 

SW1=VFS1 ( 1 )*( 1-VFS1 ( 1 ) ) /DENNl { 1 ) 

-1U1--.U1  VDENMX IX 
PRX=PRN1<1) 

D£A!X=?J)Omi..a  1 
VSX=VFS1(1 ) 

BXslil. 

CX=C1 

EX=E  1 
JLIXHaL^ 

10  DAV=(DX+D2 )/2. 


UAV=UX/2. 

X=DXN/P1 

X==X+J _ 

PRX=PRN1 ( I ) + ( PRN1 ( 1+1 1-PRN1 (!))*( (X+l )-I ) 

m.+  t.PfcftiU,  LL±1  UU  LUJL(-LA±L 

UX  =  UN1  ( l )+(UNl ( 1+1  J-U.Nl  (  I)  )*(  (X  +  l  )-I  ) 


CALL  COEF(PRX,l)ENX.V$X»BX.CX.DX.DUMl  .DUM2.EX) 

BAV=(B2+BX )/2. 

UAV=UX/2. 

EAV=(EX+F2)/2. 

IF  (REMTI2 ) .LE.O. )  GO  TO  40 

GO  TO  15 

15  PRNEW»PRX-DAV»SORT(BAV/DAV>*UX+DTN*COEAV-DTN*EAV*PUM?*HI (2 )# 

GAV* ( G1+G2 ) /2 • 

SKQ EI  =  H1»DSV»VFSH1  ) 

DC0E2«H2*DSV*VFS1 ( 1 ) 

DCMV*  I0ro£3+DCQE21/2jt 

DENNEW*DENN1U  ) +GAV* < PKNEW-PRN 1 ( 1 ) >+DTN*DC0AV 
TW2«0220EN2 


TAV»DSV*(TW1*VF$1 (1 )+TW2*VFSl (1) )/2. 
VFSNEiV«.VFSl  ( llXSAV*(DliNN£V-Q£NILLLll_l-^lV 
IF  (VFSNEW.IE.O.)  VFSNEWO. 

1 F  C ( ABSIPRNEW-PR2 I/PR2 I-O.UOOl )5.5*6 
5  I F ( (ABS(DENNEW-0tN2)/OEN2l-0.U001)7,7.6 
7  IF( AdSCVF$NEW-VS2l-o,v>wo01 18.8.6 
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6  PR2-PKNEW- 
l>tN2“DtfiNEh 
VS2=VFSNEW~ 

CAU.  C0EF|PR2»iJtN2,VS2.B2.C2»D2,G2.H2.E2  ) 
JaO-XOL  1U 
8  PW2=PR«EW 
—D^.'2=D£NSl£iiL 

T‘.v2=Pw2*  ( 1-CT  ADO*D,V2  )  /  (  Dw2»  ( 1  -ET  ADO)  ) 

...  .V.F S~2-= y£  SMLlii 
RETURN 
— END - 
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APPENDIX  VI 


REPETITIVE  FIRE  COMPUTER  PROGRAM  LISTING 

This  program  requires  input,  from  the  single  shot  program,  of  heat 
transfer  coefficient  (HP)  and  gas  temperature  (TG)  as  a  function  of  time 
after  propellant  ignition  (TIMEP)  at  a  particular  x  location.  It  calculates 
barrel  temperatures  at  that  x  location  as  a  function  of  time  for  repeated 
firings  assuming  HP  and  TG  are  the  same  for  each  cycle.  A  cyclic  time  of 
0.060  second  is  taken.  All  units  are  in  the  metric  system  (meter,  kilo¬ 
grams,  seconds,  kilocalorie,  °JCelvin) 

!*  DIMENSION  A9A(16>/3. 19422*6. 3105*9.4409* 12. 576,15. 7136,18,85191,21 

l|9iQ8lj25.JL3Q0*2A .  26952*31.4091 7*34. 54ft95*37.6Ba8.4n.w;>ft73,43.Qfi,ft7 

21*47.10875,50.24931/  ~  * - 

\*  DIMENSION  AAP< 5*301) 

>♦  COMMON  /ZL/  AAM,AAN",C0n3UC,H0*0EI-R 

!*  -DIMENSION  AK500Q) ,  A2(5000)  *A3H0*3nl)  .A4(10j 

\*  DIMENSION  aAN<20.301) 

'*  ,  READ  300*  RINiR0UT,DELR*0IFFU»CONDUC»Hn.KTMAX 

I*  30o  FORMAT (6c12.6»I5) 

■*  TAm3=300. 

:*  nn=TRout-rin>/oeuuo.oooooi 

l*  L=NN+1 

M*  1  DO  2  1=1, L 


READ(3)  NP(J),TS(J)*TXMEP(J) 

3  CONTINUE 
JKrjK*l 

IF  (JK.EQ.  2)  GO  TO  4 
98  IPT=1 
IjM£l=0t 
<S=1 

KR=1 

?-! 

KP=l 

9  IF(<T.GE,KTMAX)  SO  TO  100 
IF  (AD. fa. lT "SO  TO  41 
IF  (KR.NE,  10*0)  GO  TO  202 
4l  Ai(KO)sTSMMd) 

A2(<D)sTIME1 

SDskDm 

202  HI=(HP(KT)»HP(KT»1) )/2- 
DT=TlMEP(<T4-lT-ffMEP(KT) 

00n  TIM1SQ. 

0TfrsTDELR**2.)/(2.*DXFFU*(HX  *DEUR/C0NnUCn.*DElP/(2.*RlN) ) ) 
DTT2s<D£LR»»2.)/ ( 2 . *DI FFU* ( MO*DEU./CONDUCel • -OELR/ <  ? . *R0UT) ) > 
BtTsDTTl 

IFQTT>DTT2)5*S»6 
6  “D"Tt=3tT2 

S  lF(DTT»(DT-TIMim5*lS.16 
l6  6^T=DT-TIM1 

IS  CHNTS2.«DIFFU»PTT/(QELR»»2,) 

TNE  *4>S(i.-CHNT*lHl«.DEl.R/C0NDuCM.*0El.R/<2.*RjN) ) )*T£MMf  *)♦ 
__ICHNT»HIQELR»T»(KT)  /CQNDUC'»CHNT*(1«»DELR/(2.>btN>  2) 

X*2 

lo  TNEF t K ) s ( 1-CMNT ) *TEMM <  K ) ♦CHNy* C 0. 5*DELR/ (4.*(RiN^fK«l ) *DEL* 1)1* 
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1* 

SUBROUTINE  RUM(3. Ab J) 

2* 

commOn/rl/  aai 

_.3*_ 

a3=2?, /(7.*180, ) 

4* 

AAe=.7979/S3RT(3> 

_5*._ 

AAps1,/(8*3) 

6* 

AAS=Sl N ( <57 . ?96* ( B) -45 ) * AB ) 

7* 

AAh=SIN( <57.296*<3)*45)*A9) 

6* 

AAlsAAE* ( AAH*AaF*AaS) 

_5*. 

AA JzAA-* ( AAS-AAF*AAH) 

10* 

RETURN 

JJJL 

TNI 

_____ 

SUBROUTINE  CUMCC.2T 

2* 

COMMON  /RL/A 

3* 

*G=0. 

4* 

J=l 

5* 

3=1, 

6* 

E=l« 

7* 

1  F={ (C/2, )**(2.*J>7i7-E)*( 1-1  < 

8* 

IF  iA9S(F)  .LT.  0,0001)  SO  T< 

9* 

J=J*1 

10* 

0=0*< J**2) 

11* 

PI+HT7J) 

12* 

G=3+F 

13* 

30~fO  1 

14* 

?  Z=(2.*7./22.)*<A*(,57721>AL0( 

15* 

RETURN 

16* 

JEND 

. 1* . 

SUBROUTINE  MUM<C> 

2* 

COMMON  /RL/3 

3* 

J=1 

4* 

0=1. 

'  5* 

5=1. 

_6* 

LZ-=X< C /2 . )  *  *  <  2  *  JM  *  ( < -1 ,  )  *  *  J ) /I 

7* 

IF  <A3S<F).LT.  0,00001)  SO  TO 

9* 

0=0*1 J**2) 

JUL*. 

S=S*F 

11* 

sWT 

J12* 

2>  RETURN 

13* 

EN? 

__1* 
2* 
_  3* 
4* 
$♦ 
6* 
7* 
8* 
9* 
To* 
n* 
12* 
13* 
14* 
15* 


SU|ROuTlNEjrAI(TEMM,*,AAC»A*OJ 

COviMO'vr/'Vu/AAK 


COMMON  /RL/  A Z_ 

“DIMENSION  fEMM <  30 1 ) » A A«  <  30 1 ) 
3°  10  J=1 *301 


IF  (8  .l-F,  6)  SO  TO  20 

"crurwcTTrAYT 

_ SO  TO  30 

21)  CALL  NU^ld) 

CAu.  CuM(d.AY) 

5?  =Az/AlC^Ar7*iO 

10  continue 
return 


ENO 
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50* 

SI* 

52* 

53* 

54* 

55* 

55* 

57* 

5B* 

59* 

6Q1 

61* 

A2* 

63* 

J&i* 

65* 

M* 

67* 

68* 

69* 

10* 

71* 

72* 

73* 

1M 

75* 

16* 

77* 

7fl* 

79* 

so* 

81* 

82* 

83* 

Ml 

85* 

Ml 

87* 

AM 

89* 

SO! 

91* 

!i* 

93* 

Ml 

95* 

Ml 

97* 

Ml 

99* 

-iflfll 

181* 

ana 

103* 

_kan 

•m* 

106* 

107* 

01311 

109* 

O101 

111* 

112* 

113* 

AIM 

115* 


lTfcMVit  «+1)*CHNT*(0.5-0£LR/(4.*(RIN+U-.1)*DE.UK)))*T£MM(  K-l) 
A3C3=T^£W<K-1)»TWEW(K) 

I*  {A3CD  ,LT.  0.1)  60  TO  8 
<=K»l 

"TFuiGT.NN)  60  TO  80 
SO  TO  10 

So  TNEW<K)=(l.-CHMT*(HO*DELR/COMOuC*l»-a£LR/(2.*ROUT3))*TEMHn?T+CHNT* 
1H0»3ELR*TAv>3/CQ'sI3UC»CHNT»  C 1 , -3ELR/  ( 2 .  *R0t )T ) >  *TEMM (K-t  3 
8  TIM2=TIM1+DTT 
TIM1=TIM2 

IF  (TMEW(H  ,$T.  TEMMtl))  60  to  912 
KP=<P+1 

J  PJL*!L  *iii«_3L_so_lo  912 

apkr=tew<i3 


A2<<3)sTIME1*TIM2*OTT 

K0s<3»1 

IF  < IPT  ,NE.  10*KE>  SO  TO  912 
00  913  UsiiL 
A3{<t*L<)=TEMM(L<> 

913  CONTlNjJl 

A4(<E)sTI4E1*TI'12-3TT 


912  30  7  vi=l»K 

7  TEM4(M) sTNPWtR) 

1003  IF(TI^2.UT.DT)  SO  TO  5 

KTs<T+l 

<Rs<ft+l 

TlM£lsTIM£l»OT 
60  TO  9 

910  WRITE  (6*90)  APKR 
TSAr=TSUTMAX) 

DO  996  LKsi'L 

Riiv)  =  (W-ll*3ELR±PIN. 

TE4'mx)=TE^(L<3-300, 

30. 9|Sl-M=UJS 
AAp(«Y.U)=0. 


TI4£2s0.06*IPT-TIME1 
AAs=3IFFU/( (0.0153**23 
DO  221  <2=1*16 
TlM£3=TlMg2/3. 
TI^sTXME! 

AA4sQ. 

9AsA3A(KZ) 


CAU  PADA(R.BA,TERM*<Z) 
AArs(9A**2)*AAS 


aatsaam*exp(-aar*ti*m:3) 


AAp(<Y.J)SAAPUY»J>*AANUZ.J3*AAT 
997  eOMTlNilE 

TmE3=TIR£3*TlME4 


221  CONTINUE 

00  914  <YS1»3 

IF  < AAP<<Y,301)  *6T»  0.3  80  TO  lit 
9P5F- -AAPtF-Y*313JLl 
80  TO  113 

99  9P<S0. 

3  Ai(<3)sAAP(KY*13*300,*8PK 
A2(<3)=TmEl»KY*TmE4 
K0s<0*l  jjo 


VfrJv.  if 


116* 
117JL 
118* 
_1 1 9* 
12C* 
-1 21* 
122* 
123* 
124* 
_125*. 
126* 
127* 
128* 
.129* 
130* 
.1.31* 
132* 
133* 
134* 
135* 
136* 
137* 
138* 
139* 
140* 
141* 
142* 
143* 
144* 
145* 
146* 
J47*. 
148* 
149* 
150* 
151* 
TSS* 
153* 
"154* 
155* 
156* 
157* 


914  CONTINUE 
9o._PORMATjlH  .llEtl.5) 

'JO  998  U=iVl 

998  TEMM ( L< )  =AA?i  3_»_US*300 .  +3P< 

TlMEl=TIMEl+TIM£2 
90  TO  11 

901  JTTl  =  OEL«**2.)/(2,*JIFFU*<HI  *0eLR/CONDUC+1,+DElR/(2.*RIN) ) ) 

JH2£iD£.LS**2,)/(2,*Ji'FFU*(H0*,0ELR/CON''UC  +  l,-JELR/(?.*4OUT) ) ) 
JTT=0TT1*,9 
If (JTT-JTT2 )  25 * 25, 26 

26  DtT=3TT2*.9 

2 5.  CH.NT=2.*PIfFU«PTT/<3ELp«*?.1 

TNEA(l)  =  a#-CHMT*{HI*3E!-R/C0'^L)C  +  l,fDELR/(2.*RlN))  )*TE'1M(  l>* 

i^M!Hl*JEl«R*TSAF _ /C0MJUC+CHNT*(1  .*DELR/(2.*RlN)  )*TEMM(  2) 

2 1 0-JNSaL* K )  =  ( 1  “CHWT MTp v [(  JO  +CHNT* ( 0 ,5+PELR/ ( 4 . *  <  Rl N+ ( K-l ) *0EL« )  J ) * 

1 T e.MM <  K  +  l )  *CHNT*<n, 5-DELR7T 4 .  *  ( R I N+T<-1 ) * D£LRT)T*> eW(  ~kLTJ - 

A3C3=TNEW(K-1)-.TN£W(<) 

IF  (A3CD  ,lT •  .011  GO  TO  28' 

_K=jO-l _ 

IP(<»ST.NN)  SO  TO  280 
SO  TO  210 

28o  TNEW<K)=(l.-CHNT*tHO*JELR/C0N0uC*l,-OELR/(2T*ftbUT)))*tPMM{K)*CHNT* 
lH0*^ELR*TAM9/CQNJUC4-CHNT»(l.»JFLR/(9.»4flim  )*TEMM(K-1 1  ~ 

28  DO  2?  vi=l,K 

27  TEMM(M)=TNEW(M) 

TImeIstimITi+dtt 

<P=<0»1 

99,  irTTTMEl-IPT*O,06)  90i»il»ll 

ll  IPTsIPT^l _ 

WRIT£  (6.90)  TIM£lfTEMM(l) 

WRITE  (6,90)  TEMM<301) 

tTl=lTIME(lT2»lTjr 

IP  (1T2  .Gt.  1.40S+6)  SO  TO  994 

TPnrrPT.LE.io6>  so  Tirgi 

99a  WRITE  (11)  IPT.TIWEl,Kn»KE»<R 

'WHITE  TTTT  Ul(JTTA?TTjVJ51,KD) 

WRITE  (11)  ((A3(I«J)>Jsl,301)fI=l»ln) 

WHirt  (IT)  ( A4(I)iI=l,iO) 

WRITE  (11)  (TEMM(I) ,l:i,301) 

Wffi  ft  "Til)  (  (~AT\j(X»  J)  r  J=l»301 ) ,  I  si » 15 ) 

9999  EN; 


_1* 

2* 

-34 

4* 

_S* 

6* 

JS* 

8* 

_a* 

10* 

11* 

12* 

13* 

14* 

15* 


SUBROUTINE  PAPA (TEMM, A. SIGMA ■<*) 

COMMON  /  Yl/AAK 

-COMMON. /2L/  AAM,aAN»C0NP*H0»3ElR 

DIMENSION  TENM(30l)»AA<(301)iA«N(20rTbl) ,SIGMA{301) 

iAA=-(a.9?9l/S3RT(**A> 

AA3=3,/(8*2*») 

AA=57.?96»3nA 

A3=22,/(7. v!80. ) 


-MC=i?PN0j»A/!LC15)*AAA*lSINX(AJA-45)*A9)+AA8*SIN((AA*45)*AB))*H0*(- 
IAAA)  *(SIN(  ( AA*45)"*A3)  *(AA3/3, )  *SlN(  (AA-45)*A6)  , 

_AjJ3=i.C0NP*A/,Q.15).*AAA*(SIN((AA-l35J*A3)*AAB*SlN((AA-45)*A8))*HO*(- 
1 A A A ) * ( s I N ( ( A  A-4  5 ) * A3 ) * ( -A  A3/3  >  *SI N ( (*AV»5 ) *ABf) 


CALL  TAI (TEMM, A, AAC, AAO) 
GO  TO  11 
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16* 

12* 

16* 

AS* 

2C* 

21* 

22* 


5  DO  1C  J=l»301 
B-A*TEMM( J1/.015 
AAE=.7979/S0RT<3> 
AAFs1./(B*3) 

AAg=SINU57.296*<3>-45>*AB> 
AA-u5Im(  (57.296  *  (_3J±4  5  >  *  A  8 ) 
AAI=AAE*(AAH*AAF*AAS) 


23*  AAj=AAE*(AAS-AAF*AAH) 

j»4*  AAkj[J)=AAI/AAC-AAJ/AAO 

25*  lC  CONTINUE 

2M>  1 1  AAl=(T£MM(  301 >  * *  2 /2jJL*  ( A  A  K  ( 30l  )**2  )*(!+(  (HO*.Ol5)/(CON3»A)  )  »*2)»C 

27*  ITEMMi 1) **2/2. )*I AA<(1)**2) 

2.6*  po.  3P-2=lLagA 

25*  AAN(K2f J)=AA<(J)/5aBT<AAL) 

3a*  _3fL.COM!  Amue. 

31*  RETURN 

3g*  END 


1* 

11 

3* 

.4* 

5* 

61 

7* 

2* 

9* 

.IQ* 

11* 

JL21 

13* 

14* 


SUBROUTINE  DADA(T£mm,SISMA,K?) 

COMMON  /ZL/  AAM.AAN^QnD. ROGEIR 

DIMENSION  TEMM(301)  tSlG'MAf301)  ,A'An(20,301)  rAAA(30l) 
P0_l0_j=1.30! 

AAA(J)=TEMM(J>*AAN(K2»J)*SIGMA(jy 
In  CONTINUE 

AAm=( 1./2. ) *( AAA( 1)+AAA(301) ) 

J0.2P.J=2>.iaQ. 

AAm=AAm+AAA( J) 

20_C.PNIINJJE 
AAm=DELR*AAM 
90  FORMAT ( H  « 2E10.4) 

RETURN 

END 
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